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ABSTRACT 

(Distribution  Limitation  Statement  B) 


Flexible  and  rigid  pavement  test  sections  were  constructed  and  tested  to  obtain 
data  to  validate  present  criteria,  to  establish  modifications  to  present  criteria, 
or  to  develop  new  criteria  for  the  evaluation  and  design  of  airfield  pavements 
to  be  subjected  to  multiple-wheel  heavy  gear  loads  (MWHGL) .  The  basic  CBR  design 
method  was  expanded  and  modified  to  obtain  a  method  of  design  for  flexible  pave¬ 
ments  subject  to  traffic  by  MWHGL.  The  recommended  method  reflects  a  reduction 
of  existing  U.  S.  Army  Corps  of  Engineers  (CE)  thickness  requirements  that  is 
especially  significant  for  multiple-wheel  assemblies  in  the  higher  operational 
level.  Current  CE  evaluation  and  design  methods  for  rigid  pavements  are  based 
on  stress  in  the  concrete  pavement  as  calculated  from  the  Westergaard  analysis; 
extrapolations  to  the  existing  criteria  were  found  to  be  valid  for  MWHGL  assem¬ 
blies  insofar  as  pavement  thicknesses  were  concerned.  Results  of  traffic  test¬ 
ing  of  the  rigid  pavement  indicate  that  current  jointing  recommendations  allowing 
keyed  construction  joints  may  be  unconservative  when  MWHGL  assemblies  will  be 
trafficking  a  pavement  resting  on  a  low-strength  subgrade. 
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CONVERSION  FACTORS,  BRITISH  TO  METRIC  UNITS  OF  MEASURFMENT 


British  units  of  measurement  used  in  this  report  can  be  e  rted  to  metric 
uni-s  os  follows: 


inches 

feet 

square  inches 
square  feet 
pounds 

pounds  per  square  inch 
pounds  per  cubic  inch 
pounds  t>er  cubic  foot 
inch-pounds 
miles  per  hour 


_ s* _ 

:>.5k 
0.3048 
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0.011521 
1.609344 


_ To  Obtain _ 
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kilograms 
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grams  per  cubic  centimeter 
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SECTION  I 
INTRODUCTION 


1*  PURPOSE 

The  purpose  of  this  investigation  was  to  validate  present  criteria,  to 
establish  modifications  to  present  criteria,  or  to  develop  new  criteria  for 
the  evaluation  and  design  of  both  flexible  and  rigid  airfield  pavements  to  be 
subjected  to  multiple -wheel  heavy  gear  loads  (MWHGL).  MWHGL  as  used  in  this 
report  denotes  aircraft  with  gross  loads  exceeding  600  kips.1 

2.  BACKGROUND 

The  multiple -T*heel  gears  cf  large  new  aircraft  (such  as  the  C-5A  and 
Boeing  7^7)  may  impose  loads  o r*  pavements  that  are  radically  different  from 
those  previously  encountered.  Extensions  to  the  existing  criteria  for  pave¬ 
ment  evaluation  and  design  are  necessary  to  evaluate  the  effects  of  these 
loads.  Data  are  also  required  to  determine  the  relative  destructive  effects 
of  new  and  proposed  aircraft  on  pavement  performance.  The  Army,  Air  Force, 
and  Federal  Aviation  Administration  (FAA)  jointly  sponsored  this  investiga¬ 
tion.  The  investigation  was  conducted  by  the  U.  S.  Army  Engineer  Waterways 
Experiment  Station  (WES),  Vicksburg,  Miss.  Overall  supervision  of  the  tests 
and  all  details  pertaining  to  the  flexible  pavement  portion  of  the  tests  were 
provided  by  WES.  The  rigid  pavement  testing  was  directed  by  the  U.  S.  Array 
Construction  Engineering  Research  Laboratory,  Champaign,  Ill. 

3.  SCOPE 

This  investigation  was  accomplished  by  constructing  and  testing  a  specially 
designed  test  section  consisting  of  both  flexible  and  rigid  pavements  as  de¬ 
scribed  herein.  Testing  consisted  of  instrumentation  measurements  of  deflec¬ 
tion,  strain,  and  stress  resulting  from  applied  static  and  dynamic  (slowly 
moving)  loads;  nondestructive  vibratory  testing  to  determine  wave  velocity  and 
stiffhess;  and  traffic  testing  with  multiple-  and  single-wheel  gear  assemblies. 
This  report  presents  a  sunsnary  of  the  investigation  of  various  theoretical  con¬ 
cepts  of  pavement  behavior  with  advantages  and  limitations  of  each  concept, 
fallowing  a  discussion  of  all  these  concepts,  a  review  of  the  development  of 


table  of  factors  for  converting  British  units  of  measurement  to  metric  units 
is  presented  on  page  vii. 
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the  current  U.  S.  Arm*  Corpi  of  Engineers  (CE)  methods  of  design  for  rigid 
flexible  pav*  is  presented.  Tnen,  the  analyst  of  data  from  the 

test  section  with  resulting  criteria  is  given.  Data  used  in  th*  analysis  and 
criteria  development  were  selected  from  Volumes  II  and  in  of  this  series. 

The  MWHGL  study  represented  such  an  extensive  effort  that  the  report  of 
the  study  was  divided  into  the  following  volumes 1 

I  -  Basic  Report  (background,  summary  of  entire  study,  conclusions, 
and  recommendations) 

II  -  Design,  Construction,  and  Behavior  Under  Traffic 

III  -  Presentation  and  Initial  Analysis  of  Stress-Strain-Deflection  and 
Vibratory  Measurements 

A.  Instrumentation 

B.  Data  and  Analysis 

IV  -  Analysis  of  Behavior  Under  Traffic 
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SECTION  II 

INVESTIGATION  OF  CONCEPTS  OF  PAVMENT  BEHAVIOR 

1.  FLEXIBLE  PAViyENT 

One  purpose  of  the  MWHGL  tests  at  the  WES  was  to  determine  the  behavior  of 
conventionally  constructed  flexible  pavements  subjected  to  heavy  multiple- 
wheel  aircraft  loadings.  The  study  was  to  provide  design  and  evaluation 
criteria  for  airfield  pave  ^ents  for  very  large  aircraft. 

The  analysis  was  to  provide  the  best  flexible  pavement  design  and  evalua¬ 
tion  criteria  presently  possible.  Within  the  past  several  years,  many  differ¬ 
ent  theories  and  computer -oriented  concepts  have  been  under  study  for  use  in 
flexible  pavement  design.  Reviews  of  the  most  prominent  of  these  for  possible 
application  to  flexible  pavement  design  are  presented  herein  in  synopsis  form. 
Coincidentally,  the  existing  semiempirical  concepts  were  investigated  and  re¬ 
viewed  for  the  purpose  of  possible  design  usage. 

Any  pavement  behavior  system  must  involve  a  behavior  system  and  a  set  of 
parameters  to  cause  the  behavior  to  be  similar  to  that  of  the  prototype  being 
represented.  Since  these  two  components  of  a  system  act  as  independent  vari¬ 
ables,  two  different  systems  can  yield  identical  numerical  answers.  There¬ 
fore,  any  philosophically  correct  behavior  concept  can,  tlirough  suitable  ad¬ 
justment  of  these  two  components,  be  made  to  produce  the  desired  output.  The 
problem  area  in  design  formulation  revolves  around  tue  validity  of  the  be¬ 
havior  concept  when  translated  to  a  different  system  of  different  inputs  and 
unknown  output.  The  problem  of  translation  exists  in  a  solely  theoretical  or 
empirical  system  (or  some  combination  of  both)  as  applied  to  the  prediction  of 
the  complicated  structural  behavior  of  flexible  pavements. 

Recent  concepts,  theories,  and  computer  programs  were  examined  in  detail 
to  determine  their  applicability  to  the  MWHGL  flexible  pavement  test  sections. 
All  applications  were  made  using  data  from  the  MWHGL  static-load  instrumenta¬ 
tion  tests.  Each  program  in  its  present  form  was  studied  to  evaluate  the 
soundness  of  theory,  co  rrectness  of  mathematics  and  boundary  conditions,  JJid 
compatibility  of  computer  programs  with  the  WES  computer  system.  The  follow¬ 
ing  is  a  brief  discussioi  of  this  work. 

a.  Layered  System  -  Shell  Laboratories 


A  computer  program  for  design  was  developed  by  the  Shell  Laboratory  to 
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calculate  stresses,  strains,  and  displacements  tor  a  multilayered,  three- 
dimensional,  elastic,  homogeneous,  isotropic  media  (reference  l).  The  basic 
theory  was  derived  from  Boussinesq’s  and  Burmister's  concepts  on  one-  and 
two-layer  elastic  systems  (reference  2).  A  perfectly  rough  interface  condi¬ 
tion  was  assumed  to  exist  between  layers,  with  the  lower  layer  providing  con¬ 
tinuous  support  to  the  layer  above  it.  The  basic  study  has  been  expanded 
from  this  background  to  include  any  number  of  layers  and  loads.  For  adapta¬ 
tion  for  use  in  the  NWHGL  analysis,  changes  were  made  in  certain  phases  of 
the  program.  The  maximum  number  of  loads  was  increased  from  10  to  30,  and 
the  number  of  terms  in  the  Bessel  functions  was  altered  to  afford  the  accu¬ 
racy  desired.  The  program  computed  all  orthogonal  values  of  stresses, 
strains,  and  deflections,  as  well  as  principal  values  and  directions  of  total 
stresses  and  strains,  strain  energy,  and  strain  energy  of  distortion.  The 
preliminary  investigation  indicated  an  applicability  of  a  program  for  layered 
systems  under  normal  surface  loads  to  the  flexible  pavement  portion  of  this 
test. 

Initial  computer  output  from  this  program  reflected  good  agreement 
between  two  computed  values  of  deflection  for  a  multiple -wheel  gear  in  com¬ 
parison  with  measured  data.  A  typical  input  statement  and  format  are  shown 
in  figure  1,  and  a  typical  output  statement  is  given  in  figure  2.  A  compara¬ 
tive  plot  of  computed  and  measured  (from  Volume  III)  values  of  vertical  stress 
versus  depth  for  15,000-  and  30,000-lb  static  loadings  is  shewn  in  figure  3- 
This  computer  program  is  relatively  large,  but  can  be  reduced  greatly  by  re¬ 
ducing  the  number  of  layers  or  less  greatly  by  decreasing  the  number  of  loads. 
Trial  computations  were  made  with  various  material  properties  as  input  values 
in  an  attempt  to  duplicate  the  measured  values.  Although  this  could  be  ac¬ 
complished  for  one  set  of  loading  conditions,  the  results  could  not  then  be 
translated  to  other  loading  conditions  for  duplication  of  field  results. 

Also,  the  method  used  two  values  of  limiting  strain,  one  at  the  bottom  of  the 
surface  course  and  one  at  the  top  of  the  subgrade,  and  accounts  for  the  ef¬ 
fects  of  traffic  repetitions  on  the  pavement  structure  through  selection  of 
suitable  strain  values. 

The  Shell  multilayer  program  results  were  compared  with  instrumented 
readings  of  stress,  deflection,  and  strain.  In  item  4,  the  computed  surface 
elastic  deflection  for  center-line  12-wheel  30,000-lh  loading  was  0.226  in. , 
and  the  equivalent  instrumentation  measurement  wes  0.272  in. ;  at  the  6-in. 
depth,  the  computed  deflection  was  0.226  in.,  and  the  equivalent  value 
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Figure  1.  Typical  Computer  Input  Statement  and  Format 
for  Shell  Multilayer  Program 
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interpolated  from  instrumentation  measurements  was  0.266  in.  These  computed 
values  were  based  on  16  loaded  areas  (12  load  wheels  plus  4  load  cart  out¬ 
rigger  wheels)  and  7  layers.  The  Young's  modulus  of  elasticity  E  values 
were  determined  for  each  layer  by  multiplying  the  measured  CBR  values  by 
1500  to  arrive  at  E  values  in  pounds  per  square  inch.  This  is  an  empirical 
relationship  developed  from  dynamic  triaxial  test  results  correlated  with  CBR 
data.  It  should  be  noted  that  E  ,  when  measured  by  a  vibratory  method,  is 
higher  than  secant  or  tangent  E  obtained  from  other  testing  procedures. 
Poisson's  ratio  v  was  assumed  to  be  0.3  for  the  cohesionlesu  materials  and 
0.4  for  the  cohesive  materials.  The  values  of  E  determined  from  the  vibra¬ 
tory  testing  were  not  used  because  they  were  not  available  when  this  program 
was  run.  In  most  instances,  the  vibratory  values  of  E  ,  later  determined, 
were  higher  than  the  above-mentioned  values. 

b.  layered  System  -  Chevron  Oil  Company 

A  computer  program  developed  Chevron  Oil  Company  (reference  3)  was 
based  on  an  elastic,  N-layered,  two-dimensional,  homogeneous,  isotropic  con¬ 
cept  similar  to  the  Shell  program.  Deflection  computations  agreed  closely 
with  results  of  the  Shell  program.  Only  limited  work  beyond  initial  investi¬ 
gation  was  accomplished,  primarily  because  of  the  constraint  of  a  one-load, 
two-dimensional  input.  This  would  require  an  equivalent  load  to  replace  the 
multiple -wheel  loading  configurations.  No  applicability  to  flexible  test 
sections  existed  for  multiple -wheel  configurations.  However,  this  program 
was  used  for  single-wheel  computations.  This  program  seems  feasible  for  the 
determination  of  stress,  strain,  or  deflection  on  flexible  pavements  subjected 
to  single-wheel  loads.  The  program  is  relatively  short. 

c.  Energy  Method 

The  energy  method  of  analysis,  as  applied  to  beams  and  bars,  was 
studied  to  determine  its  applicability  to  this  project.  Energy  is  the  abil¬ 
ity  to  do  work  and  is  measured  by  the  amount  of  work  a  body  is  capable  of 
doing  in  changing  from  some  given  condition  to  some  standard  condition.  In  a 
defined  system,  the  external  work  can  be  equal  to  the  change  in  internal 
energy  of  the  pavement  structure,  or  the  sum  of  energies  applied  to  the  system 
is  equal  to  the  energy  used  plus  the  energy  transferred  outside,  if 
equilibrium-stability  exists.  This  internal  energy  can  he  computed  using  the 
principle  of  minimum  potential  energy  and  related  to^  the  external  work  done 
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by  appliea  loadings.  This  principle  of  minimum  potential  energy  states  that 
if  a  system  is  in  equilibrium  both  internally  and  externally  and  if  a  very 
small  change  in  displacement  is  imposed,  then  the  external  work  done  by  the 
loads  must  equal  the  change  in  strain  energy  6U  ,  where  6  symbolizes  a 
small  change  and  U  is  strain  energy.  The  strain  energy  (area  under  force - 
displacement  plot)  versus  depth  was  computed  (figure  4)  by  the  Shell  multi¬ 
layer  program  and  used  in  conjunction  with  the  energy  method. 

The  principle  is  stated  in  terms  of  displacements  or  strains  and  ap¬ 
plied  loads  on  deformable  bodies  (references  4-7).  The  external  work  W 
equals 

|  £UT  ^  U  =  U  (1) 

where 

T 

U  -  U  transposed 

Kj,  =  stiffness  matrix  derived  from  energy  methods 

.  Also 


du 


(2) 


where 

=  reactions 
u  =  displacements 

2 

Furthermore,  the  stiffness  matrix  K  is  defined  as  K<  .  =  a u/(du.  du.)  , 
which  expresses  a  force  induced  by  a  unit  displacement.  For  the  linear  case, 
this  matrix  is  composed  of  constants. 

A  nonlinear  relationship  may  be  considered  and  the  relationship  be¬ 
tween  force  F  and  change  of  length  x  is  as  follows,  for  exanple: 

F  =  kx  +  ax 
F  =  kx  -  ax' 


where 

k  and  a  =  constants 

When  a  =  0  ,  this  matrix  would  be  composed  of  first-order  values. 


9 


inilllsiMiiisiiisHilifiHiiiiitsiiiiiiHj 

iissiin8sis;gs^ii^^yis>is4i^isia:§i 


iHimssniHu 


lismiiRsu 


feliiiiimisisitifiiimstftititimiijfsiiiiiiiiiiiiiifiiit 

ilSLIlIiSSiiiiBHHiSiUltSSHmSMSSHiSIlSIHSiSmilHI 


a  *  *  i  *  *  »  •  *£5P 

88Basi8i8i'898»^aaaftasa8888aas89ia8  8au8<lt8i»8a<Mifi»tiiatitaf8aM84t 
S118»II9il^SlM9M81V88USI8lffi88l9lf«lU8MIHHUIII«llKUCittlUt1 

mnnnnnMnmnmmmunnmmunmnmuuummJi 

*,-  »§*•  <  v  «*  .  m  -wYjtiiB  *-*«:  &  *  ■*  *  *  b-«  ii--*  *  i»mc  iu>  i  ri«4ii»» 


.••  •  ,  •  •  •  ?  .  S  -  - •  '  »  »  -  --  ^  •«  ■’  *  •  ■!  «»  •  f  =■*  *  «  \ 

iiiiHtiihmiMiiiitimmmiimmitii 

'•i.isn'.r  »a  ■  ■  •  -  r-*«  **;»*?*  t »»  £ 

*5 smisisnin  <'**CkftJkk.»s*»«i 


RM  <M  }  ct(>  t  *  *>M  t  ^  .  *WI  »  ■;  «t«  «  [»..  !MI»  tl 

-?***Ea*‘Sea  ft  5«**'k>a$;&a  k  3  »*ft 

>  .  *  v*  *  a*  *:**  '  »»  *  ■  s  «  t  *  *  «  a  a  r  *  *  *  F 

fiiiiiiiiimiRiliiwiifinmtniayfiuiUiHmiiiUUim 

iisiiiiuiiiiHiiiiiiiiniiiimiiiiitiiHUuuiiiiiiiHintttt 


ArVL-TR-70-113 


I  i 

i 

i 

Figure  U.  Typical  Plot  of  Strain-Knorgy  Vorouo  Depth 
Item  Ut  15,000-lb  Single-Wheel  Static  Loetf 


Other  terms  A  /  in  +he  system  axe  the  work  i  tems  resulting  from 
applied  loads  or  potential  energy  of  applied  loads.  A  typical  strain- 
energy  versus  depth  plot  for  item  4,  MMrL  test  section,  15 ►COO -lb  single- 
wheel  load  is  shown  in  figure  4.  The  concept  of  ■virtual  work,  i.e.? 


A(U  +  y)  »  a(tJ  -  W)  -  0  (4) 

was  reviewed  and  considered  applicaole,  but  no  adaption  was  made  to  this  test 
section  because  of  time  constraints  and  limitations  imposed  in  developing  the 
nonlinear  constitutive  equations.  Another  area  of  difficulty  revolves  around 
ti  e  fact  that  the  input  energy  level  p^r  stress  repetition  can  be  quite  dif¬ 
ferent  from  that  in  seemingly  similar  situations.  The  application  of  this 
method  to  the  analysis  was  questionable  based  on  the  type  of  support  assumed 
and  on  the  lack  of  sufficient  data  to  describe  the  stiffness  matrix  of  the 
structure,  bat  in  principle  the  method  seems  applicable  for  analysis.  No 
computer  program  was  available  for  this  method  at  the  time  of  the  study. 

d.  Stiffness  Method  -  Structural  Analysis  (References  8-10) 

This  method  ds  incorporated  into  the  finite  difference  and  finite 
element  methods  as  applied  to  pavement  structures  ar.d  is  shown  in  detail  for 
basic  explanation. ' 

A  plane  truss  c*  frame  may  b*  considered  be  be  composed  of  members 
(elements)  connected  at  the  joints  (nodes).  The  force  components  and  moments, 
such  9 a  Xi  .  •  lA  ,  M.  ,  acting  at  fill  the  nodes  constitute  the  generalized 
force  vector  X  .  The  deflection  components  and  notations,  such  as  u^  . 
vi  5  ^i  5  ^  all  ^ne  constitute  ^he  displacement  column  vector  |u^  . 

Therefore 

•Vs*.  .  '  . r  - . 

{xj  -  {vVVVVV  (5) 

HvK'VVvvVm..«1.v»1} 

The  forces  and  displacements  are  related  through  functions  depending 
upon  the  geometric  and  material  properties  of  the  elements  of  the  structure. 

In  the  stiffness  method  of  matrix  structural,  analysis,  the  relationship  be¬ 
tween  forces  and  displacements  et  the  nodes  mey  he  expressed  as: 
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M  -  [k]  {„} 


(7) 


In  {x}  -  [k]  {u}  ,  K  is  shown  as  the  stiffness  matrix.  The  ele¬ 

ment  .  of  the  stiffness  matrix  is  the  force  aloi;g  i  corresponding  to 
unit  displacement  along  j  .  The  concept  of  the  stiffness  matrix  may  be  dis¬ 
cussed  in  relation  to  springs. 


X  »  ku  ,  where  k  is  the  stiff¬ 
ness  of  the  spring 


'f1  U4 

(l)  Stiffness  matrix  for  an  element. 


xi  »  ui  > 


©  © 


X.  ,  u. 
J  3 


X. 

1 


-X . 


©Ui®  ® 


H=  £  ■  W  -  ft 


Consider  a  spring  IJ  with  J  fixed.  Force  at  J  due  to  u^  alone, 


,  and  X..  =  -X.,  =  -ku.  .  Due  to  u.  alone,  (u.  =  0)  , 

i  l  i 


(Uj  —  0}  ,  X±  —  ku^  ,  wu.NA.  ^ 

X.  =  ku.  ,  and  X.  =  -X.=  -ku.  . 

0  0  i  0  3/  *  t  \ 

The  stiffness  relation  between  <X  >  and  <u[  may  be  written 

KL  f k  -k 

(Xj)  -k  k 

The  stiffness  matrix  for  the  spring  element  ic,  therefore 

k  -k 


as 


a 


-k 


(8} 


(9) 


where 

k  =  the  stiffness  constant  of  the  spring 

(2)  Stiffness  matrix  for  an  assemblage.  For  an  assemblage  of  two 
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springs  kith  stiffness  constants  k&  and  k^ 


©■^©vvv 

xx  — —  — - 

®„©  ©  <D 

i  i  _  v 


Due  to  u^  alone, 

*1  =  Vl  ’  *2  =  *Vl  *  X3  =  0 


©  dT©'  ® 

X,  -*m — *W*w*^^ X-. 

si 

0  ©  a  ©  b  (D  ©' 

- - #*-  X^ 


Due  to  Ug  alone, 

x  .  ,t  „  h’*A*  Va 

1  *  1  *3  ■  ‘Vs 

Due  to  u^  alone, 

X1  =  0  ’  X2  -  -V3  ’  X3  -  V3 


Stiffness  relation  for  tho  spring  assemblage  is 


hi  ka  -ka  0  \\ 

*>}=  -ka  <ka  +  V  -\  <u2 

h)  0  -“b  (U3 


The  stiffness  matrix  for  the  assemblage  may  itself  be  assembled 
from  the  stiffness  matrices  for  the  individual  elements,  expanding  each  of 
them  to  conform  to  the  order  and  composition  of  the  total  stiffness  matrix. 
Thus,  for  the  spring  a  ,  the  element  stiffness  matrix  is 


Hip  - 


k  -k 
a  a 


0  0 


Similarly,  the  stiffness  matrix  for  the  spring  b  ,  is 


(3X3) 


h> 

h 


(3X3) 
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Superposing  the  two  expanded  matrices  gives  the  stiffness  matrix  of  the 
assemblage: 


M  =  [k]12  +  [4 


-ka  <ka  +  V  -S, 
0  0 


as  before. 


For  a  long  chain  of  elements,  the  superposition  of  the  element 


stiffhess  matrices  leads  to  a  banded  matrix  as  indicated  below. 


M- 


The  stiffness  matrix  it*  also 
symmetric  with  , 

in  accordance  with  the  law 
of  reciprocal  deflections 
(Maxwell’s  law). 


(Deflection  at  i  due  to  a  unit  force  at  j  is  equal  to  the  deflection  at  j 
due  to  a  unit  force  at  i  .) 

In  the  stiffness  method,  the  nodal  displacements  due  to  a  set  of 
specified  forces  are  generally  the  primary  unknowns.  But  X  =  Ku  for  the  en¬ 
tire  assemblage  cannot  be  inveited  at  this  stage  because  the  total  K  matrix 
is  singular.  Physically,  the  situation  corresponds  to  rigid  body  motion, 
which  must  first  be  prevented  by  prescribing  support  (boundary)  constraints  * 
before  a  solution  is  possible. 

The  displacements  u  are  separated  into  the  prescribed  set  u^ 
and  the  unknown  set  .  The  force  components  may  be  correspondingly 
divided;  X  will  then  be  the  reactions  at  the  prescribed  displacements,  and 
Xp  will  be  the  specified  load  components  along  the  unknown  displacements  u^  , 
The  partitioning  of  the  K  matrix  along  the  same  lines  leads  to  the 
following: 


=  _ ; _ 


This  represents  the  two  matrix  equations,  and 


twr-  *  ?S 


(15) 
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M-  ("«]  fu4*  [k“b]{u4 

M'  [S.]{"4*  [lt»e]{J4 

Known  Known  Unknown 

From  the  latter,  the  unknown  displacements  |  ug  |  **or  the  prescribed  dis¬ 
placements  ju^  and  forces  w  may  be  determined.  Once  M  is  known, 
the  unknown  forces  at  constraints  can  b».  calculated.  (The  solution 

technique  described  is  only  symbolic,  as  the  prescribed  displacements  |uoc| 
may  be  introduced  into  the  matrix  equation  X  =  Ku  without  partitioning.) 

The  stiffness  method  of  matrix  analysis  corresponds  to  the  clas¬ 
sical  method  of  slope  deflection  for  (indeterminate)  beams  and  frames  in  which 
the  moment  at  one  end  of  a  member  is  expressed  as  a  function  of  the  rotations 
at  both  ends  and  the  relative  deflection.  The  stiffness  method  is  also  known 
as  the  displacement  method.  (The  other  important  method  of  matrix  structu  :al 
analysis  is  the  flexibility  method  or  force  method.) 

e.  Finite  difference  Method 

The  finite  difference  method  is  the  representation  of  the  derivative 
in  discrete  form  (reference  'Ll).  In  general.,  the  n^*1  difference  is  defined 
as 


An  An-1  An-1 

Ayk  = A  yk+i  - A  yk 


(16) 


where 

£ 

A  =  n  difference 

yfc  *  the  mate  of  a  finite  set  of  arguments  x^ 

This  method  is  developed  from  the  definition  of  the  derivative  as  the  limit, 
as  distance  approaches  zero,  of  the  difference  between  two  values  divided  by 
the  distance  between  them.  For  points  spaced  relatively  close  together,  a 
continuous  function  can  be  represented  by  values  taken  at  discrete  points. 

In  this  method,  it  is  assumed  that  the  behavior  at  each  grid  point  represents 
the  behavior  in  the  region  surrounding  the  point.  Region  differential  equa¬ 
tions  represent  the  behavior  at  a  point  in  the  field.  By  representing  the 
derivatives  as  finite  differences  and  substituting  these  expressions  in  the 
field  equations,  a  system  of  differential  equations  is  converted  to  a  system 
of  linear  difference  equations. 

The  nature  of  finite  difference  schemes  requires  that  four  properties 
of  the  solution  be  investigated.  These  four  properties  are  convergence, 
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consistency,  stability,  and  error.  Hie  scheme  must  be  so  formulated  that  the 
error  growth  with  each  step  is  tolerable.  This  method,  similar  to  finite 
elements  analysis,  produces  solutions  with  discretization  and  rounding  error. 

A  proper  balance  must  be  made  between  the  use  of  large  numbers  of  points  or 
elements  with  corresponding  high  roundoff  error,  and  small  numbers  of  points 
or  elements  and  a  correspondingly  high  discretization  error. 

The  finite  difference  method  is  a  type  of  numerical  technique  that 
deals  with  a  means  of  numerical  approximation  of  field  equations.  The  similar 
ity  of  the  technique  to  finite  element  analysis  precluded  any  further  investi¬ 
gation  of  this  method. 

f.  Finite  Element  Method 

The  term  "finite  element  method  of  analysis"  refers  to  a  method  of 
solving  boundary  value  problems  in  which  the  distribution  of  stress  and  strain 
throughout  a  solid  mass  under  the  influence  of  either  internal  or  external 
loadings  is  desired  (reference  12).  The  method  generally  has  been  applied 
only  to  two -dimensional  masses:  in  soil  mechanics,  plane  stress  or  plane 
strain.  (Three-dimensional  solutions  are  not  generally  available  for  use.) 

In  two-dimensional  problems,  the  mass  is  divided  into  elements  composed  of 
plane  figures,  such  as  triangles.  The  elements  are  connected  to  each  other 
at  the  vertexes  or  nodes,  thus  forming  an  interconnected  assemblage  of  finite 
elements  equivalent  in  shape  and  area  to  the  mass.  Although  the  elements  are 
connected  only  at  the  node  points,  an  additional  restriction  is  placed  on  the 
displacement  field  which  requires  that  boundaries  between  adjacent  elements 
are  in  continuous  contact  during  deformation. 

The  initial  assumption  made  in  the  finite  element  method  is  the 
character  of  i  he  displacement  field.  In  a  computer  program  used  at  WES,  it 
is  assumed  that  the  displacements  are  linear  functions  of  the  coordinate 
system.  The  assumption  that  the  displacement  field  is  linear  will  lead  to  a 
unique  strain  associated  with  each  triangular  element,  i.e.,  a  straight  line 
in  the  undeformed  element  remains  straight  after  deformation.  The  displace¬ 
ments  within  each  element  can  be  evaluated  in  terms  of  the  local  coordinates 
of  each  element;  thus,  if  the  element  were  to  approach  an  infinitesimal 
size,  the  displacements  in  th  assemblage  of  discrete  elements  would  approach 
the  idealization  of  the  continuum. 

The  nodal  displacements  are  found  by  a  procedure  often  used  in  struc¬ 
tural  mechanics.  That  Is,  the  force  on  a  member  is  equal  to  the  displacement 
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times  the  stiffness  of  the  element.  The  stiffness  is  the  force  causing  a 
unit  displacement  of  the  element,  and  the  stiffness  matrix  for  the  entire 
assemblage  of  elements  is  formed  by  relating  the  position  of  the  nodal  points 
and  the  stress-strain  relations  of  the  elements  under  investigation.  Briefly, 
displacements  are  described  in  terms  of  the  local  coordinates  of  each  element, 
and  the  strain  within  each  element  is  computed  from  the  displacements  of  the 
node  points  and  the  size  of  the  elements.  Stress  is  expressed  as  a  function 
of  strain,  and  the  force  on  the  nodes  of  the  element  is  a  function  of  stress. 
Ihus,  by  proper  mathematical  manipulation,  an  equation  of  the  type  F  =  kd 
can  be  written  for  each  element  in  which  the  stiffness  k  is  a  function  of 
the  element  coordinates  and  the  stress-strain  relation,  and  the  displacement 
d  can  be  resolved  into  its  x  and  y  components.  Since  each  (triangular) 
element  has  three  node  points  and  two  reactions,  the  equation  F  kd  can 
be  thought  of  as  representing  six  linear  equations.  Also,  since  each  ele¬ 
ment  is  connected  to  its  neighbor  at  the  node  points,  the  solutions  for  the 
nodal  displacement  in  one  element  must  be  compatible  with  the  solutions  for 
its  neighbors.  Thus,  a  total  stiffness  matrix  K  can  be  formed  by  adding 
all  of  the  individual  stiffnesses  k  ,  resulting  in  a  system  of  simultaneous 
linear  equation-  of  size  2n  ,  where  n  is  the  number  of  nodal  points.  A 
solution  to  these  equations  is  obtained  by  specifying  either  forces  or  dis¬ 
placements  at  the  boundary  and  solving  for  the  displacements  at  each  node. 

As  mentioned  previously,  when  the  nodal  displacements  are  known,  the  average 
strain  and  stress  in  each  element  can  be  computed.  The  basic  equation  for 
finite  element  analysis  relates  virtual  work  to  virtual  strain  energy  as 
follows  (references  8  and  13): 

/vOF  [»]{«„}  av  +  /vWTHdv  +  /sWT{p}ds  +  H 

=  (fv  M T  HH  A{*\  u?) 

or 

k  +  v  vx}=  DFH 


/,  M 1 M  {*.} 

/v  [*]’{«}  *  • 


dv  ,  equivalent  forces  for  initial  strains 
equivalent  forces  for  body  force  effect 
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M- 

W- 

M- 


equivalent  forces  for  surface  traction 


vector  of  actual  concentrated  nodal  forces 


/» [■] 1  [»]  w  * 


stiffness  matrix  for  the  element 


vector  of  nodal  displacements 

function  of  coordinates  matrix 

function  of  stress-strain  relation  matrix 

elasticity  matrix 

initial  strain  matrix 

body  forces  matrix 

surface  traction  matrix 


Other  element  shapes  may  be  used. 

A  typical  flexible  pavement  structure  was  analyzed  using  quadra- 
lateral  finite  elements  (references  l4  and  15),  and  the  results  were  compared 
with  Burmister's  solution  for  a  three-layered  problem.  This  comparison  is 
shown  in  figures  5  and  6. 

Computer  programs  exist  for  the  two-dimensional  case  and  can  be 
almost  unlimited  in  size,  depending  on  accuracy  desired  and  the  number  of 
simultaneous  equations.  Each  equation  would  normally  contain  six  coeffi¬ 
cients,  each  of  which  must  be  evaluated  on  the  basis  of  the  geometry  and 
stress-strain  relation  of  the  stressed  mass.  Most  computer  codes  are  derived 
from  an  elastic-stress  analysis  written  by  Professor  E.  L<  Wilson,  University 
of  California  (reference  16). 

The  axisymmetric  finite  element  program,  results  w-:rf  compared '  with . 
the  instrumentation  data.  This  comparison  indicated  .an  duplicate 

pattern  ns  the  layered  system  analysis  for  a  vertical  noma;  stress  versus 
depth  plot .  The  calculated  stress  did  not  attenuate  «s  rapidly  as  -measure¬ 
ments  and  had  a  different  shape.  The  deflection  basins  indicate  a  close 
agreement  in  both  the  finite  element  analysis  and  ife  elastic  half-space 
analysis.  However,  the  curve  shai'es  from  both' of  those  analyses-  differed 
from  the  instrumented  deflection  curves.  v  -  V 

‘The  limitations  of  two-dimensional  loading  and  pavement  response*  the 
requirement  for  measurable  values  of  E  and  v  ,  the  heglsci  of  Stress 
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rer^etitions ,  and  the  difficulty  in  accurately  describing  the  constitutive 
equations  all  limited  the  use  of  the  finite  element  program  in  accomplishing 
the  objective  of  this  study*  The  finite  element  program  results  agree 
closely  with  the  elastic  half-space  analysis  results  as  shown  in  figure  6. 
However,  the  field  instrumentation  readings  do  not  agree  with  either  of  these 
analyses.  This  is  discussed  in  detail  in  Volume  III. 

The  three-dimensional  program  eliminates  the  two-dimensional  loo  ding- 
response  limitation  but  is  still  in  the  development  stage.  The  program  is 
mammoth  in  size  and  the  band-width  increase  in  the  three-dimensional  program 
will  greatly  increase  computational  effort.  The  limitation  of  adequately 
describing  the  construction  material  characteristics  will  still  exist  in  the 
bug -free  three-dimensional  program. 

g*  Mechanistic  Models 

Mechanistic  models  consist  of  mechanical  components,  such  as  springs 
(representing  elasticity),  dashpots  (representing  viscous  behavior),  and 
Coulomb  dampers  (representing  plasticity).  Various  combinations  of  these 
components  have  been  used  to  represent  pavement  structures.  One  such  model 
was  developed  at  Purdue  University  (reference  17).  Basically,  the  theory 
used  for  the  mechanistic  model  in  the  Purdue  study  was  derived  from  a  trans¬ 
lation  of  a  Russian  book  by  V.  Z.  Vlasov  and  N.  M.  Leont’ev  (reference  l8). 

An  infinite  beam  resting  on  an  elastic  foundation  was  used  to  represent  a 
landing  mat-soil  structure. 

A  limited  investigation  of  this  study  indicated  that  a  more  general 
fonn  of  the  solution  was  available  for  use  in  studying  flexible  pavements. 

A  performance  prediction  portion  of  the  Purdue  study  was  developed  for  a 
Vlasov  and  Leont’ev  model  with  an  infinite  beam  and  a  hyperbolic  stress 
distribution.  A  limitation  of  this  program  was  the  restriction  of  two- 
dimensional  analysis  treating  single-  or  multiple -whee3  loads  that  act  only 
on  a  two-dimensional  loading  plane.  An  advantage,  however,  was  the  ability 
to  account  for  a  nonsymmetric  loading  condr  in.  The  first  model  in  this 
study,  which  was  based  on  elastic  theory,  was  shown  to  be  capable  of  dupli¬ 
cating  a  given  loading  system.  It  is  noteworthy  that  the  model  behavior  was 
extremely  sensitive  to  the  magnitude  of  the  subgrade  modulus  k  ,  which  was 
found  to  decrease  with  traffic.  However,  this  decrease  in  k  may  be  in¬ 
fluenced  by  the  effect  of  construction  technique. 

The  second  mechanistic  model  developed  at  Purdue  Was  capable  of 
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simulating  actual  loading  sequences  and  could  also  provide  a  means  of  esti¬ 
mating  residual  deformations  of  a  surface.  This  process  involved  an  itera¬ 
tive  method  of  duplicating  the  loaded-deflection  basins  from  actual  mea¬ 
sured  input  data.  This  model  handles  an  asymmetric,  load  that  can  treat  a 
viscoelastic  model,  which  predicts  a  response  to  a  simulation  of  actual 
sequence  of  traffic  loads.  This  model  is  a  two-layer  system  with  the  sur¬ 
facing  layer  elastic  and  tue  subgrade  layer  viscoelastic.  A  computer  pro¬ 
gram  for  determining  deflection  of  two  layers  was  available  but  was  not 
used  in  conjunction  with  this  project.  The  program  is  relatively  small. 

h.  Viscoelastic  Models 

(1)  General.  Freudenthal  and  Lorsch  (reference  19 )  were  among  the 
first  investigators  to  employ  a  linear  viscoelastic  analysis  to  infinite 
beams.  In  their  study,  the  soil  support  was  replaced  by  a  series  of  Kelvin, 
Maxwell,  or  Standard  linear  solid  elements.  For  each  of  these  elements,  the 
authors  were  able  to  develop  relationships  for  the  deflection  of  an  infinite 
beam  subjected  to  time -invariant  concentrated  and  uniform  loads.  Other 
investigators  have  continued  similar  work,  and  development  has  continued 
through  the  work  of  Barksdale  (reference  20).  Barksdale  developed  one  creep 
compliance  response  for  an  asphalt  mixture  and  another  for  a  clay  subgrade  in 
his  analysis  of  a  viscoelastic -layered  pavement  system.  A  3  imitation  of 
using  viscoelastic  models  is  the  difficulty  in  assigning  representative 
numerical  values  to  the  parameters  for  even  the  simplest  mode. 

(2)  WES  study.  A  viscoelastic,  multilayered,  half-space  system 
under  stationary  and  moving  axisymtiotric  loads  has  been  analyzed  in  another 
study  (reference  21).  Solutions  have  been  computed  for  normal  stress,  radial 
stress,  tangential  stress,  shear  stress,  vertical  deflection,  and  radial 
displacement  at  any  point  within  the  half- space.  This  viscoelastic  solution 
was  based  on  the  elastic -viscoelastic  correspondence  principle  and  was  ob¬ 
tained  by  applying  the  Laplace  transformation  to  replace  the  time  variable 
with  a  transformed  variable.  The  time  dependency  was  changed  to  an  associ¬ 
ated  elastic  problem.  This  program  has  not  been  written  for  computer  usage 
and  therefore  was  not  used  in  conjunction  with  this  project. 

i.  Discontinuous  Qrthotropic  Plates  and  Pavement  Slabs 

The  theory  involved  in  this  University  of  Texas  computer  program 
(reference  22)  primarily  consists  of  a  recursion  solution  of  a  system  of 
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beamu  resting  on  an  elastic  subgrade  represented  by  %  system  of  springs »  The 
model  consists  of  discrete.'  elastic  beams  connected  at  their  nodal  points  by 
elastic  blocks*  Torsion  is  accounted  for  by  r«vls  connecting  midpoints  of 
opposite  beams.  This  model  accounts  for  a  two-layer  system.  One  stipulation 
of  the  model  is  the  requireLenr,  for  a  specified  ratio  of  slab  'width  to  slab 
length  grid  system.  Applicability  cf  this  program  to  flexible  pavement  analy¬ 
sis  is  limited  because  of  the  restriction  to  two  layer**,  requiring  a  deriva¬ 
tion  of  equivalent  relationship  for  multilayers.  The  plate  analysis  assumes 
that  the  pavement  is  a  flexural  menf  er  extending  in  two  dimensions  'wich  the 
significant  generalized  forces  >f  transverse  force  and  two  moments  at  the 
joint.  The  computer  program  is  available  and  is  relatively  large  in  capacity* 


less  Theory 


In  the  Odemark  equivalent  thickness  theory  (reference  23),  the  upper 
layer  of  thickness  is  replaced  ty  an  equivalent  thickness  of  the  lower  mate¬ 
rial  according  to  a  relaticn  with  thickness  of  layer  S  (modulus  of  elas¬ 
ticity)  and  v  (Poisson’s  ratio).  The  specific  equations  followtheBms- 
sinesq  solut.’ on  and  are  purely  elastic.  *  In  the  Gdemark  method,  the  stresses 
and  displacements  are  evaluated  by  considering  the  paveaesub  layers  to  behave 
as  a  slab  resting  on  the  subgrade  soil.  The  deformation  moduli  are  deter¬ 
mined  by  plate  load  tests.  The  design  criterion  used  is  to  limit  the  max iflwsa 
curvature  of  the  deflected  pavement  surface.  Results  derived  froc  the  equi¬ 
valent  layer  theory  agree  closely  with  those  from  Burmister  (reference  2), 
assuming  the  same  E  and  v  parameters.  Computer  programs  exist  for  the 
application  of  this  theory,  but  the  programs  were  not  used  in  this  project 
due  to  the  questionabler.es s  of  layer  equivalency. 


k.  Elastic  Plate  on  Dense -Liquid  Bottom  Layer 


Application  of  the  Westergaard  theory  (reference  24)  to  flexible  pave¬ 
ment  design  has  been  investigated  only  for  the  case  of  a  vertical  deflection 
immediately  under  the  c*  iter  of  a  uniformly  loaded  circular  area  for  two 
layers*  Both  layers  ware  assumed  to  be  homogeneous  and  isotropic  and  of  infi¬ 
nite  horizontal  extent.  When  applying  the  Westergaard  theory,  it  is  assumed 
that  the  bottom  layer  material  is  a  dense  liquid.  Westerga&rd  analyzed  the 
system  of  an  elastic  plate  on  a  dense-liquid  subgrade  using  mathematical 
series  expansion.  The  usefulness  of  this  theory  is  limited  by  the  basic 
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assumptions  and  is  also  limited  to  two  layers;  therefore,  it  was  eliminated 
from  further  study  in  flexible  pavement  design. 

1 •  Elastic  Plate  on  Elastic  Solid  Subgrade 

flbgg  and  Holl  (references  25  and  26,  respectively)  independently  ana¬ 
lyzed  an  elastic  plate  on  an  elastic  solid  subgrade  by  expressing  the  sub¬ 
grade  reaction  as  a  Fourier-Bessel  transform.  The  solution  is  valid  only 
for  a  symmetrical  leading.  Also,  the  determination  of  stresses,  strains, 
deflections,  and  subgrade  .reactions  may  be  computed  only  under  the  center 
line  of  the  applied  load.  The  solution  of  the  elastic  layered  system  does 
not  have  some  of  the  restrictions  inherent  xn  the  plate  theory.  This  work  is 
limited  to  use  in  a  two-layered  system.  Finally,  essentially  the  same  infor¬ 
mation  obtained  by  this  solution  can  be  obtained  from  Burmister's  solution  cf 
the  elastic  layered  system.  Therefore,  this  theory  was  not  considered 
further. 


m.  Barenberg  Shear  Layer 

Though  it  is  common  practice  to  use  an  elastic  layered  system  analy¬ 
sis  for  pavement  structures  composed  of  granular  materials,  the  application 
of  such  an  analysis  is  obviously  not  valid.  Tensile  stresses  are  developed 
at  the  bottom  of  a  layer  supported  by  a  weaker  layer.  Since  the  E-values 
of  most  granular  materials  are  greater  than  the  E-values  for  subgrades, 
tensile  stresses  will  develop  and  the  equilibrium  assumed  in  setting  up 
basic  equations  for  the  system  will  be  violated.  Barenberg*  s  approach  (refer¬ 
ence  27)  was  to  examine  conditions  for  a  continuously  supported  shear  layer 
that  has  no  flexural  resistance.  The  development  of  a  critical  shear  layer 
may  occur  when  the  resistance  to  shear,  usually  provided  by  the  subgrade,  is 
l.^st.  No  computer  program  was  available  for  the  Barenberg  shear  layer  con¬ 
cept,  bub  the  approach  seems  applicable  to  certain  conditions  which  develop 
in  flexible  pavement  structures. 

n.  Peattie  Method 

The  Peattie  method,  which  is  still  in  the  process  of  development, 
uses  two  design  criteria  (reference  28).  The  vertical  stresses  on  the  sub¬ 
grade,  as  well  as  the  radial  tensile  strain  in  the  surfacing  layer,  are  kept 
within  certain  allowable  limits.  The  deformation  moduli  of  pavement  layers 
are  determined  in  the  field  by  vibrational  techniques.  This  method  is 
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essentially  a  part  of  the  multilayered  elastic  method  and  is  not  a  separate 
procedure . 


The  one-layer  program  that  assumes  a  pavement  is  elastic,  homogeneous, 
isotropic,  weightless,  and  of  semi -infinite  extent  lias  been  programmed  for  the 
WES  computer  system.  This  program,  which  was  derived  from  Boussinesq*  s  work 
(reference  13),  may  be  modified  to  allow  the  deflection  to  reach  zero  at  any 
predetermined  distance  from  the  loaded  area.  The  assumptions  are  circular 
load,  one  layer,  uniform  load  distribution,  and  no  lateral  surface  stresses. 

With  proper  assumption  of  Poisson’s  ratio  and  modulus  of  elasticity, 
the  program  can  closely  duplicate  measured  deflection  data  at  any  given  point. 
Also,  the  change  in  the  program  to  allow  an  offset  cutoff  in  deflection 
values  can  be  used  to  describe  the  equivalent  single-wheel  load  causing  the 
deflection. 

A  plot  of  maximum  elastic  deflection  versus  depth  is  shown  in  fig¬ 
ure  7.  This  plot  shows  a  comparison  of  measured  data  with  predictions  derived 
from  the  single -layer  elastic  theory  program.  The  measured  data  represented 
static-load  test,  30, 000 -lb -per -wheel,  12-wheel  assembly,  item  3*  The  Young’s 
modulus  of  elasticity  E  was  determined  by  inverse  computations  using  the 
measured  deflection  at  the  surface.  The  findings  from  other  studies  (refer¬ 
ence  29)  also  indicate  that  the  stress  distribution  at  normal  temperatures 
in  conventionally  constructed  flexible  pavements  follows  relatively  closely 
the  pattern  predicted  by  the  elastic  theory. 

The  measured  vertical  stress  profile  with  depth  was  different  from 
the  elastic  single-layer  case,  and  an  attempt  was  made  to  duplicate  these 
measurements  by  changing  the  concentration  factor  n  .  The  value  of  the  con¬ 
centration  factor  (reference  30)  parameter  n  depends  on  the  properties  of 
the  soil  mass.  For  n  =  3  ?  in  the  following  relation, 


where 

<3  =  vertical  stress 

Zk 

n  =  concentration  factor 

P  ~  load 

Z  «  depth 

B  -  smallest  angle  with  vertical  ordinate  and  component  of  force  in 

direction  of  radius  vector 
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Figure  7*  Comparison  of  Predicted  Maximum  Elastic  Deflection  with 
Maximum  Elastic  Deflections  Measured  During  Static  Load 
ing  of  Item  3  Using  360,000-Ib  12 -Wheel  Assembly 
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The  vertical  stress  component  is  the  same  as  that  obtained  by  Boussinesq  for 
a  semi- Infinite  body  satisfying  the  assumptions  of  the  theory  of  elasticity. 
Using  this  concept,  the  instrumented  vertical  normal  stress  basin  at  a  dc^  th 
of  2.75  ft  under  a  12 -wheel  static  load  was  duplicated  using  an  n  =  6.5 
value.  The  same  procedure  was  attempted  using  deflection  measurements,  but 
this  could  not  be  accomplished  in  a  simple  manner  because  the  equations  of 
compatibility  (reference  13)  could  not  be  satisfied. 

p.  Discussion 

The  preceding  brief  summaries  represent  a  major  portion  of  the  theo¬ 
retical  treatment  afforded  the  background  analysis  of  this  project.  Other 
methods,  such  as  Laplace  transformations,  Canadian  method,  and  Navy  method, 
were  reviewed  for  possible  usage  but  have  not  been  reported. 

The  previously  discussed  theories  allow  a  better  understanding  of 
the  behavior  of  continuous  soil  media  from  a  simplified  viewpoint.  In  the 
absence  of  knowledge  for  a  specific  condition,  these  theories  all  contribute 
to  a  general  enlightenment  regarding  soil  structure  behavior.  There  are 
many  limitations  and  constraints  imposed  by  each  of  these  concepts,  coupled 
with  advantages.  There  is  no  such  thing  as  an  exact  theory  to  explain  an 
aspect  of  the  physical  world.  All  scientific  knowledge  is  fundamentally  em¬ 
pirical  and  is  based  on  observed  physical,  behavior.  A  theory  may  be  used  to 
explain  or  predict  physical  behavior  within  the  present  limitations  of  know¬ 
ledge  or  ability  to  observe.  A  precise  mathematical  derivation  of  a  theory 
often  deludes  the  observer  into  believing  that  bis  theory  is  exact,  but  this 
is  not  true  as  observational  powers  are  limited  by  precision  of  available 
measuring  equipment  (reference  31 )  •  The  MWHTrL  test  represents  controlled 
experimentation  necessary  to  accomplish  specific  objectives.  The  following 
limitations  would  hinder  a  purely  theoretical  treatment  for  the  MWHGL  project. 

(1)  The  theories  are  based  on  a  static  loading  condition  and  thus 

do  not  account  for  the  detrimental  effect  of  repetitive  design  loadings  or  for 
the  effect  of  mixed  loadings  on  the  performance  of  the  structure. 

(2)  Most  concepts  are  an  extension  of  classical  structural  material 
analysis  in  which  material  properties  are  reasonably  exact  quantities.  For 
this  analysis,  the  determinations  of  two  typical  input  parameters,  Young1  s 
modulus  of  elasticity  E  and  Poisson’s  ratio  v  >  are  at  best  an  estimate 
for  different  soils.  The  determination  of  whether  E  is  in  tension  or  com¬ 
pression  may  not  have  a  significant  effect  on  the  behavior  of  the  pavement 
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layer  under  load,  but  becomes  quite  important  when  evaluating  stresses  in 
layers  based  on  strain  values, 

(3)  Body  forces  (motion,  Inertia)  are  generally  ignored  because  they 
are  higher  order  terms  in  comparison  with  forces  distributed  over  the  surface 
of  the  body  (surface  forces).  Th^v  may  or  may  not  have  an  effect  on  the  pave¬ 
ment  structure,  which  is  known  to  exhibit  certain  degrees  of  deterioration 
even  in  the  absence  of  applied  loading. 

(4)  Layers  must  be  composed  of  isotropic,  homogeneous  materials  for 
which  Hooke's  linear  law  may  be  applied  to  linear  systems.  Also,  the  upper 
layers  must  be  weightless  and  infinite  in  extent  in  the  horizontal  direction, 
but  must  be  of  finite  thickness. 

(5)  Certain  boundary  conditions  must  be  imposed  to  identify  the 
specific  problem.  Type  of  support  and  conditions  at  the  layer  interfaces 
must  be  specified.  Normally,  the  top  layer  must  be  free  of  normal  and  shear¬ 
ing  stresses  outside  the  loaded  area;  in  the  lower  layer,  the  displacement  and 
stress  must  approach  zero  as  the  depth  approaches  infinity.  Similarly,  it  is 
necessary  that  the  layers  be  in  continuous  contact  with  each  other. 

(6)  Results  from  each  program  depend  on  the  constitutive  equations 
that  have  been  developed  for  the  particular  material.  At  present,  there  are 
no  realistic  constitutive  equations  for  the  materials  employed  in  construc¬ 
tion  of  flexible  pavements  or  to  describe  their  property  changes. 

(7)  Variations  of  vertical  subgrade  stress  occur  with  variations  in 
pavement  temperature,  subgrade  moisture,  and  vehicle  speed,  and  as  a  result 
of  deep  frost  action.  The  change  in  the  moisture  content  in  the  subgrade  is 
accompanied  by  a  change  in  the  rheologic  properties  of  the  soil.  No  program 
handles  these  changing  conditions  caused  by  environmental  effects  or  those 
related  to  previous  stress  history. 

(8)  The  time  and  space  locations  of  loadings  are  almost  random  be¬ 
cause  of  the  manner  in  which  an  aircraft  operates.  Also,  the  distribution  of 
load  repetitions  is  difficult  to  define  for  aircraft  prior  to  actual  operation. 
This  random  transfer  of  coordinate  system  complicates  the  assumptions  required 
of  these  programs  and  presents  another  problem  area  in  the  mathematical  ex¬ 
pression  of  this  loading  (for  example,  point  loading  or  uniformly  distributed 
loads), 

(9)  Numerous  arbitrary  failure  criteria  presently  eiaployed  are  diffi¬ 
cult  to  incorporate  into  a  generalized  theory.  Most  theories  seemingly  refuse 
to  operate  a3  a  preselected  failure  criterion  is  approe  ed.  A  limiting 
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quantity  of  stress,  strain,  deflection,  or  other  parameter  is  currently  not 
established  as  an  acceptable  unique  value  of  failure.  The  assumption  that 
the  system  obeys  Eulerian  conditions  regardless  of  the  stress  level  is 
questionable. 

(10)  Any  representation  of  the  soil  structure  as  a  series  of  springs, 
dashpots,  dense  liquids,  etc.,  is  limited  by  this  representation  and  by  the 
assumptier  that  the  system  is  not  altered  by  changes  in  the  system.  The  main 
disadvantage  of  this  type  of  simulation  is  the  difficulty  in  assigning  param¬ 
eter  values  that  correlate  with  various  actual  conditions,  such  as  quality 
of  materials  from  a  viewpoint  other  than  strength,  i.e.,  drainage,  intended 
use  of  structure,  etc. 

The  above  limitations,  coupled  with  the  inability  to  apply  any  one  theoretical 
concept  to  the  total  problem  of  the  design  and  evaluation  of  flexible  pavement 
airfields,  led  to  the  review  of  the  present  Corps  of  Engineers  CBR  method  and 
investigation  to  determine  its  adaptability  to  the  analysis  of  this  test 
section. 

2.  RIGID  PAVEMENT 

A  wj.de  variety  of  algorithms  are  currently  available  to  predict  pavement 
response  to  loads.  These  solution  methods  are  generally  applicable  only  to 
static  or  quasi-static  loads.  The  basic  differences  among  the  solution 
methods  lie  in  the  assumptions  of  subgrade  and  pavement  structure  character¬ 
istics.  The  analysis  methods  usually  applied  to  rigid  pavement  structures 
involve  the  assumptior  of  either  a  dense-liquid  subgrade  or  an  elastic  solid 
subgrade  together  with  an  assumption  of  one  or  several  elastic  layers  (plates) 
over  the  suograde.  Considerable  work  has  been  done  in  computerizing  these 
solutions.  Of  particular  interest  has  been  the  use  of  finite  element  tebh- 
niques  to  more  accurately  describe  the  structure  and  to  provide  for  more 
realistic  boundary  conditions.  The  computerized  solutions  including  the 
finite  element  procedures  are,  nevertheless,  subject  to  the  same  assumptions 
listed  above;  i.e.,  elastic  behavior  is  assumed  for  all  components  and/or 
subgrade  behavior  is  characterized  by  a  dense  liquid,  and  the  loads  are  static. 
The  algorithms  generally  applicable  to  rigid  pavements  are  discussed  in  the 
following  paragraphs.  (Details  already  given  in  the  flexible  pavement  discus¬ 
sions  are  not  repeated  here,.) 
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a.  E^ugtic  Plate  on  Dense -Li quid  Bottom  Layer 

Th*s  method  is  commonly  referred  to  as  the  Westerga&rd  analysis  (refer¬ 
ence  32).  The  pavement  slab  is  assumed  to  be  an  elastic,  homogeneous,  iso¬ 
tropic,  uniform  plate  supported  by  a  dense  liquid  of  infinite  depth.  The  sub¬ 
grade  is  represented  by  a  series  of  independent  linear  springs  that  react  in 
direct  proportion  to  the  slab  deflection  at  the  point.  The  slabs  are  assumed 
to  cover  a  semi-infinite  half- space.  Application  of  the  Westerga&rd  analysis 
to  rigid  pavement  design  is  given  in  Section  IV. 2. 

b.  Elastic  Plate  on  Elastic  Solid  Subgrade 

In  this  method,  the  pavement  slab  is  assumed  to  be  an  elastic,  homoge¬ 
neous,  isotropic,  uniform  plate  supported  by  an  elastic  solid  of  variable 
thickness  (references  25  and  26).  The  elastic  solid  is  characterized  by  a 
rigidity  factor  that  is  a  function  of  the  modulus  of  elasticity  and  Poisson* s 
ratio  of  the  subgrade  material.  One  advantage  of  this  method  is  that  the 
subgrade  materials  are  not  independent,  i.e.,  a  deformation  occurring  at  one 
point  influences  adjacent  points  to  some  degree.  One  disadvantage  of  this 
analysis  method  is  defining  the  input  parameters  for  the  elastic  solid. 

c.  Layered  Systems  (References  1-3) 

In  this  type  of  analysis,  the  pavement  structure  is  assumed  to  be 
composed  of  a  variable  number  of  elastic  layers  resting  on  an  elastic  solid. 

The  analysis  has  the  advantage  of  being  capable  of  solving  problems  of  com¬ 
posite  pavements  and  of  examining  conditions  in  a  wide  variety  of  locations 
in  the  pavement  structure.  One  disadvantage  is  defining  the  pseudo-elastic 
constants  for  the  various  layers  based  on  standard  testing  procedures  in 
which  the  appropriate  boundary  conditions  are  not  reasonably  approximated. 

d.  Discontinuous  Orthotropic  Plates  and  Pavement  Slabs 

This  method  was  developed  by  the  University  of  Texas  and  is  a  finite 
element  approach  to  the  problem  (reference  22).  Slabs  are  represented  by  a 
grid  of  beams  connected  by  elastic  blocks  ~nd  torsional  springs  all  supported 
by  vertical  springs.  The  output  presents  the  deflections,  binding  moments, 
and  twisting  moments  at  all  nodal  points  in  the  slab.  The  output  describes 
conditions  in  the  slab  only;  no  information  is  produced  on  the  subgrade  con¬ 
ditions.  The  solution,  if  programmed  with  sufficient  nodal  points,  tends  to 
yield  the  same  results  as  the  testergaard  method.  Tie  finite  elements 
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representation  allows  more  freedom  in  describing  boundary  conditions,  such  as 
cracks,  discontinuous  subgrade  support,  nonuniform  thickness,  etc. 
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SECTION  III 

REVIEW  OF  CURRENT  CE  METHODS 

1.  DEVELOPMENT  OF  CBR  METHOD  OF  DESIGN  OF  FLEXIBLE  PAVEMENTS  (REFERENCE  33) 

During  1928  and  1929 >  the  California  Division  of  Highways  made  an  inten¬ 
sive  investigation  of  pavement  failures  throughout  California.  Areas  that  had 
failed  were  investigated  to  determine  the  local  drainage  conditions  and  other 
factors  affecting  the  stability  of  the  pavement.  A  test  method  was  sought 
that  would  be  satisfactory  for  establishing  the  density  that  should  be 
achieved  in  the  construction  of  subgrades  as  well  as  the  shear  strength  that 
was  required  to  support  the  anticipated  traffic.  A  quick,  simple  test  was 
sought  that  could  be  used  to  evaluate  the  quality  of  subgrade,  subbase,  and 
base  course  materials.  Static- load  field  tests  were  unsuccessful,  mainly  be¬ 
cause  of  the  influence  of  the  elastic  and  plastic  deformation  of  the  material 
beneath  the  load  and  also  because  saturated  conditions  were  not  considered. 

A  California  bearing  ratio  (CBR)  test  was  devised  that  determined  the  re¬ 
sistance  of  a  specimen  to  lateral  displacement  (reference  3*0  >  thus  measuring 
the  combined  influence  of  cohesion  and  internal  friction.  Samples  were  first 
thoroughly  consolidated  to  the  density  obtainable  unde^  good  construction 
methods  or  to  approximately  the  density  ultimately  produced  by  traffic  in  good 
subgrade  materials.  A  plot  of  total  thickness  versus  CBR  values  resulting 
from  this  study  is  shown  in  figure  8.  Curve  A  was  developed  on  a  basis  of  ex¬ 
perience  available  in  1929  and  was  believed  to  be  adequate  for  average  traffic 
conditions.  Curve  B  represents  values  derived  from  the  original  survey. 

During  the  early  19**0's  the  responsibility  for  the  design  and  the  construc¬ 
tion  of  military  airfields  was  assigned  to  the  Corps  of  Engineers,  U.  S.  Army. 
Later  the  Corps  was  delegated  the  same  function  pertaining  to  civil  airfields 
required  for  military  occupancy.  After  deliberation  and  investigation,  the 
principles  used  by  the  California  Highway  Department,  known  as  the  CBR  method, 
were  adopted  tentatively  for  the  design  of  flexible  pavements.  Among  the  rea¬ 
sons  for  adoption  were:  (a)  The  CBR  method  had  been  correlated  to  the  service 
behavior  of  flexible  pavements  and  construction  methods  and  had  been  used  suc¬ 
cessfully  by  the  State  of  California  for  a  number  of  years.  Two  other  states 
were  known  to  have  methods  of  a  similar  nature  that  had  been  successful;  (b) 

It  could  be  adapted  for  immediate  use  in  airfield  pavement  design;  (c)  It  was 
thought  to  be  as  reasonable  and  as  sound  as  any  of  the  other  methods 
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Figure  8.  Original  California  Pavement  Design  Curves  Relating 
Total  Pavement  Thickness  to  the  CBR  Value  of  the 
Subgrade  (Reference  35) 


AFWL-TR-70-113 

investigated;  (d)  Tne  subgrade  strength  value  could  be  determined  either  in 
the  laboratory  or  in  the  field  with  simple  portable  equipment.  Testing  could 
be  done  on  samples  of  soil  in  the  condition  representative  of  the  foundation 
moisture  state  under  most  pavements. 

To  evaluate  figure  8  curves  in  terms  of  aircraft  traffic,  the  relationship 
among  (a)  relative  weights  of  loaded  vehicles,  (b)  characteristics  of  airplane 
and  automobile  taxes,  and  (c)  lateral  distribution  of  traffic  were  considered. 
On  the  basis  of  these  factors  and  recognizing  the  need  for  verification  tests, 
curve  A,  figure  8,  was  assumed  to  represent  a  12,000-lb  aircraft  wheel  lead, 
and  curve  B,  a  7000- lb  aircraft  wheel  load.  Service  behavior  records  of  ade¬ 
quate  pavements  had  indicated  that  it  was  necessary  for  elastic  deformation  to 
govern  pavement  performance  over  an  extensive  period  of  use.  Accordingly,  the 
Office,  Chief  of  Engineers  (OCE),  decided  to  develop  empirical  curves  by  ex¬ 
trapolating  the  original  data  on  the  basis  of  the  elastic  theory.  Since  all 
bearing  tests  (even  those  such  as  the  CBR  west  conducted  on  confined  specimens) 
are  essentially  shear  tests  and  since  shear  deformation  must  be  eliminated  in  a 
satisfactory  pavement,  shear  stresses  were  used  as  a  guide  in  making  the 
extrapolation. 

A  review  of  airplane  tire  data  indicated  that  a  uniform  tire  pressure  of 
60  psi  could  be  selected  for  the  entire  group  of  planes  in  use.  Wheel  loads 
of  25,000,  40,000,  and  70,000  lb  were  selected  to  cover  the  range  of  existing 
and  future  heavy  aircraft  loads.  Contact  areas  were  computed  from  wheel  loads 
and  tire  pressures.  Circular  areas  were  used  for  ease  of  computation  and  also 
because  the  shear  stresses  in  the  base  course  and  subgrade  did  not  vary  materi¬ 
ally  when  elliptical  areas  were  used.  Shear  stresses  were  computed,  as  shown 
in  figure  9,  by  use  of  stress  tables  (reference  36).  The  thicknesses  of  base 
course  and  pavement  corresponding  to  bearing  ratios  of  3>  5,  7,  and  10  were 
located  on  the  stress  curve  for  the  12,000-lb  load  curve,  and  the  stresses 
corresponding  to  these  thicknesses  were  noted.  On  the  basis  that  these 
stresses  should  not  be  exceeded  for  other  wheel  loads  in  order  to  retain  a 
uniform  standard  of  design,  the  stress  values  were  located  on  the  curves  for 
the  25,000-,  40,000-,  and  70,000-lb  wheel  lceds  (figure  9)«  (The  values  indi¬ 
cated  by  the  curve  in  figure  8  for  the  12,000-lb  wheel  load  were  comparable  to 
the  California  values  for  medium  heavy  highway  traffic.)  The  thicknesses  cor¬ 
responding  to  these  stresses  were  transferred  to  the  graph  of  thickness  versus 
bearing  ratio,  and  the  curves  similar  to  those  shown  in  figure  10  were  drawn. 

The  preliminary  design  curves  obtained  by  extrapolation  were  developed  at 
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Figure  9.  Extrapolation  of  Highway  Pavement  Thicknesses  by  Elastic 
Theory  (Tire  Pressure,  60  psi).  From  Reference  37 


CALIFORNIA  BEARING  RATIO  (CBR),  IN  AT  0.1-INCH  PENETRATION  FOR  COMPACTED  AND  SOAKED  SPECIMEN 
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Figure  10.  Tentative  Design  of  Foundation  for  Flexible  Pavements 
(Reference  37) 
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the  meeting  of  consultants  held  in  February  19^2  in  Washington,  D.  C.,  which 
included  T.  A.  Middlebrooks  and  others  from  OCE;  0.  J.  Porter,  who  developed 
the  California  method  for  design  and  evaluation  of  highway  pavements;  and  Pro¬ 
fessor  Arthur  Casagrande  of  Harvard  University.  The  consultants  each  made  in¬ 
dependent  calculations  to  extrapolate  the  basic  curves*  Those  of  Mr.  Porter 
were  based  on  an  allowable  deformation  for  the  various  loads,  whereas  thc^e  of 
Professor  Casagrande  were  based  on  the  relationships  between  the  relative  size 
of  the  loaded  areas.  The  three  sets  of  computations  (OCE,  Porter,  Casagrande) 
were  in  substantial  agreement  when  the  data  were  reviewed  by  the  group.  It 
was  decided  that  the  average  thicknesses  shown  by  the  three  extrapolations 
were  reasonable  for  the  low  CBR  values;  however,  the  majority  agreed  that  the 
less  conservative  values  should  be  chosen  for  higher  CBR  values.  Figure  10 
expresses  the  best  judgment  of  the  engineers  and  consultants.  This  later  be¬ 
came  capacity  operational  criteria. 

Test  sections  and  evaluation  of  field  performance  continued,  and  in  April 
19^9 >  a  relationship  for  the  effect  of  less-than-capacity  operations  was  of¬ 
fered  (reference  33)*  Also, the  tern  "cycles  required  to  produce  one  coverage," 
200 

defined  as  -=rr-  >  where  N  =  number  of  main  wheels  and  W  =  tire  width  in 
NW 

4nches  ,  was  introduced.  Full-scale,  multi-item  tests  were  conducted  through 
19>5  (references  38-^6). 

In  1958,  an  analysis  of  all  available  se. -vice-behavior  data  from  test  sec¬ 
tions  a. id  prototype  airfields  (reference  39)  indicated  that  the  CBR  design 

criteria  for  single-whee:  ids  could  be  expressed  in  the  two  basic  parameters 

t  CBR  , 

-7=.  and  -  (reference  ’{).  These  parameters  could  be  reduced  to  a  single 

V  A  Pe 

plotted  curve  that  separated  service-behavior  data  with  regard  to  failures  and 
no:  failures,  but  without  regard  to  repetitions.  The  relationship  for  the 
capacity-operational  level  is  shown  in  figure  11  and  is  expressed  mathemat¬ 
ically  as 


(19) 


where 

t  =  design  thickness,  in. 

P  =  single-wheel  load,  lb 
CBR  =  measurement  of  supporting  layer  strength 
A  =  measured  tire  contact  area,  sq  in. 

In  1959*  the  effects  of  less-than-capacity  design  repetitions  and 
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multiwheel  configurations  were  combined  into  a  single  relation  that  Is  essen¬ 
tially  a  combination  of  all  developed  CBR  design  criteria  resulting  from  CE 
research  (reference  48) «  The  equation  for  relating  flexible  pavement  design 
in  the  lower  CBR  regions  (approximately  <15  CBR)  is 


t.  =  f 


ESWL 

IcbrT 


A 

TT 


(20) 


where 

t  =  design  thickness,  in. 

f  =  percent  of  design  tiiickness  (0.23  log  C  +  0.15) 

ESWL  =  equivalent  single-wheel  load  (load  on  a  single  tire  that  produces 
the  same  vertical  deflection  on  '  ’  :  supporting  medium  as  that  par¬ 
ticular  multiple-wheel  assembly  -with  the  same  single-wheel  tire 
contact  area) 

C  =  coverage  (sufficient  wheel  passes  to  cover  every  point  of  a  traffic 
lane  once) 

CBR  =  measure  of  supporting  layer* s  strength 

A  =  measured  contact;  area  of  loaded  tire,  sq  in. 

This  formulation  represents  the  required  thickness  to  prevent  excessive  de- 

t* 

flection  of  underlying  pavement  layers  resulting  from  repetitions  of  load. 

Construction  control  procedures  are  thoroughly  treated  in  appropriate  CE 
specifications  and  are  a  part  of  the  design  that  accompanies  these  thickness 
requirements.  These  controls  are  intended  to  ensure  that  a  pavement  is  not 
subject  to  differential  settlement  of  materials  in  the  underlying  layers  or 
lateral  displacement  of  the  subgrade  resulting  from  water  absorption.  The  in¬ 
fluence  of  mixed  aircraft  traffic,  climatic  effects,  and  aircraft  lift  and 
speed,  etc.,  are  also  automatically  incorporated  in  this  method.  The  judgment 
of  the  design  engineer  as  to  intended  use  of  facility  and  distribution  of 
traffic,  channelized-nonchnnnelized  traffic,  etc.,  are  also  incorporated  in 
the  CE  design  concept. 

This  is  the  current  CE  procedure  for  the  design  of  each  layer  in  a  flexi¬ 
ble  pavement  structure <  Concurrently,  each  layer  must  be  compacted  suffi¬ 
ciently  to  prevent  detrimental  consolidation,  as  shown  in  reference  49,  and  the 
wear-  and  weather-resistant  surface  must  be  designed  >  as  not  to  be  displaced 
by  traffic  (reference  50). 

This  brief  background  investigation  indicated  that  the  CBR  method  was  both 
reasonable  and  sound  and  had  been  successfully  correlated  with  service  behav¬ 
ior.  It  was  decided  that  the  CBR  method  offered  the  best  available  means  of 
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providing  design  and  evaluation  criteria  and  would  teat  reflect  the  rinding  cf 
the  Mk'HGL  test  section.  This  method  also  represents  tne  best  experiment  at  ion 
and  semierrpirical  technique  of  analyziiig  pavements  with  consideration  of  theo¬ 
retical,  materials  testing,  and  construction  central  aspects.  The  next  phase 
involved  adapting  the  present  CBR  method  of  behavior  prediction  to  the  spe¬ 
cific  results  of  traffic  from  the  test  section. 

2.  DEVELOPMENT  OF  CURRENT  RIGID  PAVEMENT  DESIGN  PROCEDURE 

a.  Plain  Concrete 

The  development  of  a  design  procedure  for  any  structure  to  resist  load 
requires  that  the  most  likely  mode  of  failure  be  determined  for  the  weakest 
nonredundant  element  comprising  the  structure.  The  mode  of  failure  assumed  by 
the  current  CE  design  method  for  a  rigid  pavement  subjected  to  repeated  roll¬ 
ing  loads  is  that  a  crack  is  initiated  at  the  bottom  of  the  slab  due  to  flex¬ 
ural  action  and  migrates  upward  through  the  slab  under  continued  traffic.  Re¬ 
peated  application  of  traffic  results  in  further  cracking  of  the  slab  in  the 
same  manner  until  the  slab  is  cracked  into  many  pieces.  This  action  continues 
until  the  integrity  of  the  slab  is  destroyed  to  the  point  where  it  cannot  ef¬ 
fectively  distribute  load  through  beam  action  and  excessive  permanent  deforma¬ 
tions  occur  because  of  excessive  shear  forces  in  the  foundation  materials.  As 
the  slab  deteriorates,  deflections  due  to  load  increase,  causing  a  working  of 
the  concrete  at  the  cracks.  This  working  results  in  spalling  and  faulting  at 
the  cracks  and  joints,  which  produces  a  rough  riding  surface  with  debris 
formed  by  the  spalling  of  the  concrete.  The  formation  of  the  crack  at  the 
bottom  of  the  slab  is  assumed  to  be  a  function  of  the  magnitude  of  the  stress 
and  the  fatigue  effects  of  the  repeated  loading.  The  rapidity  at  which  the 
crack,  once  initiated,  migrates  through  the  slab  is  found  to  be  a  function  of 
the  foundation  strength  and  number  of  load  repetitions.  This  mode  ef  failure 
has  been  observed  in  bo\h  operational  and  experimental  pavements.  It  is  rec¬ 
ognized  that  failures  of  the  materials  supporting  the  concrete  slab,  such  as 
pumping,  frost  action,  swelling,  shrinkage,  etc.,  can  precipitate  cracks  in 
the  pavement  slab;  current  desig:.  procedures  attempt  to  prevent  these  effects 
through  the  proper  design  and  construction  oi  the  supporting  medium.  In  gen¬ 
eral,  a  pavement  slab  will  fail  in  flexure  before  sufficient  distress,  such  as 
shear  failure,  can  be  generated  in  the  supporting  medium.  Only  under  excep¬ 
tional  conditions  are  cracks  initiated  in  the  surface  cf  the  slab  and  migrate 
downward  due  to  traffic. 
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Current  CE  evaluation  and  design  methods  are  based  on  stress  in  the 
concrete  pavement  as  calculated  from  the  Westergaard  analysis  (reference  32). 
Westergaard  assumed  the  slab  to  be  an  elastic,  homogeneous,  isotropic  uniform 
plate  supported  on  a  dense  liquid.  The  p^ *te  was  either  infinite  in  horizon¬ 
tal  extent  or  semi-infinite  (one  free  edge)  in  horizontal  extent.  Calcula¬ 
tions  and  numerous  tests  have  indicated  that  the  stresses  produced  by  loading 
near  a  free  or  jointed  edge  are  equal  to  or  greater  than  stresses  produced  by 
loading  in  the  interior  of  a  slab.  Thus  the  edge-loading  assumption  was 
adopted  as  the  basis  for  the  CE  design  method.  In  practice,  come  means  of 
load  transfer  is  provided  in  the  joints  between  individual  slabs  or  the  edges 
of  the  slab  are  thickened  to  resist  the  higher  eoge  stresses.  Load-transfer 
mechanisms  commonly  used  are  keys,  dowels,  or  aggregate  interlock  th,?t  devel¬ 
ops  from  the  crack  through  the  cui^rete  thickness,  such  as  at  the  weakened- 
plane  (dummy  groove)  joints.  Laboratory  model  tests  and  prototype  aircraft 
loading  tests  in  ine  field  have  been  employed  to  determine  the  amount  of  load 
that  can  be  transferred  from  the  loaded  slab  to  the  adjacent  unloaded  slab 
through  the  load-transfer  mechanism.  I’nder  normal  conditions,  it  has  been 
found  that  the  load  transfer  can  vary  from  a  maximum  of  nearly  50  percent  to  a 
minimum  of  25  percent.  The  minimum  value  of  25  percent  has  been  selected  for 
use  in  the  CE  criteria,  and  when  computing  the  stress  due  to  a  load  tangent  to 
a  joint  having  a  load- transfer  mechanism,  the  applied  load  is  reduced  by  25 
percent.  This  value  varies  depending  upon  such  factors  ps  the  joint  opening, 
subgrade  strength,  warping  and  curling  within  the  slabs,  degree  cf  distress  in 
the  slabs,  and  condition  of  the  load -transfer  device  constructed  in  the  pave¬ 
ment;  however,  longtime  performance  studies  have  indicated  that  25  percent  is 
a  reasonable  value. 

Westergaard' s  equation  for  calculating  free-edge  stress  is  as  follows: 
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where 


a  =  frec-edge  stress  due  to  single  wheel  tangent  to  edge,  psi 
u  =  Poisson's  ratio  of  slab  material 
P  =  load,  lb 


*  rs 

*K 


AFWL-TR-70-1 L3 


h  «  thickness  of  slab,  in. 

E  «  modulus  of  elasticity  of  slab  materials,  psi 
k  =  modJilus  of  subgrade  reaction  lb/in. ‘"/in. 

a,b  a  semimajor  and  semimincr  axes,  respectively,  of  elliptical  contact 
area,  in. 

t  =  radius  of  relative  stiffness  of  slab  and  subgrade,  in,,  where 

i 

For  multiple-wheel  gears,  the  equations  developed  by  Westergaard  be¬ 
come  rather  complex.  A  graphical  solution  of  these  equations  in  the  form  of 
influence  charts  was  developed  by  Pickett  and  Hay  (reference  5-1) «  Current  CE 
practice  utilizes  results  found  from  these  influence  charts  fc'  multiple-wheel 
assemblies.  The  calculated  values  are  again  adjusted  to  account  for  load 
transfer  across  joints. 

Fatigue  or  repeated  load  effects  are  accounted  for  in  the  design 
method  by  establishing  an  allowable  flexural  stress  for  a  specified  number  of 
load  repetitions  based  on  design  factor  DF-versus-traffic  volume  relation- * 
ships.  The  design  factor  is  defined  as  the  ratio  of  concrete  flexural 
strength  to  the  maximum  edge  stress  as  calculated  by  use  of  equation  21,  with 
allowances  given  for  load  transfer  across  joints. 

From  the  results  of  full-scale  accelerated  traffic  tests  utilizing 
single-wheel  gear  loadings  ranging  from  25,000  to  150,000  lb,  twin- wheel 
gear  loadings  ranging  from  60, 00C  to  100,000  lb,  and  twin- tandem-wheel  gear 
loadings  ranging  from  150,000  to  325 >000  lb,  DF-versus-traffic  volume  relation¬ 
ships  have  been  developed  for  three  separate  conditions  of  failure  in  plain 
concrete  pavements.  These  failure  conditions  have  been  designated  as  ’'ini¬ 
tial,"  "shattered  slab,"  and  "complete."  The  DF-versus-traffic  volume  for 
each  of  these  failure  conditions  is  shown  in  figures  12  through  l4.  In  these 
figures,  the  traffic  volume  has  been  expressed  in  terms  of  coverages,  a  cover¬ 
age  being  defined  as  when  each  point  in  the  selected  traffic  area  width  has 
been  subjected  to  a  maximum  stress  repetition  by  the  gear  in  question.  Since 
the  aircraft  gears  are  composed  of  a  variety  of  combinations  of  wheels,  wheel 
sizes,  wheel  spacings,  and  arrangements  of  the  wheels  and  gears,  it  is  obvious 
that  the  traffic  volume  (number  of  aircraft  passes)  represented  by  any  one 
coverage  level  is  variable,  depending  upon  the  aircraft  gear  configuration. 

Thus  to  convert  from  coverages  to  aircraft  passes  or  vice  versa  requires  the 
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COVERAGES 


Figure  12.  Desigxi  Factor  Versus  Coverages  for  Initial  Failure  Condition. 

Pass  per  Coverage  Ratios  for  12-Wheel  and  Twin-Tandem  Assem¬ 
blies  =  1.34  and  3.30,  Respectively 


COVERAGES 


Figure  13.  Design  Factor  Versus  Coverages  for  Shattered- Slab  Failure 
Condition.  Pass  per  Coverage  Ratios  for  12-Wheel  and  Twin- 
Tandem  Assemblies  =1.34  and  3-30,  Respectively 
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determination  of  a  pass  per  coverage  (p/c)  ratio  that  is  constant  for  any  one 
aircraft.  The  development  of  the  p/c  ratio  is  illustrated  in  Section  IV  lc(2) 
of  this  report.  In  this  section,  the  development  of  the  p/c  ratio  is  specific¬ 
ally  for  flexible  pavements  wherein  a  coverage  is  defined  as  sufficient  con¬ 
tiguous  wheel  passes  that  every  point  in  a  traffic  lane  has  been  touched  one 
time  by  a  moving  wheel.  The  p/c  ratio  is  the  same  for  both  types  of  pavements 
for  single  wheels,  for  multiple  wheels  where  the  wheels  are  arranged  abreast, 
and  for  multiple  wheels  where  the  wheels  are  arranged  in  tandem  when  the  tan- 
den  spacing  is  sufficiently  large  to  result  in  separate  maximum  stresses  in 
the  concrete  under  each  tandem  wheel.  However,  when  the  tandem  spacing  is 
small,  only  one  maximum  stress  occurs  in  rigid  pavements  for  each  pass  of  the 
gear,  while  in  flexible  pavement,  the  pavement  is  subjected  to  two  moving 
wheel  passes.  For  these  situations,  the  p/c  ratio  for  rigid  pavements  is  two 
times  the  p/c  raLio  for  flexible  pavements.  Figures  15  and  16  show  strain 
profiles  in  the  concrete  as  the  12-wheel  and  twin-tandem-wheel  gears  pass  over 
an  instrumented  point  in  the  pavement.  For  the  12-wheel  gear,  it  can  be  seen 
that  there  are  four  excursions;  however,  only  two  are  maximum.  Thus  for 
rigid  pavements,  a  pass  of  the  12-wheel  gear  produces  two  maximum  stresses, 
while  for  flexible  pavements,  a  point  on  the  pavement  is  subjected  to  four 
wheel  passes.  For  the  twin-tandem  gear,  it  can  be  seen  that  there  are  two 
strain  excursions  of  almost  equal  magnitude  with  very  little  relaxation  be¬ 
tween  the  wheels.  Thus  for  one  pass  of  this  gear,  there  is  only  one  maximum 
stress  repetition  in  rigid  pavements,  while  in  flexible  pavements,  there  are 
two  wheel  passes.  Therefore,  for  both  the  12-wheel  and  twin-tandem-wheel 
gears  used,  the  p/c  ratio  for  rigid  pavements  is  twice  the  p/c  ratio  for  flex¬ 
ible  pavements. 

Figure  12  represents  the  DF-versus-traffic  volume  relation  used  for  CE 
plain  concrete  pavement  design  criteria.  The  initial  failure  condition  repre¬ 
sented  is  the  point  at  which  the  crack,  which  originates  at  the  bottom  of  the 
slab  and  migrates  upward  to  the  surface,  commences  to  spall  and  ravel,  which 
produces  debris  on  the  surface.  This  failure  is  more  a  functional  failure 
than  a  structural  failure  since  the  pavement  slab  will  continue  to  support  the 
loaded  aircraft  after  the  Initial  failure  condition  has  been  reached.  Test 
data  result  in  a  family  of  curves  depending  upon  the  modulus  of  subgrade  re¬ 
action;  however,  the  range  of  the  DP's  for  subgrade  moduli  from  25  to  200  lb/ 

3 

in.  is  small  and  for  practical  purposes  is  represented  by  a  single  curve. 

Figure  13  represents  the  DF-versus-traffic  volume  relation  for  a 
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failure  condition  designated  as  shattered  slab  (sometimes  referred  to  as  three 
or  four  cracks  per  slab).  This  relation  is  used  primarily  for  rigid  pavement 
evaluation  for  operational  conditions  where  greater  structural  deterioration 
can  be  tolerated  to  accommodate  overloading. 

Figure  lJi  represents  the  DF-versus-traffic  volume  relation  for  a  fail¬ 
ure  condition  designated  as  complete.  The  curves  represent  a  condition 
wherein  the  concrete  lias  been  cracked  into  pieces  of  about  20  to  25  sq  ft  in 
area.  Test  data  show  that  shortly  after  Uns  condition  has  been  achieved, 
shear  failures  and  excessive  deformations  begin  to  occur  and  danger  of  a  punch¬ 
ing  shear-type  failure  is  imminent.  This  relationship  forms  the  basis  for  the 
CE  nor~igid  overlay  design  criteria.  That  is,  when  a  nonrigid  overlay  is  used, 
the  base  pavement  is  allowed  to  crack  during  traffic  to  about  the  condition 
represented  by  this  figure. 

b.  Nonrigid  Overlay 

Current  CE  nonrigid  overlay  thickness  requirements  for  strengthening 
existing  rigid  pavements  are  completely  empirical,  having  been  developed  from 
full-scale  accelerated  traffic  tests.  Several  test  tracks  have  been  con¬ 
structed  and  traffic  tested  to  failure  wherein  the  load,  gear  configuration, 
thickness  of  existing  rigid  pavement,  and  thickness  and  types  of  nonrigid  over¬ 
lay  pavement  have  been  varied.  From  these  tests,  it  has  been  found  that  when 
the  nonrigid  overlay  materials  are  of  high  quality  and  are  well  constructed 
traffic  can  be  safely  applied  until  the  rigid  base  pavement  has  broken  into 
pieces  of  about  20  to  25  sq  ft  in  area.  S^me  reflection  cracking  in  the  non- 
rigid  overlay  induced  by  the  movement  at  joints  and  cracks  in  the  base  pave¬ 
ment  may  be  experienced  which  creates  minor  maintenance  requirements;  however, 
the  system  is  structurally  safe  and  will  provide  a  relatively  smooth  operat¬ 
ing  surface  until  the  base  pavement  has  deteriorated  to  the  above-described 
condition.  Continued  traffic  beyond  this  point  will  result  in  further  crack¬ 
ing  and  deterioration  of  the  base  pavement  and  eventually  exc  -ssive  deflec¬ 
tions,  faulting,  and  even  punching  shear-type  failure.  From  traffic  test  re¬ 
sults,  it  was  determined  that  the  required  thickness  of  nonrigid  overlay  could 
be  related  to  the  deficiency  in  plain  concrete  thickness  that  would  be  re¬ 
quired  to  carry  the  load  and  volume  of  traffic.  The  required  thickness  of  non- 
rigid  overlay  can  be  expressed  as  2.5  times  the  difference  between  the  existing 
rigid  pavement  thickness  and  the  thickness  of  plain  concrete  pavement  that 
would  be  required  for  the  loading  and  traffic  volume  based  upon  the  complete 
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failure  concept.  This  relation  is  expressed  as  foil  *s: 


4- 


2.5 (Fhd 


(22) 


where 


t  -  required  thickness  of  nonrigid  overlay,  in. 

F  -  factor  expressing  the  percentage  reduction  in  required  plain  concrete 
thickness  from  the  initial  to  the  complete  failure  condition.  This 
factor  is  dependent  upon  both  subgrade  modulus  and  volume  of  traffic. 

h^  =  required  thickness  of  plain  concrete  for  initial  failure  condition, 

C  =  factor  dependent  upon  the  structural  condition  of  the  existing  rigid 
pavement  prior  to  overlay.  C  =  1.0  when  the  existing  slabs  have 
only  initial  cracking  and  0.75  when  existing  slabs  have  two  or  more 
cracks . 
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SECTION  IV 

ANALYSIS  OF  MWHGL  DATA 


FLEXIBLE  PAVEMENT 


Comparison  with  Existing  Criteria 


The  analysis  of  the  MWHGL  test  section  data  was  based  on  the  as- 
constructed  values  shown  in  table  1,  the  traffic  test  data  shown  in  table  2, 
and  the  vertical  deflection  measurements  obtained  from  instrumentation  in 
items  3  and  4.  Tables  1  and  2  are  from  Volume  II. 

A  plan  layout  is  shown  in  figure  17 .  The  test  section  was  designed 
based  on  medium- load  pavement  design,  consisting  of  a  design  loading  of  100 
kips,  twin  wheels  spaced  37  in.  center-to-center,  and  a  contact  area  of  267 
sq  in.  Items  1,  2,  3,  4,  and  5  were  designed  for  8,  40,  200,  200,  and  1000 
coverages,  respectively,  based  on  a  subgrade  of  4  CBR.  The  reason  for  this 
design  stemmed  from  the  requirement  to  determine  if  a  medium- load  pavement 
was  adequate  for  operation  of  C-5A  aircraft.  The  existing  design  equation, 


was  used  to  predict  (see  equation  20)  thickness  requirements  for  the  above 
coverage  levels.  The  equivalent  single-wheel  load  (ESWL)  for  the  12-wheel 
gear  was  determined  using  the  single-layer  elastic  theory  with  a  20-radii  cut¬ 
off  modification  incorporated  (see  figure  18).  This  modification  adjusts  the 
deflection  based  on  a  subtraction  of  the  prescribed  constant  deflection  factor 
at  the  desired  offset,  20  radii.  The  ESWLfs  for  other  wheel  gear  configura¬ 
tions  were  computed  using  the  single-layer  elastic  theory  with  no  cutoff  modi¬ 
fication  (i.e.,  20-radii  cutoff  for  C-5A  and  single-layer  no  cutoff  for  Boeing 
747).  These  two  methods  of  determining  ESWL  represented  design  procedure  at 
the  time  that  the  test  section  was  designed. 

The  as-constructed  subgrade  CBR  varied.  The  as-constructed  thick¬ 
nesses  were  the  same  as  design  thicknesses  except  In  item  5.  Therefore,  the 
design  was  adjusted,  based  on  the  as-constructed  values  (table  1).  Using  the 
design  equation  and  as-constructed  values  of  CBR  and  thickness,  the  design 
coverages  became  9>  79*  206,  248,  and  810  for  items  1,  2,  3*  4,  and  5,  respec¬ 
tively.  A  comparison  of  existing  criteria  with  the  performance  of  the  MWHGL 
test  section  is  shown  in  table  3. 
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189.22 

186.13 

187.29 

186.19 

185-23 

4.0 

2.7 

3.1 

2.7 

2.0 

31.8 
’4 . 1 

31.5 

£3.5 

:-6 

87.5 

83.5 

84.1 

98.1 

98.6 

88 

84 

88 

98 

99 

100 

100 

96 

100 

100 

kaat 

naneuver 

area 

Cr  atone 

0*30 

192.20  ' 

0.63 

7.5 

33 

32 

41 

1  ■> 

138.9 

140 

99 

Subbaaa 

191.57 

1.91 

23.0 

16 

8.9 

126.5 

126 

100 

Mirada  (CL) 

185.66 

-  '■ 

8 

20.5 

103.4 

.  103 

100 

a  Tayer  thlcteaim  Matured  at  tart  pit  io^ationa. 

**  total  actual  tMctoeaaea  derived  fmi  eroaa  tact  Iona  taken  during  coortruction  at  the  toy  of  #a:h  pavanant  eleneot. 
t  Laboratory  denaitiaa  ahown  la  tbla  coluwrt  ara  tha  (3  55  naxiMun  danaltla*  at  optima  water  content  for  the  eruahad-etoa*  bate  and 
tha  gravelly  seed  aubbaaa  Material  and  tha  CS  55  dennlty  eorreaponding  to  tha  field  In-place  water  content r  for  the  robgrade 
Material. 
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CSWL,  %  ASSEMBLY  LOAD 

Figure  18.  Comparison  of  Measured  and  Computed  Equivalent  Single- 
Wheel  Loads  Versus  Depth  for  12-Wheel  Assembly  (Data 
from  Volume  III,  Item  3>  Static  Load) 
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Table  3 

Comparison  of  Existing  Criteria  with  Performance 
of  Multiple -Wheel  Heavy  Gear  Load  Test  Section 


ESWL* 

Thickness  t,  in. 

Failure 

Coverages 

Identification 

CBR 

lb 

Actual 

Indicated** 

Actual 

Predicted 

Lane  1,  item  1 

3 .7 

52,200 

15 

14.5 

8 

9 

item  2 

4.4 

63,000 

24 

25.1 

104 

79 

item  3 

3-8 

74,900 

33 

42.6 

1500 

206 

item  4 

4.0 

74,900 

33 

41.5 

1500 

248 

item  5 

4.0 

84,250 

4l 

48.8 

3850+ 

810 

lane  2,  item  1 

3*7 

50,000 

15 

13.1 

6 

10 

item  2 

4.4 

50,000 

24 

24.6 

200 

170 

Lane  3A,  item  1 

3.7 

30,000 

15 

19.0 

120 

33 

item  2 

4.4 

30,000 

24 

20.8 

450t 

1438 

Lane  3B>  item  3 

3.8 

124,800 

33 

32.6 

40 

43 

item  4 

4.0 

124,800 

33 

31.8 

4o 

49 

item  5 

4. 0 

i44,ooo 

4l 

47.0 

280 

U3 

*  Lane  1  ESWL  computations  made  using  elastic  theory  with  20-radii  cutoff 
modification.  All  other  lanes  ESWL  computed  us^ag  elastic  theory  with  no 
modification. 

**  Thickness  indicated  by  number  of  coverages  at  failure, 
t  Nonfailure . 
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After  testing,  the  failure  coverages  were  used  in  the  design  equation 
to  determine  the  thicknesses.  These  values  are  shown  in  column  5  of  table  3 
(indicated  Thickness) ,  A  comparison  of  indicated  and  actual  thicknesses  re¬ 
flects  the  close  agreement  between  predictions  based  on  the  present  criteria 
and  the  actual  behavior  of  the  thinner  items  at  low  coverage  levels.  The 
agreement  is  not  as  good  for  lane  1,  item  5,  for  which  present  criteria  indi- 
cate  that  8l0  coverages  (t  =  4l.O)  should  have  failed  the  item,  but  the  item 
actually  withstood  ?v°/50  coverages  (t  ~  hi  in.)  without,  developing  a  failure. 

Also  in  lane  3B,  item  5*  present  criteria  indicate  a  required  thickness  of 
U7.0  in.  for  280  coverages,  but  actually  hi  in.  sufficed.  A  plot  of  thickness 
versus  failure  coverages  for  the  12-wheel  loading  shown  in  figure  19  shows 
this  divergence  at  higher  coverage  levels.  Curve  A  represents  actual  perform¬ 
ance  data,  and  curve  B  represents  design  thickness  requirements  determined  by 
the  existing  criteria. 

b.  Equivalent  Single-Wheel  Load 

The  computation  of  ESWL  for  use  in  predicting  the  behavior  of  uhe 
test  items  with  traffic  was  investigated  next.  Using  the  relation 

t  -  (0.23  log  e  +  0.15)^^^  -  I”  (23) 

the  ESWL  was  computed  by  rearranging  the  behavior  pattern  to  the  following 
form: 

ESWL  -»  8.1  CBR  (  -  +  - — - 5  )  (2i) ) 

\n  (0.23  log  C  +  0.15)  / 

The  ESWL  was  then  calculated  using  field  data  and  the  present  flexible  pavement 
design  method.  The  results  of  this  computation  for  five  coverage  levels  are 
shown  in  table  .  These  computations  show  a  decrease  in  the  32-wheel  ESWL  be¬ 
ginning  at  shallow  depths  as  shown  in  the  plot  of  t  versus  ESWL  for  5000  cov¬ 
erages  in  figure  20.  This  curve  also  shows  a  discontinuity  at  shallow  depths. 
The  indicated  decrease  in  ESWL  with  depth  is  impossible  by  definition.  The  ra¬ 
tio  of  the  deflection  under  a  multiple-wheel  assembly  to  the  deflection  under 
any  single -wheel  assembly  must  approach  unity  with  depth  as  the  deflections  un¬ 
der  each  assembly  approach  the  same  value;  therefore,  the  present  flexible  pave¬ 
ment  design  equation  does  not  accurately  predict  ESWL  for  use  in  this  analysis. 

Deflections  determined  through  use  of  the  elastic  theory  yield  values 
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Figure  19 •  Comparison  of  Actual  Data  (Curve  a)  with  Existing  Design 
Curve  (Curve  B)  for  Thickness-Failui e  Coverage  Relation¬ 
ship,  3.2-VJheel  Gear 
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Table  4 

Computed  ESWL's  and  Depths 


25,000 

Coverages 

5000 

Coverages 

1000 

Coverages 

200 

Coverages 

- Iio - 

Coverages 

t. 

ESWI, 

t 

ESWL 

t 

ESWL 

t 

ESWL 

t 

ESWL 

in. 

!■ 

in. 

i 

in. 

$ 

in. 

JL 

in. 

$ 

5.8 

35 .8 

5 

47.3 

5.04 

49.5 

4.76 

51.0 

5.2 

>*3.5 

7.0 

OJ 

OO 

OJ 

t 

34 1 

5.88 

35.8 

5.44 

40.0 

6.24 

32.2 

8.1 

24.5 

7 

28.2 

6.72 

29.5 

6.12 

33.0 

7.80 

25.5 

9-3 

21.8 

8 

24.8 

7.56 

26.0 

6.80 

29.2 

10.4 

20.0 

10.5 

19.8 

9 

22.2 

8.4 

23.5 

8.16 

24.5 

13.0 

17.6 

11.6 

18.7 

10 

20.65 

10.08 

20.5 

10.2 

20.2 

15.6 

16.5 

13.9 

17.2 

12 

18.45 

12.6 

17.8 

13.6 

17.3 

18.2 

15.5 

17.4 

15.8 

15 

16.75 

16.8 

16.0 

17.0 

15.9 

20.8 

15.0 

23.2 

14.6 

20 

15.18 

21.0 

15.0 

20.4 

15.1 

23.4 

l4.6 

29.O 

14.0 

25 

14.38 

25.2 

14.4 

23.8 

14.5 

26.0 

14.3 

34.8 

X3-^ 

30 

13.88 

29.4 

14.0 

27.2 

i4.o 

28.6 

14.0 

40.6 

13.1 

35 

13.5 

33.6 

13.55 

30.6 

13.8 

31.2 

13.7 

46.4 

12.7 

40 

13.25 

37.8 

13<3 

34.0 

13.5 

33.8 

13-5 

45 

12.9 

42.0 

13-0 

37.4 

13.3 

36.4 

13.4 

50 

12.5 

46.2 

12.8 

40.8 

13.2 

55 

12.25 

44.2 

12.9 

60 

11.80 

47.6 

12.7 

65 

11.50 

70 

11.20 

Note:  The  present  flexible  pavement  equation  was  used  in  confutations  of  t 
and  ESWL. 
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in  an  elastic,  isotropic,  homogeneous  half-space  under  a  uniformly  distributed 
circular  load,  and  the  computations  will  always  indicate  some  deflection  at 
any  distance  from  the  loaded  area.  In  reality,  the  deflection  basin  atten¬ 
uates  to  insignificantly  small  values.  In  an  effort  to  duplicate  the  theo¬ 
retical  ESWL  that  most  closely  represents  the  behavioral  pattern  based  on  the 
results  of  actual  performance,  different  offset  distances  in  radii  of  negli¬ 
gible  effect  were  assumed,  and  ESWL’s  were  computed.  These  ESWL's  were  com¬ 
puted  using  modified  deflection  curves  in  which  the  theoretical  deflections 
for  the  12 -wheel  assembly  attenuated  to  zero  at  the  respective  radii  distances. 
Figure  1 8  shows  the  result  of  these  computations.  Past  experience  has  shown 
the  need  for  this  nonaceumulation  of  small  deflections  at  given  offsets  from 
the  loaded  area.  This  is  especially  important  when  a  large  number  of  wheels 
are  involved  and  is  one  method  of  representing  field  performance  with  a  modi¬ 
fied  theory.  This  modification  is  accomplished  by  a  subtraction  of  a  constant 
deflection  factor  with  depth  at  the  selected  offset  distance.  The  following 
deflection  factor  modification  constants  were  used  for  10-,  12-,  15- ,  and  20- 
radii  cutoffs,  respectively:  0.0077,  0.0053,  C.003^,  and  0.0019.  The  de¬ 
flection  modification  factors  were  computed  as  follows: 


wheie  ~  ; 

F  =  deflection  modification  factor 

W  =  vertical  deflection,  in. 

E  =  elastic  modulus,  psi 

p  =  average  contact  pressure,  psi 
0 

r  ~  radius  of  loaded  circular  area,  in. 


A  computer  grid  and  coordinate  syslem  was  selected  by  an  iterative 

a- 


process  to  select  the  location  of  maximum  deflections  with  depth.  Migration 
of  deflection  with  depth  is  shown  in  figure  21.  The  theoretical  maximum 
deflection  occurred  under  the  tire  at  the  surface  (point  A).  The  maximum 
deflection  actually  never  reached  the  centroid  of  the  six.  wheels  (point  B) 
within  the  dep+h  of  12  ft ,  but  would  "with  increased  depth  and  would,  at 
greater  depths,  approach  the  centroid  of  the  12  Wheels  (point  C) .  The  instru¬ 
mentation  results  for  vertica3.  deflections  with  depth  for  item  3  are  shown  in 
figure  22,  and  the  ESWL’s  were  computed  from  this  plot.  The  ESWL's  were 
determined  from  figure  22  in  the  folic  wing  manner. 
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LEGEND 


O  PTA,  CENTROID  OF  TIRE 

A  PT  B,  CENTROID  OF  THE  FRONT 
SIX  WHEELS 

□  PT  C,  CENTROID  OF  TWELVE  WHEELS 

Figure  21.  General  Flow  of  Maximum  Vertical  Deflection  with  Depth 
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For  example,  at  9- in.  depth 

ESWL,  $  assembly  =  x  ^  x  100  =  18.28^ 

The  ESWL’s  are  shown  in  figure  18  for  comparison  of  the  measured  and  computed 
values. 

Of  the  relationships  of  ESWL  with  depth  shown  in  figure  18,  the  10- 
radii  cutoff  values  agree  best  in  the  thickness  design  relationship, 

t  -  (0.23  log  C  +  0.15)  £  (23> 


as  shown  in  the  following  tabulation. 


Item 

CBR 

Coverage 

ESWL 

lb 

A 

sq  in 

Thickness 

Calculated 

(Approx) 

,  in. 

Actual 

1 

3.7 

8 

44,000 

285 

13.3 

15 

2 

4.4 

io4 

51,500 

285 

22.6 

24 

3 

3.8 

1500 

59,400 

285 

37.7 

33 

4 

4.0 

1500 

59,400 

285 

36.7 

33 

5* 

4.0 

>3850 

67,250 

285 

43.4 

4i 

*  Nbnfailure. 


The  calculated  t  is  approximate,  since  the  ESWL  is  dependent  on  the  calcu¬ 
lated  t  ♦  and  the  value  is  shown  only  for  first-order  aoourR.cyr 

A  similar  type  of  approach  was  then  directed  toward  the  data  from 
traffic  by  the  twin-tandem  gear  spaced  44  by  58  in.  on  test  items  of  lane  3B. 

A  plot  of  ESWL  (percent  of  assembly  load)  versus  depth  from  computations  using 
elastic  theory,  different  modifications  of  elastic  theory,  and  instrumentation 
deflection  values  is  shown  in  figure  23.  it  is  noteworthy  that  the  instrumen¬ 
tation  ESWL  curve  for  the  twin-tandem  gear  loading  more  closely  approximates 
the  values  derived  from  the  e1 astic  theory  than  those  derived  from  the 
modified  elastic  theory.  However,  the  12 -wheel -gear-loading  ESWL  curve  is 
best  approximated  by  the  modified  elastic  theory  20-radii  cutoff. 

Analysis  of  traffic  test  data  with  the  two  gear  configurations  did  not 
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yield  values  of  F-SWL  consistent  wi'jh  the  elastic  theory  or  any  one  modifica¬ 
tion  thereof.  To  prevent  the  necessity  of  assigning  different  methods  of 
computing  the  ESWL  values  for  different  gear  conf i jurations ,  the  elastic 
theory  without  modification  was  selected  as  best  representing  the  relation¬ 
ship  of  ESWL  for  all  multiple-wheel  assemblies.  This  decision  was  also  sup¬ 
ported  by  related  work  (reference  2*9).  The  difference  in  the  behavior  of  the 
nonmoving  single-wheel,  twin-tandem,  and  12 -wheel  assemblies  can  be  related 
to  their  different  configurations,  spacing,  and  consequently,  different  load 
distributing  characteristics.  However,  the  effects  of  total  number  of  wheels, 
method  of  application  of  traffic,  and  load  repetitions  must  also  be  considered. 

c.  Traffic  Study 

(l)  Repetitions  effect.  The  factor  f  ,  which  equals  (0.23  log  C 
+  0.15),  employed  in  the  current  design  of  flexible  pavements  accounts  for 
the  known  effect  of  traffic  repetitions.  This  factor  was  developed  from 
previous  studies  (reference  38).  Using  the  MWHGL  test  section  data  and  the 
elastic  theory  for  ESWL  determination,  the  following  analysis  was  made.  The 
flexible  pavement  equation  was  rewritten  as 


t 


a 


4 


ESWL  A 
8.1(CBR)  ’  rr 


(27) 


then 


t 


a 


or 


4 


ESWL  A 
8.l(CBR)  '  it 


where 

a  -  load  repetitions  factor 
t  -  actual  test  section  thickness,  in. 
ESWL  =  ESWL,  lb,  at  the  depth  t 
CBR  =  strength  value  of  supporting  layer 
A  =  measured  tire  contact  area,  sq  in. 


(28) 


The  a  factors  or  the  effect  of  repetitions  on  the  total  thickness  require¬ 
ments  were  computed  and  are  shown  in  table 

Figure  2^1  shows  a  semilog  plot  of  the  a  values  versus  coverages. 
For  comparison,  the  existing  relationship  of  f  versus  coverages  is  also 
shown  in  figure  2k,  The  existing  criteria  for  low- intensity  traffic  seems 
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Table  5 

Comparison  of  a  and  f  with  Coverages 


Identification 

Failure 

Coverages 

Actual 
Thickness 
t  in. 

ESWL,* 

lb 

CBR 

A 

sq  in. 

at 

f 

Lane  1,  item  1 

8 

15 

62,280 

3.7 

285 

0.336 

0.357 

item  2 

104 

24 

77,760 

4.4 

285 

0.525 

0.614 

item  3 

1500 

33 

94,320 

3.8 

285 

0.605 

0.880 

item  4 

1500 

33 

9^,320 

4.0 

285 

0.621 

0.880 

item  5 

3850** 

4l 

108,360 

4.0 

285 

0.719 

0.975 

Lane  2,  item  1 

6 

15 

50,000 

3.7 

285 

0.378 

O.329 

item  2 

200 

24 

50,000 

4.4 

285 

0.662 

0.679 

Lane  3A,  item  1 

120 

15 

30,000 

3.7 

285 

0.497 

0.628 

item  2 

450** 

24 

30,000 

4.4 

285 

0.875 

0.760 

Lane  3B,  item  3 

4o 

33 

124,800 

3.8 

290 

0.524 

0.518 

item  4 

40 

33 

124,800 

4.0 

29c 

0.536 

0.518 

item  5 

280 

4l 

144,000 

4.0 

290 

0.621 

0.713 

5 


*  Obtained  from  figure  18  or  23,  as  appropriate.  All  computed  using  elastic 
theory  without  modification. 

**  Nonfailure . 
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adequate  for  all  wheel  gears  shown  (reference  52).  However,  as  noted,  with  an 
increase  in  coverages,  the  present  load  repetition  factor  is  too  large  and 
hence  yields  too  great  a  thickness.  Curves  1,  2,  and  3  in  figure  24  represent 
failure  envelopes  for  the  single-wheel,  twin-tandem,  and  12 -wheel  gears,  respec¬ 
tively.  An  additional  point  was  obtained  fbr  curve  3  for  the  12 -wheel  as¬ 
sembly.  This  was  obtained  by  extrapolating  the  nonfailure  of  item  5  to  fail¬ 
ure  at  15,000  coverages  as  shown  in  curve  A,  figure  19. 

A  statistically  weighted  curvilinear  regression  (references  53 
and  54)  was  performed  on  the  MWHGL  failures,  given  in  table  5>  to  determine 
the  best-fit  curve  for  coverages  and  y  .  The  results  are  shown  in  figure  25 
and  table  6.  Th  ’  results  showing  standard  deviations  and  comparison  of  esti¬ 
mates  with  prediction  equations  and  actv » 1  data  are  shown  in  table  6.  The 
quadratic  equation,  curve  3,  represents  the  best-fit  curve.  Curve  A  is  the 
existing  design  curve.  Curve  B  is  a  good  estimate  for  the  twin-tandem  assembly 
but  requires  toe  much  thickness  for  the  12 -wheel  gear.  If  the  curve  B  were 
shifted  downward  a  constant  amount  to  include  12-wheel  daua,  it  would  not  re¬ 
quire  enough  thickness  for  the  twin-tandem  assembly.  This  would  indicate  a 
need  for  a  separave  y  versus  traffic  volume  factor  plots  for  each  of  the 
assemblies  used. 

Using  the  relationships  for  saown  in  figure  24  and  earlier 
conclusions,  the  following  design  pattern  was  developed  for  thickness  require¬ 
ments  to  protect  supporting  layers. 

orqiJ  =  equation  of  curve  1,  (23a) 

bW  figure  24 

ttmm  =  equation  of  curve  2,  (23b) 

u  figure  24 

Q'touu  =  equation  01  curve  3,  (23c) 

12WH  figure  24 

Using  the  above  relationships,  computations  were  made  to  check  predictions  of 
the  pattern  with  performance  of  the  test  section  (table  7).  As  shown  in 
table  7,  this  method  of  predicting  required  thickness  is  in  agreement  with 
the  performance  of  the  test  pavements. 

(2)  Traffic  distribution  studies.  During  the  early  stages  of  develop¬ 
ment  of  criteria  for  the  design  and  evaluation  of  pavements,  it  became  apparent 
that  a  method  was  needed  to  account  for  repetitions  of  traffic.  To  simply 
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Figure  25.  Comparison  of  Existing  Design  Curve  (A)  and  Best-Fit  Curve  (b) 
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Tubir  t 

Kj-a  .  !  Curvi . ..  e*«?  It  t-Redi 


1 .  First -Order  Equation 


Tern 

i> 

1 

Coefficient  tandard  Devlat 

.*09017  S .  li>7  3  tf — ? 

.110179  .02380;’ 

Standard  Deviation  ■  5 . 8bto92R-2 

,  t 

X  Value 

a  Value 

Y  Estimate 

Real dual 

.90  309 

.336 

.1*08719 

-7.2T19>*E-2 

2.01703 

.525 

.531676 

-6.67561E-3 

3.17609 

.605 

.659612 

-5.1*6ll5E-2 

3.17609 

.621 

.659612 

-3.8611c--2 

.776151 

.378 

.  o  1*929 

-1.6928TB-: 

2. 30103 

.662 

.5'-' 302  3 

9.89771E-2 

2.07918 

.1*97 

.538535 

-V.1535kE-2 

1.60206 

.521* 

.1*85871 

3.01288E-2 

1.60206 

.538 

.1*95871 

5.21268E-2 

2.1*1*716 

.621 

579152 

U.18I.76E-2 

II.  Quadratic 

Equation 

Ter* 

0 

1 

2 

Coefficient 

.125569 

.326C92 

-5.38l5**E-2 

Landard  Deviation 

Standard  Deviation 

9.30189E-? 

9.9519V>-2 

2.l»3l*53E-2 

•  6.81 3lftE-2 

X  Value 

Y  Value 

Y  Estiaate 

Real dual 

.90309 

.336 

.376169 

-V.01695E-? 

2.01703 

.525 

.56V 36k 

-3.936V1 t-2 

3.17609 

.605 

.6l8«*02 

-1.3V018E-2 

3.17609 

.621 

.6181*02 

2.598211-3 

•TT8151 

.378 

.3V6732 

.031268 

2.30103 

.662 

.590978 

T.10215E-2 

2.07918 

Ml 

<57093 

-7.39302E-2 

1.60266 

•  52k 

.509866 

.0lkl3k 

1.60206 

.538 

.509866 

.02813k 

2.*i,i716 

.621 

.CCI29 

1.97O90E-2 

III.  Cubic  Equation 

Term 

Coefficient 

Standard  Deviation 

0 

.257267 

.320773 

1 

6.96VV5E-2 

.603V99 

2 

8.88VT5E-2 

.  33.1537 

3 

-.023712 

5.l*936kK-2 

Standard  Deviation  ■  5.1199l**-^ 


X  Value 

Y  Value 

Y  Fat  1 mate 

Residual 

.90309 

.336 

.375159 

-3.9159VE-? 

2.01703 

.525 

.56V628 

-3.96277E-2 

3.17609 

.60* 

.61501 

-1 . 000963-2 

3.17609 

.621 

.61501 

5.990V3E-3 

.778151 

.378 

.35V087 

2.39126E-2 

2.30103 

.662 

.59905k 

.06291*6 

2.07918 

.V97 

.573028 

-7.602T6E-2 

1.60206 

.32k 

.1*99378 

2.16223E-2 

1.60206 

.530* 

.*•99378 

3.86223E-2 

2.kk?l6 

.621 

.612269 

8.73069E-3 

Motet  X  » 

Y  - 

aircraft  volume  traffic 
load  repetition  factor, 

factor,  coverages. 

or  . 

i 

■5 


!  1 
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Table  7 

Comparison  of  Predicted  Thicknesses  Versus  Test 
Section  Thicknesses 


«< 

f  ESWL  A* 

Thickness 

t  ,  in. 

Identification 

Factor 

\  8.1  GBR  *  n 

Predicted 

Actual 

Lane  1,  item  1 

0.340 

44.5 

15.1 

15 

item  2 

0.550 

45.7 

25.1 

24 

item  3 

0.661 

54.5 

36.0 

33 

item  4 

0.661 

53.1 

35.1 

33 

item  5** 

O.692 

57.0 

39.4 

4l 

Lane  2,  item  1 

0.360 

39.7 

14.3 

15 

item  2 

0,670 

36.2 

24.2 

24 

Lane  3A,  item  1 

0.625 

30.2 

18.9 

15 

item  2-** 

0.740 

27.4 

20.3 

24 

Lane  3B,  item  3 

0.510 

62.9 

32.1 

33 

item  4 

0.510 

61.3 

31.3 

33 

item  5 

0.640  - 

66.0 

42.2 

41 

*  BSWL  determined  using  elastic  theory  without  modifi cations . 

**  Nonfailure. 


-  8MPMt*^:sc»8rt«JW$Qags^^  , 
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count  the  number  of  aircraft  using  an  airfield  is  not  adequate.  The  incre¬ 
mental  detriment  to  a  pavement  reeulting  from  a  particular  aircraft  wheel  at 
a  specified  location  on  pavements  is  influenced  by  many  factors.  Some  of  the 
more  important  factors  are:  (a)  number  of  wheels,  (b)  wheel  configuration, 

(c)  tire  contact  area,  (d)  tire  inflation  pressure,  and  (e)  location  of  wheel 
on  pavement. 

In  an  attempt  to  normalize  these  factors  so  that  one  number  could 
be  obtained  to  reflect  their  collective  influence  on  the  total  system  of  de¬ 
sign  and  evaluation,  the  concept  of  coverage  was  introduced  where  coverage  is 
defined  as  sufficient  contiguous  wheel  passes  so  that  every  point  in  a  traffic 
lane  has  been  touched  one  time  by  a  moving  wheel. 

In  the  past,  traffic  on  test  sections  was  programmed  so  that 
successive  wheel  paths  were  not  overlapping,  and  an  accurate  value  of  cover¬ 
age;?  resulted.  It  was  assumed  that  coverages  on  a  test  section  could  be 
related  directly  to  coverages  on  airfield  facilities  (runways,  taxiways). 
However,  the  random  traffic  on  runways  and/or  taxiways  can  only  be  conven¬ 
iently  counted  as  aircraft  passes.  The  general  equation  currently  used  for 
determining  pass  per  coverage  ratios  (sometimes  called  operation  per  coverage 
ratios),  given  as  equation  7  in  reference  55  >  is: 

E  -  _£ - -  (29) 

c  0.75  N  Wt  K  y) 


where 

p/c  =  pass  per  coverage  ratio 
W^.  =  tire  width,  in. 

S  =  wheel  spacing,  in. 

W  =  wander  width  (40  in.  assumed  for  channelized  traffic  of  military 
w  aircraft  on  a  taxiway) 

N  =  number  of  wheels  in  traffic  lane 


In  accordance  with  the  traffic  pattern  used  on  the  C-5A  test 
sections  (Volume  II )  where  the  traffic  was  programmed  in  such  a  manner  that 
22  passes  of  the  load  cart  generated  32  coverages,  p/c  equaled  22/32  or 
0.68.  In  comparison,  when  the  dimensions  for  the  C-5A  gear  are  substituted 
in  equation  29  (Wt  -  14.3  in.,  S  =  121  in.,  and  N  =  12),  then  p/c  equals 
1.36.  Thus  the  pass  per  coverage  ratio  determined  from  the  current  method 
for  aircraft  operating  on  runw&ya  and/or  taxiways  does  not  agree  with  the 
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actual  pass  per  coverage  ratio  for  the  test  section.  In  addition,  there  has 
been  considerable  confusion  and  some  misunderstanding  relative  to  the  use  of 
coverage  as  a  means  of  accounting  for  traffic  repetitions  and  some  indication 
that  the  current  method  of  converting  coverages  into  passes  (or  vice  versa) 
may  contribute  to  some  of  the  difficulty  encountered  in  analysis  of  test  sec¬ 
tion  data.  Tor  these  reasons,  a  revised  method,  based  on  statistical  analysis, 
for  determining  pass  per  coverage  ratios  has  been  developed.  The  revised 
method  and  its  application  are  described  in  the  following  paragraphs.  The 
backipround  and  details  of  the  development  and  application  of  the  method  are 
given  in  reference  55* 

(3)  Development  of  revised  traffic  distribution  concepts.  Statis¬ 
tical  analysis  combined  with  actual  traffic  distribution  data  has  been  used  to 
develop  a  new  method  for  determining  pass  per  coverage  ratios.  A  fundamental 
assumption  is  made  that  airfield  traffic  is  normally  distributed  instead  of 
uniformly  distributed  as  was  assumed  in  development  of  equation  29.  Tor  a 
large  number  of  aircraft  passes,  the  lateral  placement  of  the  center  line  (of 
the  aircraft)  will  be,  statistically,  as  shown  in  figure  26.  The  general 
shape  is  assumed  to  be  normal,  while  the  specific  shape  depends  on  the  stand¬ 
ard  deviation  cr  .  The  specific  shape  may  also  be  described  by  prescribing  ^ 
the  wander  width  Ww  ,  which  is  defined  as  that  width  over  which  the  center 
line  of  the  aircraft  traffic  is  distributed  73  percent  of  the  time.  It  has 
been  determined,  bn  the  basis  of  an  analysis  of  a  small  amount  of  actual  mili¬ 
tary  aircraft  traffic  distribution,  that  wander  widths  of  kC  and  80  in.  should 


Figure  26.  Theoretical  Normal  Distribution  of  Aircraft  Traffic 
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be  used  in  determining  pass  per  coverage  ratios  for  taxiways  and  runways, 
respectively.  These  \ulues  represent  the  best  values  obtainable  from  exist¬ 
ing  data  and  are  subject  to  change  if  and  when  additional  actual  traffic 
distribution  data  are  obtained. 

The  definition  for  coverage  was  revised  for  the  traffic  distribu¬ 
tion  study.  It  was  considered  that  coverage  represents  the  maximum  number  of 
tire  prints  or  partial  tire  prints  applied  to  the  pavement  surface  at  that 
point  where  maximum  accvmilation  occurs.  In  figure  26,  the  distribution  curve 
also  represents  the  distribution  of  the  center  line  of  one  wheel  on  the  air¬ 
craft.  When  the  wheel  center  line  is  at  a  =  -(W^2)  ,  wheel  passes  accumu¬ 
late  at  0.  Similarly  when  the  wheel  center  line  is  between  a  =  -(W^/2 )  and 
b  =  W.^/2  ,  there  will  be  accumulations  at  0  equal  to: 


As  can  be  seen,  this  is  approximately  equal  to  (C^)(W^.)  .  It  is  assvmed  in 
equation  30  that  the  effect  of  the  edge  of  a  tire  at  0  is  as  detrimental  as 
the  effect  of  the  tire  center  line  at  0.  This  is  not  necessarily  true  and  a 
further  refinement  could  be  made  through  use  of  a  location  weighting  factor 
as  discussed  in  detail  in  reference  55*  However,  for  the  present  this  refine¬ 
ment  will  not  be  used  and  the  simplified  definition  for  coverage,  C 
»  (C^HW^.)  ,  will  be  used.  This  method  can  be  extended  to  aircraft  having 
many  wheels  by  graphical  addition  of  any  number  of  single -wheel  traffic  distri¬ 
bution  curves,  such  as  that  shown  in  figure  26.  The  area  under  the  cumulative 
traffic  distribution  curve  will  then  represent  the  total,  number  of  wheel 
passes.  Pass  per  coverage  ratios  based  on  an  assumed  normal  traffic  distri¬ 
bution  can  be  developed  as  described  in  the  following  paragraphs. 

The  standard  normal  distribution  (SND)  curve  is  shown  in  fig¬ 
ure  27*  Properties  of  the  SND  curve  are  tabulated  in  numerous  standard  refer¬ 
ences  for  statistics  and  tables  of  probability  functions,  and  these  tabulated 
values  may  be  used  to  determine  the  properties  of  a  general  normal  distribu¬ 
tion  (<2fD)  curve,  figure  28.  In  the  SND,  standard  deviation  a  *  1  ,  area 
under  the  curve  A  =  1  ,  maximum  ordinate  C  =  0.399  >  and  75  percent  A  lies 
between  -1.15  <  z  < +1.15  •  The  GND  curve  may  be  related  to  the  SND  curve 
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Figure  2?.  Standard  Normal  Distribution  (SND)  Curve 


Figure  28.  General  Normal  Distribution  (GND)  Curve  for  Aircraft  Traffic 


using  the  following  substitution. 


where 


(31 


z  *  a  variable  in  SND 
x  =  a  variable  in  GND 
u  »  mean  value  in  GND 
<x  «  standard  deviation  in  (20 

Then  where  the  wander  width  Ww  equals  to  in.  or  2x  ,  X  =  20  in.,  and 
u  ~  0  ,  by  substitution,  1.15  -  ft  ,  or  a*  for  the  GND  =  17.4  .  Now 
the  SND  curve  has  the  equation:  x 
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J^z  =  -i-  e  2 

V2S 


(32) 


and  the  G1JP  equation  has  the  equations 

i  -ifH 


Pt(x)  =  - —  e 

t  -  ^2tt 


X 


(33) 


where  in  substitution  of  z  =  -  becomes: 


and 


Pt(x)  =  ±-{<h) 


x 


(34) 


c  =  —  c 

x  a  z 

X 

where  Cx  is  the  maximum  ordinate  on  the  GND  curve.  Also,  the  area  under 
the  GND  curve  A  is  equal  to: 


Pt(x)dx 


(35) 


/' 

>00 

\  f(.0U))axdz 

J  v 


(36) 


■  /%( 


z)dz  =  1.00 


Note,  however,  that  as  wander  width  is  changed,  must  be  recalculated 
as  indicated  above. 

The  GND  curve  can  now  be  applied  to  the  study  of  real  aircraft 
traffic  distribution.  For  example,  using  an  aircraft  with  a  single-wheel 
tricycle  landing  gear  (figure  29)  and  neglecting  the  nose  gear,  plot  the  gen¬ 
eral  normal  distribution  of  each  wheel,  both  on  the  same  plot,  as  shown  in 
figure  29*  Graphically  adding  these  curves  results  in  a  cumulative  curve 
identical  (for  all  practical  considerations)  with  the  GND  curve.  The  total 
area  under  the  cumulative  curve  is  equal  to  two  wheel  passes  (wp)  or  one 
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MAIN  GEAR  WHEEL 


-SINGLE-WHEEL  SPACING 
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Figure  29*  General  Normal  Distribution  for  Nonoverlapping  Wheels  for  the  Typical 
Single -Wheel  Tricycle  Landing  Gear  Shown  Above 
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aircraft  pass  (ap),  and  the  maximum  number  of  coverages  C  is:  C  =  (C  )(W. ) 

X  t 

and 


E 

c 


(37) 


As  the  wheel  spacing  S  becomes  smaller,  the  GND  curve  for  each 

single -wheel  overlap  and  the  distribution  pattern  shown  in  figure  30  are 

obtained.  In  this  case,  graphical  addition  of  the  individual  single-wheel 

curves  results  in  a  cumulative  distribution  curve  with  a  maximum  ordinate 

greater  than  that  of  either  single -wheel  curve.  The  area  under  the  cumulative 

curve  is  still  2.00.  Now,  however,  the  maximum  number  of  coverages  is  the 

maximum  area  under  the  cumulative  curve  with  a  width  of  W^  .  This  value  may 

not  necessarily  occur  vnder  one  wheel  and  its  location  must  be  determined  by 

trial  and  error.  For  example,  the  maximum  area  may  be  as  shown  in  figure  30 

(cross-hatched  area).  Such  a  refinement  is  not  considered  warranted,  and  for 

simplicity,  the  value  (C^c)(W^)  will  be  used  as  the  maximum  value  for 

coverages*  Values  for  C  developed  in  this  study  have  been  obtained 

graphically.  The  determination  of  C  could  easily  be  a  programmed  deter- 

xc 

mination  using  a  computer,  but  this  has  not  been  done  because  the  determina¬ 
tion  need  be  made  only  once  for  each  aircraft.  Thus 


As  an  aid  in  determining  the  magnitude  of  C  ,  figure  31  has 

been  prepared.  This  figure  shows  C  versus  wheel  spacing  for  40-  and  8o-in. 

xc 

wander  Tridths.  For  a  wander  width  of  40  in.,  C  =C  =  0.0229  when  the 

XC  X 

wheel  spacing  is  greater  than  60  in.  (nonoverlapping  wheel  paths).  When  the 
wheel  spacing  is  zero,  C  =  2C  =  0.0458  (tandem  assembly).  When  the  wheel 

XC  X 

spacing  is  greater  than  zero  and  less  than  60  in. C,_rt  can  be  read  from  fig- 

Xv. 

ure  31  •  For  a  wander  of  30  in.,  C  =  C  =  0.01145  when  the  wheel  spacing 

XC  X 

is  greater  than  156  in*  (nonoverlapping  wheel  paths)*  When  the  wheel  spacing 
is  zero,  C  =  2C  =  0.0229  (tandem  assembly).  When  the  wheel  spacing  is 

XC  X 

greater  than  0  and  less  than  156  in.,  C  can  be  read  from  figure  31*  The 

xc 

values  read  from  figure  31  are  applicable  to  the  following  landing  gear 


configurations  without  exception:  single  conventional,  single  tricycle. 
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-  TREQUENCT  WHEEL  C£*T£ *-(.*€  <MSS«  Pf*  *CN  iWW,  P^(X) 


.  cum/t  a  T've  curve 

C.r/*  ^ 


AREA  UNDER  EACH  CURVE  *  1.00 
AREA  UNDER  CUN  JuATtVE  CURVE  =  1.00 

NOTE:  NOSE  GEAR  NEGLECTEO. 


LATERAL  PLACEMENT  OT  WHKl,  CENTER  LINE,  IN. 


Figure  30.  General  Normal  Distribution  Curve  for  Overlapping  Single  Wheels 
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Figure  31.  Maximum  Ordinate  on  Cumulative  Traffic  Distribution  Curve  for  Two 
Wheels  Versus  Wheel  Spacing 
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2  2 

single-tandem  tricycle,  and  twin  bicycle  and  to  the  following  landing  gear 
configurations  when  the  tread  is  greater  than  60  in.  (for  Ww  =  60  in.)  or 
156  in.  (for  W  -  80  in.):  twin  conventional,  twin  tricycle,  twin-tandem 

p  *  2 

tricycle,  and  twin-twin  bicycle.  For  other  landing  gear  configurations, 

C  can  be  determined  by  plotting  the  individual  wheel  distribution  and 
determining  the  cumulative  distribution  curve  by  graphical  addition  as  demon¬ 
strated  in  figure  32. 

The  versatility  of  this  method  is  demonstrated  with  the  deter¬ 
mination  of  p/c  ratios  for  the  C-5A  and  Boeing  747  landing  gear  wheel  con¬ 
figurations  and  a  wander  of  40  in.  as  shown  in  the  following  paragraphs. 

(a)  p/c  ratio  for  the  C-5A.  The  C-5A  wheel  configuration  is 
shown  in  figure  33*  The  normal  distribution  curves  for  wheel  groups  A,  B, 
and  E  are  shown.  In  figure  32  where  the  cumulative  distribution  curve  was 
obtained  by  graphical  addition  of  individual  wheel  distribution  curves.  By 
inspection,  it  can  be  seen  that  the  maximum  ordinate  for  Group  A  plus  Group  C 
(or  Group  B  plus  Group  D),  on  the  cumulative  curve,  is  the  maximum  ordinate 

fbr  the  assembly.  Thus  C  =  2(0.0500)  =  0.1000  •  For  the  entire  assembly, 

xc 

total  number  of  wheel  passes  »  24.  However,  24  wheel  passes  *  1  aircraft  pass 
Thus, 


1 

(o.ioooHwJ 


Or  as  before  (equation  19) >  p/c  *  .  Therefore,  for  the  C-5A 

where  *  15.1  in.  ,  p/c  =  =  0.662  ,  which  is  in  general 

agreement  with  the  operations  per  coverage  ratio  determined  from,  actual  test 


traffic . 


(b)  p/c  ratio  for  the  747  >  The  Boeing  747  wheel  configuration 
is  shown  in  figure  34.  The  pass  per  coverage  ratio  for  the  Boeing  ,?47  can 
be  determined  by  the  graphical  procedure  used  for  the  C-5A  aircraft.  How¬ 
ever,  since  the  individual  twin-tandem  gear  assemblies  are  spaced  far  enough 
apart  to  prevent  any  overlap  between  twin-tandem  assemblies  for  a  1*0- in. 
wander,  the  value  of  C  can  be  determined  through  use  of  figure  31*  For  & 

Av 

44- in.  twin-wheel  spacing,  the  maximuin  ordinate  for  a  44-in.  spacing  is  0.024 

and  C  =  (2) (0.024)  =  0.048  .  Therefore,  for  the  Boeing  747,  where 
xc 


2Note  that  for  these  configurations,  Cxc  will  be  twice  the  value  read  from 
figure  31  because  there  are  two  wheels  in  line  in  the  direction  of  travel. 
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Figure  33*  C-5A  Landing  Gear  Configuration 
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Wt  “  **•  >  p/C  -  (0.0^)02.91  =  1>62  ‘ 

Pass  per  coverage  ratios  for  the  various  landing  gear  wheel 
configurations  used  in  the  MWHGL  investigations  reported  herein  and  for  other 
similar  data  obtained  from  related  studies  have  been  determined  through  use  of 

equation  36  or  37  a$'d  are  shown  in  tables  8  and  9*  Values  for  C  and  C 

X  xc 

were  read  from  figure  31  or  were  determined  by  the  graphical  procedure  just 
demonstrated* 

Using  these  p/c  factors,  the  test  section  failure  coverages 
were  converted  to  actual  field  facility  failure  passes*  This  conversion  is 
shown  in  table  8.  The  Oj,  factors  versus  passes  were  plotted,  as  shown  in 
figure  35,  and  limiting  failure  envelopes  were  drawn  for  the  12 -wheel,  twin- 
tandem,  and  single-wheel  gear  configurations. 

(4)  Related  studies.  Previous  flexible  pavement  failures  (references 
39-46)  in  which  insufficient  thickness  could  be  determined  as  having  an  effect 
on  pavement  behavior  are  listed  in  table  9*  These  previous  failures  were 
plotted  as  cr^  versus  passes  to  form  a  plot  as  shown  in  figure  36.  Differ¬ 
ences  in  test  objectives,  failure  criteria,  methods  of  rating  strength,  rate  of 
field  measurements  and  observations,  construction  techniques  and  materials, 
and  methods  of  applying  traffic  account  for  some  of  the  variations  shown  in 
figure  36.  These  variables  were  minimised  in  testing  the  MWHGL  test  section; 
therefore,  the  data  shown  in  figure  35  are  more  consistent  and  will  be  used 
with  weighted  emphasis  in  this  analysis. 

(5)  Summary.  The  results  of  the  investigations  discussed  previously 
indicate  that  load  repetition  factors  decrease  at  a  reduced  rate  with 
repetitions  for  the  single-wheel,  twin-tandem,  and  12 -wheel  configurations.  A 
better  performance  dependent  on  a  larger  number  of  wheels  may  be  explained  in 
part  by  the  interior  soil  confinement  afforded  by  a  larger  number  of  perimeter 
wheels.  An  additional  advantage  may  be  that  the  partial  stress  reversals 
imposed  within  the  elastic  domain  are  actually  of  benefit  to  the  performance 
of  the  pavement  structure.  The  decrease  in  the  rate  of  thickness  increase 
with  repetitions  reflects  the  past  and  recent  testing  experience  in  protect¬ 
ing  supporting  layers  from  shear  deformation.  Figure  37  shows  a  plot  similar 
to  that  in  figure  35,  with  a  logarithm  extension  of  the  test  pattern  shown  to 
encompass  different  numbers  of  wheels. 

The  design  curves  to  be  used  to  obtain  have  been  developed 
using  the  load  on  one  landing-gear  assembly  of  an  aircraft.  In  actual 
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Table  8 

Conversion  of  MMHGL  Teat  Section  Coverages 
to  Actual  Airfield  Facility  Passes 


Identification 

Assembly 

A 

3a  in. 

Failure 

Coverages 

«i 

Factor 

p/c 

Factor 

Failure 

Passes 

Lane  1,  item  1 

12  wheel 

285 

8 

0.340 

0.66 

5 

item  2 

12  wheel 

285 

104 

0.550 

0.66 

69 

item  3 

12  wheel 

285 

1500 

0.661 

0.66 

990 

item  4 

12  wheel 

285 

1500 

0.661 

0.66 

990 

item  5 

12  wheel 

285 

3850* 

0.692 

0.66 

2541 

Lane  2,  item  1 

Single  wheel 

285 

6 

0.360 

2.9 

17 

item  2 

Single  wheel 

285 

200 

0.670 

2.9 

580 

Lane  3A,  item  1 

Single  wheel 

285 

120 

0.625 

2.9 

348 

item  2 

Single  wheel 

285 

450* 

0.740 

2.9 

1305 

Lane  3B,  item  3 

Twin  tandem 

290 

40 

0.510 

1.62 

65 

item  4 

Twin  tandem 

290 

40 

0-510 

1.62 

65 

item  5 

Twin  tandem 

290 

280 

0.640 

1.62 

454 

*  Nonfailure. 


....  HiAs&tg Ujgjpjg 
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Table  3 

Previous  Related  St  idles 


Reference 

No. 

Assexhbly 

A 

sq  in. 

Failure 

Passes 

*1 

Factor 

38 

Single  wheel 

1501 

190 

0.646 

Single  wheel 

1501 

2,l40 

0.924 

Single  wheel 

1501 

13 

0.550 

Single  wheel 

1501 

76 

0.686 

Single  wheel 

1501 

455 

0.704 

Single  wheel 

1501 

1,890 

O.98O 

Single  wheel 

1501 

1,640 

0.959 

Single  wheel 

250 

11,650 

0.813 

Single  wheel 

2ri0 

11,650 

0.893 

41 

Twin  wheels 

330 

5,250 

0.787 

45 

Twin  tandem 

150 

4io 

0.758 

150 

135 

0.475 

150 

2,250 

0.870 

4l 

Twin  tandem 

262 

1,310 

0.693 
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AIRCRAFT  TRAFFIC  VOLUME  FACTOR,  PASSES 

Figure  35.  Load  Repetition  Factor  Versus  Passes  for  MWHGL  Study 
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practice.,  however,  designs  are  to  be  based  upon  the  load  on  all  main-gear 
tires.  Therefore,  the  versus  passes  curves  (figure  37)  show  not  only  the 
landing  gear  type  used  to  develop  the  curves  but  also  the  total  number  of  air¬ 
craft  main-gear  tires  represented  by  the  gear  type.  Use  of  the  criteria, 
therefore,  is  accomplished  by  determining  the  ESWL  and  the  c for  all  main- 
gear  tires.  However,  where  it  is  shown  that  some  combination  or  grouping  of 
tires  other  than  all  main-gear  tires  will  produce  a  greater  thickness  require¬ 
ment  than  the  main  tires,  then  the  other  combination  or  group  will  be  used. 

Use  of  the  criteria  for  aircraft  having  a  number  of  wheels  other  than  shown 
requires  interpolation  between  the  curves.  The  unique  limiting  curve  shown 
from  1  to  100  passes  represents  a  composite  of  the  single,  twin-tandem,  and 
12 -wheel  curves  shown  in  figure  35*  Actually,  the  composite  curve  for  the 
very  low  operational  level  is  used  for  convenience  because  there  is  such  a 
small  difference  in  repetitions  effect  and  because  it  is  difficult  to  differ¬ 
entiate  failure  at  a  low  operational  level. 

d.  Stationary  Requirements 


The  static  requirements,  other  than  due  to  repetitions,  may  be  well 
defined  by  the  two  basic  paramet  ^rs  t/VK  and  CBR/p  .  These  parameters  may 
be  reduced  to  a  single  plotted  curve  that  separates  service -behavior  data  with 
regard  to  failures  and  nonfailures.  Such  a  representation  is  shown  in  fig¬ 
ure  11.  Results  from  the  MWHGL  test  section  and  those  from  earlier  related 
test  sections  are  separated  by  this  curve  in  figure  11,  with  the  failure 
located  above  this  curie  and  nonfailures  below.  In  the  region  of  higher 
CBR/pe  and  lower  tyVA  values,  certain  minimum  requirements  (i.e.  drainage 
and  mix  design  requirements)  are  considered  in  the  total  thickness  require¬ 
ments.  To  incorporate  these  parameters  of  design  and  extreme  cases  of  high 
tire  pressure,  light  loads  and  low  tire  pressure,  or  heavy  loads,  the  curve 
shown  in  figure  38  was  developed.  This  curve  represents  an  average  of  these 
extremes  in  this  upper  region  of  CBR/pe  and  approximately  the  same  as  fig- 

CBR  / CBR  \ 

ure  11  in  the  lower  regions  of  -j—  f-jj-  <0.05  I  .  The  statistical  equation 
of  the  best-fit  cubic  curve  is  6  e 


t  -* 

\[a  /-0.0481  -  1.1562  |log  22.  j 

-  0.6414  | 

bf)2 

,  \ 

-  0.4730 

hxf) 

(39) 
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This  cubic  equation  was  the  best-lit  curve  of  the  n  =  1.5  order  polynomials 
attempted.  The  standard  deviation  and  comparison  of  predicted  t/VA  are 
shown  in  table  10. 

e.  Design  Method 

The  cubic  equation  yielded  the  best  statistical  curve  fit  for  the 
separation  of  failures  and  nonfailures.  The  equation  would  be  used  for  a  com¬ 
puter  solution  for  the  computation  of  thicknesses  of  overlying  layers  required 
to  prevent  shear  deformation  in  supporting  layers.  With  the  cubic  equation 
representing  the  static  condition  and  the  factors  developed  previously 
representing  the  effect  of  rep  *itions,  the  thickness  requirement  equation 
based  on  the  MWHSL  test  section  was  finalized  in  the  following  form: 

t  -*at  |VT  ^-0.0481  -  1.1562  ^log 

-  0.6414  /log  SE.  j  .  0.4730  ^log  2p  (40) 

where 

t  =  total  thickness  required  above  supporting  layer. 

=  load-repetition  factor  which  varies  with  number  of  wheels  on  main 
gear  of  aircraft  considered  and  the  volume  of  aircraft  traffic,  in 
passes,  anticipated  as  shown  in  figure  37. 

A  =  measured  contact  area  of  tire,  sq  in. 

Pe  =  ESWL  or  SWL  tire  pressure,  psi.  For  multiple -wheel  gear,  pe 

=  ESWL/A  where  ESWL  is  determined  by  the  method  shown  in  reference 
33  and  for  single-wheel  gear,  p  =  SWL/A  .  This  is  an  artificial 
tire  pressure  for  multiple-wheeleloads  consistent  with  ust  of  con¬ 
tact  area  of  one  tire  and  has  no  relation  to  actual  tire  inflation 
pressure.  However,  for  single -wheel  loads  this  pressure  is  the 
actual  average  contact  pressure  and  is  nominally  the  same  as  the 
tire  inflation  pressure. 

ESWL  =  ESWL  computed  using  all  main -gear  wheels  of  the  aircraft.  However, 
where  it  can  be  shown  that  some  other  combination  or  grouping  of 
wheels  will  produce  a  greater  thickness  requirement  than  the  main- 
gear  wheels,  the  other  combination  or  Troup  will  be  used. 

2.  RIGID  PAVJEMENT 

a.  Plain  Concrete  Criteria 

(l)  12 -wheel  assembly.  A  comparison  of  pavement  performance  observed 
unier  MWHGL  traffic  and  the  predicted  performance  based  on  current  criteria 
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Table  10 

Comparison  of  Predicted  and  Actual  t/VX  Values  Using 
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has  been  made.  Due  to  the  age  of  the  pavement  at  the  time  of  testing  (approx¬ 
imately  1  year),  the  properties  of  the  pavement  structure  determined  in  the 
after-traffic  testing  program  were  considered  more  applicable  than  the  con¬ 
struction  control,  properties  for  the  comparative  analyst:.  Some  increase  was 
noted  in  the  subgrade  modulus  and  concrete  flexural  strengths  from  the  time 
of  construction  to  the  after-traffic  test  program  completion.  An  estimate  of 
strength  at  any  given  time  could  have  been  formvJLated,  but  due  to  the  scarcity 
of  available  data,  a  linear  relationship  would  have  had  to  be  assumed.  Pre¬ 
vious  experience  has  shown  that  the  relationship  between  strength  change  (sub¬ 
grade  and  concrete)  and  time  is  a  curvilinear  relation  that  tends  to  become 
asymptotic  to  a  given  strength  value  as  time  progresses.  Based  on  this  knowl¬ 
edge,  the  after-traffic  test  measurements  were  used  for  subgrade  and  concrete 
strengths  in  the  analysis  of  the  MWHGL  traffic  test  data.  The  physical  prop¬ 
erties  used  for  this  analysis  are  given  in  table  11. 

The  application  of  the  analysis  of  an  elastic  plate  supported  on 
a  dense  liquid  by  Westergaard  is  inherent  in  the  existing  CE  and  FAA  criteria. 
The  adequacy  of  the  current  criteria  will  be  examined  by  first  assessing  the 
applicability  of  the  basic  algorithm  to  the  problem  at  hand.  Some  comparisons 
of  predicted  and  measured  response  are  shown  in  table  12.  The  comparisons 
were  chosen  merely  to  be  representative  of  the  de^v  collected  under  loading. 

In  table  12,  the  computed  values  are  at  the  same  point  in  the  slab  with  re¬ 
spect  to  the  gear  geometry  as  the  measured  values  but  not  necessarily  at  the 
same  points  from  one  test  item  to  another.  Some  measured  values  had  to  be 
discarded  because  of  the  noise  level  in  the  recording;  the  criterion  used  to 
discard  a  reading  was  a  signal -to -noise  ratio  of  less  than  10:1.  These  data 
show  that  the  measured  deflection  for  the  single-wheel  and  twin-tandem-wheel 
gear  loadings  were  normally  slightly  less  than,  but  in  generally  good  agree¬ 
ment  with,  the  computed  values.  This  is  consistent  with  similar  results  from 
other  test  tracks.  The  deflection  measurements  for  the  12 -wheel  static  load 
were  very  erratic  and  inconsistent;  therefore,  the  measured  values  shown  in 
table  12  are  those  taken  during  the  dynamic  load  (slowly  movir.g  load  cart0) 
instrumentation  tests.  From  a  comparison  of  the  few  reliable  static  load 
measurements  available,  there  is  little  difference  between  those  obtained 
from  the  static  and  dynamic  load  tests.  There  :  s,  in  general,  a  fairly  good 
agreement  of  measured  and  computed  deflections;  ^wever  ,  there  is  a  strong 

3]joad  cart  traveled  at  a  rate  of  speed  of  approximately  3  mph. 
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Table  11 

Summary  of  Physical  Properties  Used  for  Analysis 
of  MtfHGL  Test  Results 


Test 

Item 

Slab 

Concrete 

Thickness 

in. 

Overlay 

Thickness 

in. 

Modulus 
of  Rupture 
psi 

Modulus  of 
Elasticity 
10-6  psi 

Subgradi 

Modulus 

lb/ln.3 

12-Wheel  Assembly 

1 

SW 

10.4 

3-9 

725 

6.6 

140 

SE 

10.0 

4.5 

725 

6.6 

140 

2 

SW 

12.0 

N/A 

800 

7.2 

78 

SE 

12.1 

n/a 

800 

7.2 

78 

3 

SW 

l4.1 

n/a 

700 

6.7 

65 

SE 

13.9 

n/a 

700 

6.7 

65 

4 

SW 

8.1 

5.3 

775 

6.3 

125 

SE 

8.1 

5.6 

775 

6,3 

125 

Twin-Tandem  Assembly 

i 

NW 

9.7 

3.9 

695 

6.5 

169 

NE 

9.6 

3.8 

695 

6.5 

169 

2 

NW 

11.5 

n/a 

700 

6.5 

111 

NE 

11.3 

n/a 

700 

6.5 

111 

3 

NW 

13.6 

n/a 

660 

6.4 

115 

NE 

13.9 

n/a 

660 

6.4 

115 

4 

NW 

7.5 

5.4 

605 

6.1 

128 

NE 

7.4- 

5.6 

605 

6.1 

128 

i 


i 

i 

I 

t 


AJWL-TR-70-113 


Table  \? 

Selected  v.oi%rUpr.«  of  Fiedlcted  and  Measured  ravenent  Hesro&se 


— noKB — —  . 

Single  Wheel  Load  single  Wheel  load  5Vtn-T*ndea  Luad  Twin-Tandcs  Load 

Strain  Inflect  ion  Strain  Deflection  gtnln  Dettertlce.  strain  Deflect  irr. 

.«  in.  in.  a  tn.  In.  t»  In.  _ In.  »  in.  _ in. 


Single  Wheel  Load  single  Wheel  load  ?vln-T*nd««  Lund  Twin-Tandem  Load 

Strain  Inflect  ion  Strain  Deflection  lit  rain  Dettertlce.  strain  Deflection 

.«  in.  in.  a  tn.  In.  a  In.  _ In.  »  in.  _ In. 


Test 

Location 

He  as- 

C*>- 

*»r  as- 

"cm* 

Meas- 

Cost- 

Meas- 

C-na- 

Haas  - " 

Cos- 

Meat- 

Con- 

Meas- 

Co»- 

Meas- 

Co«- 

Uw 

cn  Slob 

ured 

fated 

or  *4 

fUted 

ured 

rated 

ured 

ruted 

ured 

rated 

u-ed 

puted 

ured 

ruted 

ured 

rated 

1 

Longitudinal 

Joint 

3 

45 

0.009 

0 

1 

63 

ft 

0.026 

* 

28 

0 .036 

0.037 

ft 

42 

0.051 

0.055 

Transverse 

Joint 

12 

51 

0.C13 

0.02; 

?? 

77 

* 

0.032 

17 

28 

0.042 

0.027 

38 

42 

0.062 

0.041 

Center 

• 

51 

0.C13 

0.013 

29 

n 

• 

0.019 

21 

32 

•• 

0.038 

34 

48 

•• 

0.057 

2 

Longitudinal 

Joint 

0 

33 

0.009 

■>.012 

4 

50 

0.012 

0.018 

14 

22 

0.029 

0.0J3 

30 

33 

0.040 

0.049 

Transverse 

Joint 

1 

3? 

0.016 

0.016 

i 

56. 

* 

0.023 

0 

17 

* 

0.0.30 

2 

26. 

* 

0.0*0 

Center 

16 

3 9 

0.010 

0.009 

25 

57 

0.013 

0.014 

2? 

21 

•• 

0.029 

32 

32 

0.01Q 

O.OLJ 

i 

Longitudinal 

Joint 

4 

24 

0.009 

0.009 

10 

36 

o.ou 

0.013 

15 

19 

0.017 

0.027 

13 

?9 

0.019 

0.041 

Transverse 

Joint 

* 

26 

c.oofi 

0.012 

11 

4? 

0.010 

0.018 

•  ft 

“■ 

0.021 

0.020 

•  • 

•• 

0.024 

0.041 

Center 

11 

39 

0.007 

0.0 0ft 

12 

44 

0.009 

0.011 

8 

21/ 

0.022 

0.024 

c. 

41 

0.027 

0.030 

4 

longitudinal 

joint 

» 

t 

0.0*1 

0.028 

t 

t 

tt 

— 

♦ 

t 

0.l4(. 

0.074 

t 

t 

tt 

— 

Transverse 

Joint 

t 

1 

— 

♦ 

t 

It 

t 

t 

** 

** 

t 

t 

tt 

Cents r 

t 

t 

0.043 

t 

1 

tt 

— 

t 

t 

t 

0.040 

t 

t 

tt 

— 

15 ,000- lb-P«  r -yhee 1 
-heel -Gear  Load 

22,500-lb-Per-Wheel 
12-Wheel-Cear  load 

30,600-lb-Per-Wheel 
12-Wheel-Gear  Load 

St  sin 

reflection 

Strain 

be flection 

Strain 

Deflection 

a  In . 

in. 

P 

in. 

in. 

**- 

in. 

in. 

*eas- 

cZZ 

JEST" 

Co*- 

Meas- 

Coat* 

Meas- 

Co*- 

Meas- 

Con- 

Meas- 

Com- 

ured 

SUM. 

ared 

rated 

ured 

e&i 

.gal 

njtfid 

ssaa 

ured 

gated 

1 

Longitudinal 

83 

27 

0.030 

0.048 

2 

40 

0.013 

0.072 

41 

60 

•ft 

- 

Joint 

Transverse 

40 

29 

0.068 

0.028 

44 

0.049 

0,041 

ft# 

c- 

•  • 

— 

Jjint 

Center 

65 

59 

0.0?4 

0.036. 

•» 

89 

0.002 

0.054 

62 

98 

0.076 

0.086 

? 

Lrmgttudlnal 

R 

?4 

0.045 

n.ngfl 

36 

0.029 

0.042 

69 

46 

C.048 

0.056 

Joint 

Traneverse 

ft# 

20 

0.C10 

0.046 

15 

30 

0.052 

O.O69 

96 

47 

•• 

-- 

Joint 

Center 

23 

44 

0.033 

0.030 

16 

66 

0.024 

0.044 

36 

48 

0.044 

0.060 

3 

Longitudinal 

117 

14 

0.029 

0.025 

64 

.30 

0.024 

0.037 

50 

37 

0.036 

0.050 

Joint 

Transverse 

#ft 

15 

0.033 

0.040 

ft# 

23 

0.044 

0.060 

86 

4o 

0.049 

0.080 

joint 

Center 

78 

20 

0.055 

0.025 

3* 

20 

0.019 

0.037 

42 

40 

•  ft 

— 

4 

Longitudinal 

t 

t 

0.017 

O.O36 

t 

t 

tt 

— 

t 

t 

ft* 

~ 

Joint 

Transverse 

t 

t 

.* 

— 

t 

t 

tt 

— 

t 

t 

## 

— 

Joint 

Center 

t 

t 

0.123 

0.039 

t 

t 

tt 

— 

t 

t 

ft# 

... 

*  Soiae-to-algnal  ratio  excessive. 

»»  C age  inoperable  • 

*  strain  gages  not  installed  In  Iter.  No  computation*  node, 
♦t  Load  not  applied  to  lte«v  because  of  possible  overload. 

»  Magnetic  tape  saturated. 
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tendency  for  the  measured  deflection  to  be  greater  than  the  computed  deflec¬ 
tion.  This  is  opposite  to  that  which  was  experienced  with  the  single-wheel 
and  twin -tandem-wheel  gear  tests  on  this  and  other  test  tracks.  A  rather  poor 
agreement  was  obtained  between  the  measured  and  predicted  strains;  however, 
the  responses  of  the  strain  gages  were  of  doubtful  validity  considering  the 
erratic  nature  of  the  output  signal  and  the  high  loss  of  gages  due  to  mal¬ 
functioning.  A  few  of  the  strain  gages  responded  well  and  gave  good  results. 

A  comparison  of  the  measurements  from  these  gages  with  computed  values  indi¬ 
cated  that  the  measured  strains  were  generally  less  than  the  computed  strains 
for  the  single -wheel  and  twin-tandem-wheel  gear  loadings;  however,  for  the 
12 -wheel  gear  loadings,  the  reverse  was  generally  true.  Therefore,  the  com¬ 
parisons  of  the  measured  and  computed  deflection  and  strain  values  seem  to 
raise  a  question  regarding  the  applicability  of  the  Westergaard  algorithm  for 
the  multiple -wheel  heavy  gear  loadings.  This  should  not  be  surprising,  since 
the  original  algorithm  was  for  single  wheels,  later  was  extended  for  multiple- 
wheel  gears  (twin  and  twin-tandem),  and  is  now  being  further  extended  to  ac¬ 
count  for  deflections  and  strains  under  8-  and  12 -wheel  gears. 

The  Westergaard  algorithm  yields  the  concrete  stress  and  deflec¬ 
tion  under  static  load  conditions.  Airfield  pavements  are  subjected  to  roll¬ 
ing  loads,  inpact  loads,  and  vibrations,  as  well  as  static  loads.  Inasmuch 
as  these  dynamic  and  repetitive  loadings  may  result  in  loadings  excessive  of 
the  static  load  and  may  result  in  fatigue  of  the  pavement  materials,  these 
loadings  must  be  taken  into  consideration  in  the  design  and  evaluation  pro¬ 
cedures.  The  current  CE  design  method  makes  use  of  a  DF  to  account  for  the 
effects  of  rolling  repetitive -type  loadings.  The  DF  is  the  ratio  of  the  con¬ 
crete  flexural  strength  to  the  computed  flexural  stress,  and  from  the  full- 
scale  accelerated  traffic  testing,  the  relationship  between  the  DF  and  traffic 
volume  has  been  established  for  three  separate  and  distinct  levels  of  slab 
distress:  initial  crack,  shattered  slab,  and  couplet e  failure.  The  curves 
presented  in  Section  III  of  this  report  represent  an  tipper  bound  for  these 
tests. 

Although  each  test  item  was  composed  of  four  slabs,  the  traffic 
area  was  laid  out  so  that  the  100  percent  traffic  was  all  in  the  south  slabs 
and  tangent  and  parallel  to  the  longitudinal  joint  between  the  north  and  south 
slabs.  Thus,  for  the  analysis  of  the  12-wheel  gear  results,  only  the  two 
south  slabs  in  each  test  item  are  considered.  Inasmuch  as  there  are  slight 
variations  in  the  physical  properties  of  the  slabs,  each  slab  in  each  test 
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item  has  been  analyzed  separately*  This  results  in  two  data  points  for  each 
test  item*  For  computational  convenience,  use  was  made  of  the  graphical 
solution  of  the  Westergaard  algorithm  from  reference  51.  Using  the  edge¬ 
loading  influence  charts,  a  relationship  between  N  ,  the  number  of  blocks 
covered  by  the  tire  contact  areas  on  the  influence  charts,  and  I  ,  the 
radii  ,  of  relative  stiffness,  which  was  previously  defined  in  equation  21, 
was  developed  and  is  shown  in  figure  39  for  the  C-5A  gear  configuration.  As 
illustrated  in  figure  39>  only  one  of  the  two  6-wheel  arrays  of  the  12-wheel 
main  gear  was  used  for  the  N  versus  t  relationship,  since  the  other  6- 
vheel  array  does  not  change  the  N  count  from  the  influence  chart.  An  exam¬ 
ple  comparative  analysis  for  item  2  is  shown  below. 


Slab  thickness,  h 
Subgrade  modulus,  k 
Concrete  flexural  strength,  R 
Tire  contact  area,  A 
Load  per  wheel,  P 

Modulus  of  elasticity  of  concrete,  E 
Poisson's  ratio  of  concrete,  \x 


SW  Slab 
12.0  in. 

78  lb/in.3 
800  psi 
285  sq  in. 
30,000  lb 
7.2  X  106  psi 
0.20  (assumed) 


SE  Slab 
12.1  in. 

78  lb/in.3 
800  psi 
285  sq  in. 
30,000  lb 
7.2  x  106  psi 
0.20  (assumed) 


Using  the  above-listed  properties,  the  l  value  is  computed  for  each  slab  and 
the  N  value  for  each  slab  is  determined  from  figure  39*  According  to  refer¬ 
ence  51>  the  flexural  stress  a  is  computed  from  the  following  equation:  . 


a  =  -g  0*1) 

10, 000 n 

where 

q  =  contact  pressure  per  tire,  psi,  obtained  by  dividing  the  wheel 
load,  P  by  the  tire  contact  area,  A 

l  -  radius  of  relative  stiffness,  in. 

N  »  number  of  block'  from  influence  chart  (dimensionless) 

h  -  slab  thickness,  in. 

10^000  ®  ccns^an‘ts  associated  with  the  unit  width  moment  of  inertia  (6)  and 
a  scaling  factor  associated  with  the  graphical  solution  (10,000) 
(dimensionless) 


The  current  criteria  are  based  on  the  flexural  stress  developed  at  a  jointed 
edge  by  a  landing  gtur  oriented  parallel  and  tangent  to  the  jointed  edge.  The 
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load-transfer  device  in  the  joint  is  assumed  to  transfer  a  minimum  of  25  per¬ 
cent  of  ihe  load  from  the  loaded  side  of  the  joint  to  the  unloaded  slab; 
therefore,  this  load  transfer  must  be  reflected  in  equation  4l  above.  This  is 
accomplished  in  the  computation  of  q  as  follows: 


q  = 


°'75  &f|0°° ■■  79  ■»! 


The  flexural  stress  computed  using  equation  hi  yields  the  stress  at  the  edge 
of  the  slab  when  the  N  is  determined  from  the  influence  chart  for  edge  load¬ 
ing  and  the  25  percent  load  transfer  is  reflected  in  the  computation  of  q  as 
sham  above.  The  computed  edge  stress  divided  into  the  concrete  flexural 
strength  gives  the  DF  for  the  test  slab.  The  DF  is  then  plotted  versus  the 
number  of  12-wheel  gear  coverages  required  to  produce  the  initial-crack  and 
shattered-slab  conditions  on  the  current  CE  DFS  versus  coverage  relationships 
presented  in  Section  III  (figures  12-14)  to  determine  the  adequacy  of  the  cur¬ 
rent  CE  design  method.  The  following  tabulation  shows  the  development  of  the 


each  slab  of  each 

test 

item. 

Test 

l 

CTe 

R 

Coverages  for  Failure 

Item 

Slab 

in. 

N 

£Si 

psi 

DF 

Initial 

Shattered-Slab 

1 

SW 

46.2 

680 

635 

725 

1.14 

251 

592* 

SE 

44.9 

727 

694 

725 

1.04 

192 

592 

2 

SW 

6l.O 

520 

637 

800 

1.25 

4496  ' 

5008* 

SE 

61.4 

512 

625 

800 

1.28 

3963 

5008* 

SW 

70.5 

*+55 

540 

700 

1.30 

2208 

5008* 

SE 

69.5 

460 

553 

700 

1.27 

592 

5008* 

4 

SW 

39.0 

812 

893 

775 

0.87 

180 

240 

SE 

39.0 

812 

893 

775 

O.87 

181 

240 

*  Slab  did  not  reach  failure  condition  at  maximum  coverage 
level  applied. 

A  comparison  of  the  predicted  versus  actual  performance  for  each  slab  is  shown 
in  figure  40  for  the  initial  failure  condition  and  in  figure  4l  for  the 
shattered-slab  failure  condition.  An  examination  of  figure  40  shows  that 
there  is  some  scatter  of  data  points  for  the  initial  failure  condition; 
however,  with  the  exception  of  item  4,  the  data  agree  fairly  well  with  the 
relationship  currently  used  by  the  CE  criteria  for  jointed-edge  loading. 

Item  3>  especially  the  SE  slab,  failed  prematurely,  which  is  attributable  to 
the  excessive  pumping  that  occurred  with  the  resultant  loss  of  subgrade 
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support.  Both  slabs  of  item  4  performed  identically,  and,  based  upon  the  data 
available,  there  is  no  explanation  for  the  better-than-predicted  performance. 

It  is  possible  that  the  strengths  of  the  concrete  and  subgrade  were  higher 
than  the  measured  values  indicate;  however,  it  is  more  probable  that  for  such 
a  low  coverage  level  this  magnitude  of  data  scatter  is  not  unreasonable. 

Items  1  and  2  performed  almost  as  predicted,  and  overall  the  performance 
data  indicate  that  use  of  the  CE  DF  versus  coverage  relation  with  the  Wester- 
gaard  algorithm  is  reasonable  for  the  ? 2-wheel  assembly  loading. 

Only  three  slabs  were  trafficked  to  the  shattered- slab  failure 
condition  and  the  comparison  of  the  performance  of  these  slabs  with  current 
CE  DF  versus  coverage  relationship  is  shown  in  figure  4l.  All  three  slabs 
reached  the  shattered-slab  condition  in  the  low-coverage  range  and  at  a 
slightly  lower  coverage  level  than  the  existing  relationship  predicted.  Items 
2  and  3,  which  according  to  the  existing  DF  versus  coverage  relation  should 
have  been  approaching  the  shattered-slab  condition,  had  not  reached  that  con¬ 
dition  at  the  conclusion  of  traffic.  Based  upon  their  condition  at  the  con¬ 
clusion  of  traffic,  it  is  probable  that  they  would  have  exceeded  the  predicted 
values . 

(2)  Twin-tandem  assembly.  The  twin-tandem  assembly  traffic  was  per¬ 
formed  in  the  north  paving  lane  of  the  test  track.  Both  the  Westergaard 
algorithm  and  past  test  track  experience  with  twin-tandem  assemblies  indicate 
that  maximum  stresses  occur  in  the  concrete  when  the  tandem  wheels  are 
oriented  parallel  and  tangent  to  the  longitudinal  joint.  However,  since  the 
12 -wheel  gear  had  trafficked  along  the  longitudinal  joint  in  the  south  lane 
with  some  traffic  across  the  joint  in  the  north  lane,  the  effects  of  this 
12-wheel  traffic  on  the  analysis  of  results  of  tests  with  the  twin-tandem 
assembly  along  the  same  joint  would  have  been  indeterminable.  Therefore,  it 
was  decided  to  apply  the  twin-tandem  assembly  traffic  in  the  center  of  the 
north  paving  lane  to  avoid  any  influence  of  the  prior  12 -wheel  traffic.  This 
meant  that  the  maximum  stress  from  this  traffic  should,  according  to  theory, 
occur  when  the  twin  wheels  were  tangent  to  the  transverse  joints.  In  addi¬ 
tion,  theory  indicates  that  this  orientation  of  the  twin-tandem  assembly 
results  in  stresses  only  slightly  less  than  the  more  critical  orientation  with 
the  tandem  wheels  tangent  to  a  joint.  The  analysis  of  the  twin-tandem  assembly 
traffic  results  was  accomplished  in  the  same  manner  as  for  the  12-wheel  traffic 
results,  i.e.  determination  of  a  DF  for  each  slab  of  each  test  item,  which, 
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along  w5.th  coverages  required  to  produce  the  initial  crack  and  shattered  .slab 
in  the  test  item,  were  compared  with  the  current  CE  DF  versus  coverage  rela¬ 
tionship.  The  physical  properties  of  the  test  pavements  were  again  assumed  to 
be  those  measured  after  the  traffic  program  bad  been  completed  for  the  same 
reasons  as  explained  before.  These  properties  are  summarized  in  table  11. 

Only  test  items  2  and  3  could  be  analyzed  as  plain  rigid  pavements,  since 
items  1  and  4  were  nonrigid  overlay  lavements  throughout  the  twin -tandem  as¬ 
sembly  traffic  testing. 

An  example  of  the  analysis  of  pavement  behavior  under  tv- in¬ 
tandem  assembly  traffic  is  discussed  below.  A  curve  of  N  ,  the  number  of 
blocks  covered  on  the  influence  chart  for  edge  loading,  versus  l  ,  the  radius 
of  relative  stiffness,  was  prepared  for  the  twin-tandem  assembly  oriented 
perpendicular  and  tangent  to  a  jointed  edge  (transverse  joint).  This  rela¬ 
tionship  is  shown  in  figure  40.  Utilizing  this  relationship  and  equation  4l, 
the  flexural  stress  caused  by  the  twin-tandem  assembly  can  be  computed  for 
each  test  item.  A  DF  is  then  determined  which,  when  plotted  versus  the  traffic 
causing  failure  on  the  current  CE  DF  versus  traffic  volume  relationship,  shows 
the  comparison  of  current  criteria  to  actual  performance.  A  comparative  analy¬ 
sis  for  item  2  is  shown  below: 


Slab  thickness,  h 

Subgrade  modulus,  k 

Concrete  flexural  strength,  R 

Modulus  of  elasticity  of  concrete,  E 

Tire  contact  area,  A 

Load  per  wheel,  P 

Poisson’s  ratio  for  concrete,  y, 

0 

Tire  contact  pressure,  ■ 


m  Slab 

11.5  in. 

Ill  lb/in.3 
700  psi 

6.5  x  10^  psi 
207  sq  in. 
41,500  lb 
0.20  (assumed) 
150.4  psi 


m  Slab 

11.3  in. 

Ill  lb/in.3 
700  psi 

6.5  X  106  psi 
207  sq  in. 
41,500  lb 
0.20  (assumed) 
150.4  psi 


Using  the  above  listed  properties,  the  i  value  is  computed  for 
each  slab  and  the  value  for  each  slab  is  determined  from  figure  42.  The 
edge  stress  is  then  computed  using  equation  4l,  again  assuming  that  25  percent 
of  the  load  is  transferred  from  the  loaded  to  the  unloaded  slab  through  the 
aggregate  interlock  of  the  transverse  joint.  A  DF  is  then  computed  for  each 
slab  by  dividing  the  computed  edge  stress  into  the  concrete  strength  and  is 
plotted  versus  the  number  of  coverages  of  the  twin-tandem  assembly  that  pro¬ 
duced  the  initial-crack  and  shattered- slab  conditions  on  the  current  CE 
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DF  versus  coverage  relationsnip.  The  Following  tabulation  shows  the  develop¬ 
ment  of  the  DF  for  each  slab  of  each  test  item: 


Test 

Item 

31  ab 

t 

JjK 

JL 

R 

£Si 

DF 

Coverages  for  Failure 
Initial  Shattered  Slab 

2 

m 

53.0 

3**7 

664 

700 

1.05 

40 

680 

11 E 

52.0 

350 

OOO 

700 

1.03 

150 

680 

3 

m 

59.0 

320 

51*2 

660 

1.22 

260 

680* 

m 

60.0 

310 

520 

660 

1.27 

150 

680* 

*  Slat  did  not  reach  failure  condition  at  maximum  coverage 
level  applied. 


A  comparison  of  the  predicted  versus  actual  performance  for  each  slab  is  shown 
in  figure  43  for  the  ir.itiaJ  failure  condition  and  in  figure  44  for  the 
shattered-slab  failure  condition. 

As  shown  in  figure  43,  item  2  performed  about  as  predicted  for 
the  initial  failure  conditions;  however,  item  3  failed  prematurely.  There  is 
no  evident  explanation  for  this  premature  failure.  Unlike  item  3  under  the 
12-wheel-assemblv  traffic,  which  pumped  excessively  and  undoubtedly  contributed 
to  its  premature  failure,  item  3  under  twin-tandem  traffic  did  not  pump  exces¬ 
sively.  As  reported  in  Volume  II  of  this  report,  there  was  some  pumping  in 
item  3  that  was  not  believed  to  be  of  such  severity  as  to  affect  the  perform¬ 
ance  of  the  slab.  In  fact,  from  all  observations  the  performance  of  item  3 
was  entirely  normal.  This  suggests  that  the  slope  of  the  DF  versus  coverage 
relation  may  be  different  for  the  twin-tandem  assembly  than  for  the  12-wheel 
assembly;  hewever,  to  make  such  a  conclusion  from  the  limited  data  available 
is  not  warranted. 

As  mentioned  earlier,  the  twin-tendem  assembly  traffic  was  ap¬ 
plied  in  the  middle  of  the  north  paving  lane,  and  the  critical  stress  was 
selected  as  being  with  the  gear  tangent  to  the  transverse  joint.  This  meant 
that  the  maximum  stress  would  be  parallel  to  the  transverse  joint  and  the 
resulting  crack  would  be  longitudinal.  A  review  of  the  crack  development  and 
progression  described  in  Volume  II  shows  that  the  initial  crack  in  both  items 
2  and  3  was  in  the  longitudinal  direction  and  originated  at  the  tranverse 
joint.  Thus,  the  crack  development  validates  the  assumption  regarding  location 
of  critical  stress 

Eoth  slabs  of  item  2  had  reached  the  shattered-slab  failure  con¬ 
dition  at  the  conclusion  of  the  twin-tandem  assembly  traffic,  whereas  both 
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Figure  43.  Comparison  of  Actual  Performance  Data  with  Current  CE 
Design  Methodology  for  Twin-Tandem  Assembly  Initial 
Failure  Condition.  Pass  per  Coverage  Ratio  -  3.30 


coverages 


Figure  44.  Comparison  of  Actual  Performance  Data  with  Current  CE 
Design  Methodology  for  Twin-Tandem  Assembly  Shattered 
Slab  Failure  Condition.  Pass  per  Coverage  Ratio  =  3*. 
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slabs  of  item  3  contained  only  one  crack  (longitudinal).  A  ooctpariaon  cf  the 
predicted  with  the  actual  performance  for  the  shattered- slab  failure  condition 
is  shown  in  figure  44.  Both  slabs  of  item  2  reached  the  shavttered-slab  tjondi- 
tion  at  a  coverage  level  3lightly  less  than  the  existing  DF  versus  coverage 
relation  would  have  predicted.  These  two  sl&ba,  with  their  k  value  of  1X1 
Ib/in.  ,  performed  as  the  DF  versus  coverage  relation  would  have  predicted 
for  a  k  value  of  about  65  to  70  lb/in.^.  A $  shown,  both  slabs  of  iteii  3 
according  to  the  DF  versus  coverage  relation  should  have  carried  about  7000  to 
8000  coverages  before  reaching  the  shattered-siab  condition,  and  at  the  con¬ 
clusion  of  traffic,  each  contained  only  one  crack.  It  is  probable  that  both 
slabs  would  have  reached  the  shattered -slab  condition  before  the  predicted 
coverage  level. 

b.  Nonrigid  Overlay  Criteria 

The  MWHGL  test  sections  were  utilized  to  provide  information  as  to  the 
adequacy  of  criteria  for  the  strengthening  of  existing  rigid  pavements  with 
nonrigid  overlay. 

Items  1  and  4  were  overlaid  with  asphaltic  concrete  after  seme  struc¬ 
tural  damage  had  occurred  from  the  12-wneel  assembly  traffic.  This  was  done  c. 
to  simulate  an  existing  pavement  that  hai  received  traffic,  was  beginning  to 
show  distress,  and  was  to  be  strengthened  with  a  nonrigid  overlay.  The  non-¬ 
rigid  overlays  were  designed  using  the  properties  of  the  subgrade  and  concrete 
determined  irom  the  construction  control  tests.  The  design  was  for  4000  to 
4500  coverages  of  the  12-wheel  assembly.  A  condition  factor  C  of  1.0  was 
assigned  to  the  slabs  in  item  1,  since  these  slabs  showed  only  initial  crack¬ 
ing.,  and  a  C  of  0.75  was  assigned  to  the  slabs  in  Item  4  inasmuch  as  each 
slab  was  cracked  into  6  pieces.  The  same  thickness  of  nonrigid  overlay  was 
applied  to  both  the  north  and  south  slabs  of  items  1  and  4,  and,  whOii  evaluat- 
ing  the  perfc -mance  of  the  north  slabs  under  twin-sandem  assembly  traffic, 
condition  factors  C  of  1.0  were  assigned,  since  the  north  slabs  were  not 
trafficked  urior  to  the  overlay  and  did  not  contain  structural  cracking. 

In  retrospect,  the  thickness  of  the  nonrigid  overlay  Vised  was  exces¬ 
sive  for  the  12 -wheel  'ssembly  traffic  that  was  programmed  because  both  the 
modulus  of  subgrad'1  re;  ction  and  concrete  flexural  strength  were  considerably 
higher  than  the  values  selected  for  design  from  the  c  instruction  control 
testing.  As  a  result,  neither  item  1  nor  4  nonrigid  overlay  was  failed  with 
the  programmed  12 -wheel  assembly  traffic.  The  evaluation  of  item  1  nonrigid 
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overlay  was  excessive  For  the  programmed  twin-tandem  assembly  traffic  and  a 
failure  was  not  experienced;  however,  item  *4  ncnrigid  overlay  was  considered 
failed  at  the  conclusion  0**  the  programmed  twin-tandem  assembly  traffic. 

The  nonrigid  overlay  test  items  were  analyzed  by  relating  the  combined 
plain  concrete  and  asphaltic  concrete  thickness  to  an  equivalent  thickness  of 
plain  concrete  pavement  as  follows: 

he-?K  +  S3)  <h2> 

where 

h  =  equivalent  plain  concrete  thickness,  in. 

A  table  of  equivalent  thickness  computations  for  all  applicable  slabs  is  shown 
below. 


Gear 

Test 

Item 

Slab 

% 

in. 

c 

F 

k  , 
lb/ in . i 

t 

in. 

ke 

in. 

12-wheel 

1 

SW 

10.4 

1.0 

0.86 

140 

3.9 

13.9 

SE 

10.0 

1.0 

0.86 

140 

4.5 

13-7 

4 

SW 

8.1 

0.75 

0.87 

125 

5.3 

9.4 

SE 

8.1 

0.75 

0.87 

125 

5.6 

9-6 

Twin-tandem 

X 

m 

9.7 

1.0 

0.69 

169 

3.9 

16.3 

m 

94 

1.0 

0.69 

169 

3.8 

16.1 

4 

m 

7.5 

1*0 

0.71 

128 

5.4 

13.6 

m 

7.4 

1.0 

0. 1 1 

r  128 

5*6 

13.6 

Using  the  equivalent  thickness  h  ,  the  test  items  were  analyzed  in  the  same 
manner  as  for  the  pleOn  concrete  test  items.  The  N  versus  l  values  shown 
in  figures  39  and  4$  are  used  for  the  12-wheel  and  twin- tandem  assembly  traf¬ 
fic,  respectively.  The  '"Lexurai  stress  was  computed  using  equation  4l  and  was 
used  to  obtain  a  DF  by  dividing  it  into  the  concrete  flexural  strength.  The 
following  tabu3*H.on  shows  the  determination  of  the  DF  for  each  test  item: 


Tost 

10^ 

£ 

K 

Cover spec 

Gear 

Item 

Slab 

in . 

u 

t>si 

psi 

DF 

to  Fail 

12-,-chael 

x 

SW 

6.^ 

57.5 

550 

72r 

1.78 

If-Mo* 

SE 

f  .0 

rA'^S 

551 

iv',2 

725 

1.6l 

Hilo* 

- 

■■ 

HK 

6.3 

’4  V  *% 

725- 

737 

77- 

1  •  O’ 

nbif.-w 

ss 

6.3 

3 

715 

720 

775 

1.08 

huie* 

TVin-tander. 

■. .  .i. 

If! 

6.7 

6-1 .  L 

300 

383 

695 

1,82 

680* 

titv 

E. 

ho.  9 

302 

389 

69=. 

i  .73 

-jSCr 

nw 

6.1 

r.7 

330 

M5 

1.17 

-7. . .. 

C.l 

A*.  6 

330 

515 

0", 

1.17 

*  50 

*  Traffic  was-,  halted  before  failure  occurred. 


uo 


A 


a 


-■m 

.  i 


AFWUTR-70-U3 


A  comparison  of  existing  criteria  vith  actual  performance  is  made  by  plotting 
the  above  values  of  DF  and  the  coverages  to  failure  in  the  existing  CE  rela¬ 
tionship  of  DF  and  aircraft  coverages  as  shovn  by  figure  U5 - 

A&  stated  earlier,  the  thicknesses  of  nonrigid  overlay  used  were  ex¬ 
cessive  for  the  l?-wheel -traffic  volume  anticipated  because  of  a  greater-than- 
normal  gain  in  concrete  strength  and  subgrade  modulus.  This  is  evident  for 
item  1  when  the  DF  for  each  slab  is  plotted  on  the  CE  relationship.  Based 
upon  this  relationship,  coverages  of  up  to  100,000  would  have  be*n  required  to 
produce  failure.  Based  upon  the  CE  DF  versus  coverage  relationship,  item  4 
should  have  been  approaching  failure  when  traffic  was  halted.  According  to 
Volume  IX,  item  4  at  the  conclusion  of  traffic  was  approaching  failure,  but 
it  is  probable  that  the  item  would  have  exceeded  the  predicted  coverage  level. 
Item  1  for  the  twin-tandem  traffic  was  also  overdesigned  because  of  the  high 
strength  gain  in  concrete  and  rubgrade  and,  as  shown  In  figure  45,  probably 
would  not  have  failed  at  any  reasonable  coverage  level.  Item  4  was  considered 
failed  at  the  conclusion  of  the  twin-tandem  traffic  and  according  to  the  cur¬ 
rent  CE  relationship  of  DF  versus  coverages  failed  prior  to  its  predicted 
coverage  level. 

A  polypropylene  membrane  was  placed  on  the  surface  of  the  southwest 
slab  of  test  item  1  prior  to  placement  of  the  nonrigid  overlay  to  assess  the 
effects  of  the  membrane  on  the  retardation  of  reflection  cracking.  At  the 
time  of  the  nonrigid  overlay,  the  southwest  slab  had  both  a  longitudinal  and 
transverse  crack  which  subdivided  the  slab  into  four  pieces.  The  transverse 
crack  had  occurred  prior  to  the  beginning  of  traffic,  and  during  the  traffic, 
the  crack  spalled  and  raveled  badly.  Prior  to  placement  of  the  polypropylene, 
the  crack  was  cleaned  of  spalled  concrete  and  then  primed  and  filled  with  as¬ 
phaltic  concrete.  The  membrane  did  seem  to  retard  the  reflection  cracking  to 
some  degree;  however  this  cannot  be  supported  by  data,  for  both  slabs  began  to 
show  reflection  cracks  from  the  joints  in  the  base  pavement  even  prior  to 
traffic  on  the  overlay.  During  traffic,  reflection  cracks  from  both  joints 
and  structured  cracks  in  the  jointing  systems  occurred  with  little  difference 
between  the  sections  where  the  polypropylene  was  and  was  not  used. 

G  •  Performance  of  Keyed  Joints 

The  keyed  longitudinal  Joints  failed  early  in  all  test  items  under  the 
MWHGL  traffic.  The  exact  time  of  failure  is  not  known  for  the  first  evidence 
was  faulting  of  the  longitudinal  joint.  Faulting  was  noted  in  most  test  items 
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Figure  45.  Comparison  of  Actual  Performance  Data  with  Current 
CE  Design  Methodology  for  12-Wheel  and  Twin-Tandem 
Assemblies  Complete  Failure  Condition.  Pass  per 
Coverage  Ratios  for  12-Wheel  and  Twin-Tandem  Assem¬ 
blies  =  1.34  and  3.3O,  respectively. 
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after  only  a  few  hundred  coverages,  and  it  is  highly  probable  that  the  fault¬ 
ing  was  the  result  of  the  Keyed  joint  failure;  thus,  the  keyed  joint  failed 
at  an  even  lower  coverage  level.  Past  experience  indicated  that  a  keyed 
joint  would  yield  satisfactory  performance  if  the  pavement  slabs  were  ade¬ 
quate  for  the  flexural  stresses  generated  by  the  load.  Keyed  joint  design 
is  based  on  the  results  of  model  tests  and  observations  from  previous  test 
tracks  and  operational  pavements  containing  keyed  joints. 

The  failure  of  the  keyed  joint  did  not  appear  to  be  associated  with 
the  flexural  stress  due  to  the  lead.  For  example,  in  item  4,  the  flexural 
stress  was  high,  and  in  item  3,  the  flexural  stress  was  low,  with  items  1 
and  2  falling  between;  however,  insofar  as  can  be  determined,  the  keyed  joint 
failed  at  almost  the  same  coverage  level  in  each  item.  Ordinarily  a  keyed 
joint  is  not  used  in  pavements  less  than  9  in.  thick  because  of  the  dif¬ 
ficulty  in  construction  and  because  of  pa3t.  poor  performance  in  thin  pave¬ 
ments;  however,  the  keyed  joint  was  used  in  item  4  (B  in.)  to  reduce  the 
number  of  test  variables.  Nevertheless,  the  performance  of  the  keyed  joint 
in  the  8-in. -thick  item  was  as  good  as  for  the  thicker  pavements  insofar  as 
could  be  determined.  The  performance  of  the  keyed  joint  in  item  3  may  have 
been  influenced  by  the  pumping  that  occurred;  however,  it  is  the  opinion  of 
the  personnel  involved  that  the  keyed  joint  failed  before  pumping  became 
severe  enough  to  influence  the  joint  performance.  As  a  matter  of  fact,  it  is 
probable  that  the  keyed  joint  failure  accentuated  the  pumping  that  occurred. 

Because  both  the  keys  and  keyways  failed  without  a  significant 
pattern,  the  optimal  key  and  keyway  dimensions  adopted  by  the  CE  are  valid. 
The  optimal  keyed  joint  dimensions  are  those  that  will  produce  an  equal 
number  of  key  and  keyway  failures.  No  purpose  would  be  served  by  construct¬ 
ing  a  large,  very  strong  key  as  the  keyway  will  fail  prematurely,  and  con¬ 
versely,  if  the  keyway  is  overdesigned,  the  key  will  fail  prematurely.  It 
is  possible  that  a  thickened-edge  keyed  joint  would  be  satisfactory.  Trans¬ 
verse  contraction  joints  of  this  aggregate-interlock  type  (weakened  plane) 
seemed  to  pose  no  problems  under  multiple-wheel  leadings.  No  abnormal  be¬ 
havior  or  premature  failures  (except  where  pumping  was  severe)  were  experi¬ 
enced.  In  the  test  track,  the  transverse  contraction  joints  were  always 
normal  to  the  longitudinal  axis  of  the  12-wheel  and  twin-tandem  assembles 
and  therefore  probably  had  to  transfer  less  shear  than  did  the  longitudinal 
joint. 
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d.  Subgrade  Behavior 

Changes  in  subgrade  material  under  rigid  pavement  occur  due  to  en¬ 
vironmental  conditions,  mainly  water  content,  and  to  inherent  soil  properties 
such  as  thixotropy.  In  rigid  pavement  design  and  analysis,  the  subgrade 
materials  are  described  by  the  snbgr&de  modulus  or  k  value.  The  su^grade 
modulus  is,  in  essence,  a  spring  constant  that  represents  the  pressure  re¬ 
quired  to  produce  a  unit  deflection  and  is  determined  by  a  plate  bearing 
test.  Changes  in  the  subgrade  support  value  are  reflected  as  changes  in  the 
subgrade  modulus.  The  plate  bearing  test  consists  of  loading  a  nest  of  bear¬ 
ing  plates  up  to  30  in.  in  diameter  and  observing  the  load-deflection  curve 
of  the  plates  for  various  pressures.  Ihe  k  value  is  defined  as  the  ratio 
of  the  10-psi  load  to  the  deflection  that  occurred  under  the  10-psi  load. 

Some  technique  i?  involved  in  the  interpretation  of  the  load-deflection 
curves  as  the  curves  are  seldom  linear.  When  the  load-deflection  curve  is 
nonlinear,  a  straight  line  is  drawn  through  the  portion  of  the  curve  with 

the  least  curvature  and  corrected  for  seating,  and  a  value  of  deflection  is 

\ 

established  for  the  10-psi  load.  The  determination  of  the  portion  of  least 
curvature  is  a  judgment  decision  and  can  drastically  affect  the  results  of 
the  plate  bearing  test. 

Changes  in  subgrade  moduli  were  noted  on  the  rigid  test  section.  The 
load-deflection  curves  observed  on  the  test  section  were  nonlinear  and  of 
nearly  constant  curvature  (see  Volume  II),  making  interpretation  of  the  data 
somewhat  difficult-.  The  changes  in  subgrade  moduli  seemed  to  follow  the 
drainage  pattern  for  the  site.  A  low  spot  was  formed  under  item  3  and  was 
probably  the  cause  of  the  rather  severe  pumping  that  occurred  under  the  test 
item.  The  pumping  was  undoubtedly  aggravated  by  the  ponding  of  water  in  the 
low  area.  In  the  after -traffic  test  program,  no  major  variations  were  found 
in  water  content,  but  the  similarity  between  the  changes  in  subgrade  moduli 
and  the  drainage  pattern  seem  too  close  to  be  coincidence. 

The  changes  in  subgrade  moduli  from  the  time  of  completion  of  the 
embankment  until  the  posttraffic  testing  were  more  uniform  in  the  north  lane 
than  in  the  south  lane.  The  changes  follow  the  drainage  pattern  for  the 
site;  however,  no  apparent  changes  in  water  content  of  sufficient  magnitude 
to  affect  the  subgrade  moduli  to  the  extent  measured  were  found.  Densifica- 
tion  of  the  subgrade  material  due  to  traffic  was  ruled  out  as  a  possibility, 
as  marked  increases  in  subgrade  densities  were  not  found  and  the  water 
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co.  ‘.eat  was  near  saturation,  making  undrained  compaction  nearly  impossible. 

The  pumping  in  item  3  may  be  associated  with  this  drainage  process. 

The  permanent  deformations  experienced  under  many  applications  of 
the  MWHGL  assembly  were  larger  than  those  encountered  previously  under  more 
conventional  wheel  arrangements  on  other  test  tracks  and  operational  pave¬ 
ments.  Large  permanent  deflections  should  be  anticipated  on  rigid  pavements 
subjected  to  operations  of  MWHGL  aircraft  and  some  loss  of  crown  will  prob¬ 
ably  result,  particularly  on  low-strength  subgrades.  Thus,  if  various  soil 
types  are  present,  some  differential  settlements  may  occur.  No  evidence  of 
deep-seated  shear  failures  was  found.  However,  the  deflection  measurements 
indicated  appreciable  magnitudes  of  both  transient  and  permanent  deflections 
at  a  considerable  depth.  Apparently,  the  multiple -wheel  assemblies  are  quite 
efficient  in  distributing  flexural  stresses  in  the  pavement  slab,  but  require 
more  work  from  the  supporting  materials  than,  for  example,  conventional  twin- 
tandem  assemblies. 

e.  Pavement  Requirements 

(l)  Pavements  for  C-$A  aircraft.  Pavements  to  accommodate  the  C-5A. 
aircraft  are  discussed  in  two  categories:  (a)  existing  pavements  (evalua¬ 
tion)  and  (b)  new  pavements  (design).  The  existing  evaluation  method  for 
determination  of  allowable  loadings  at  various  operational  levels  is  adequate 
insofar  as  the  pavement  thickness  and  physical  properties  of  materials  are 
concerned.  The  major  problem  in  evaluation  lies  with  the  jointing  and  load- 
transfer  systems  of  existing  pavements.  Most  existing  pavements  having  the 
thickness  required  for  operation  of  the  C-5A  aircraft  will  have  been  con¬ 
structed  using  the  keyed  construction  joint  on  part  of  the  pavement  system. 
Based  on  the  results  of  the  MWHGL  traffic  tests,  such  pavements  constructed 
on  low-strength  subgrades  will  experience  rapid  joint  deterioration  with  a 
subsequent  increase  in  edge  stress  and  shortened  pavement  life.  Although  not 
demonstrated  in  the  test  track,  pavements  wit’:  keyed  construction  jcints  on 
medium-  to  high-strength  subgrade3  may  also  experience  rapid  joint  deteriora¬ 
tion  under  the  C-5A  aircraft  traffic.  Existing  pavements  resting  directly  on 
pumpable  foundations  can  be  expected  to  experience  pumping,  depending  on  soil 
type,  moisture  conditions,  drainage,  and  traffic  volume,  which  will  result  in 
loss  of  subgrade  support  to  the  pavement  slabs  and  shortened  pavement  life. 

For  new  construction,  use  of  the  Westergaard  equations  for  compu¬ 
tation  of  stresses  in  the  slab  due  to  edge  loading  (equation  21)  is  considered 
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valid.  The  DF  versus  traffic  volume  (figures  12-14)  relationships  are  con¬ 
sidered  valid  for  determination  of  the  effects  of  load  repetitions  and  en¬ 
vironment.  Therefore,  from  the  results  of  the  MWHSL  tests,  no  changes  in 
the  CE  basic  methods  for  determination  of  required  rigid  pavement  thickness 
are  indicated  for  the  C-5A  aircraft.  Use  of  the  standard  keyed  construction 
joint  is  discouraged,  and  use  of  doweled  or  thickened-edge  construction 
joints  regardless  of  subgrade  strength  is  recommended.  Adequately  designed 
filter  courses  or  stabilized  layers  between  the  rigid  pavement  and  pumpable 
subgrade  materials  should  be  used,  and  the  compaction  requirements  for  the 
foundation  materials  should  be  increased.  Recommendations  concerning  these 
factors  are  presented  in  Section  V. 

(2)  Pavements  for  Boeing  7^7  aircraft.  The  existing  CE  method  of 
evaluating  pavements  for  the  Boeing  747  aircraft  is  considered  adequate  in¬ 
sofar  as  pavement  thickness  is  concerned.  While  traffic  was  applied  with 
only  a  twin-tandem  bogie  of  the  aircraft,  inferences  can  be  drawn  from  these 
and  from  the  12 -wheel -assembly  traffic  that  would  be  applicable  to  the  Boeing 
7^7.  The  keyed  construction  joints  in  existing  pavements  are  suspect,  and 
rapid  joint  deter iorati^:  with  subsequently  increased  pavement  stresses  and 
shortened  pavement  life  should  be  anticiy  ted  under  Boeing  747  aircraft  traf¬ 
fic.  The  above  is  especially  true  of  pavements  constructed  on  low-strength 
subgrades  and,  while  not  demonstrated  by  the  MWHGL  tests,  is  probably  true 
for  pavements  on  medium-  to  high-strength  subgrades.  Pavements  on  pumpable 
foundations  will  probably  experience  some  degree  of  pumping,  depending  on 
soil  type,  moisture  conditions,  and  traffic  volumes,  with  a  long  subgrade 
support  and  shortened  pavement  life. 

For  new  construction  to  serve  the  Boeing  747  aircraft,  the  use 
of  the  Hestergaard  equations  for  computation  of  pavement  stresses  and  the 
existing  CE  relationships  of  DF  versus  traffic  volume  (fig-ires  12-14)  is  con¬ 
sidered  valid  for  development  of  rigid  pavement  thickness  requirements. 
Doweled  or  thickened-edge  construction  joints  in  lieu  of  the  keyed  construc¬ 
tion  joints,  adequately  designed  filter  or  base  courses  or  stabilized  layers 
to  prevent  pumping,  and  increased  compaction  requirements  for  the  foundation 
materials  similar  to  that  discussed  for  the  C-5A  aircraft  are  also  recom¬ 
mended  for  the  Boeing  747  aircraft. 
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SECTION  V 

CONCLUSIONS  AND  RECOMMENDATIONS 

1.  CONCLUSIONS 

a.  Flexible  Pavements 

Conclusions  based  on  the  findings  of  this  study  are  given  below. 

(l)  The  flexible  pavement  design  procedure  described  in  this  report 
is  recommended  for  use  in  the  development  of  C-5A  and  Boeing  7^+7  aircraft 
flexible  pavement  design  and  evaluation  curves.  The  procedure  reflects  a 
reduction  of  thickness  requirements  of  existing  multiple-wheel  criteria, 
especially  significant  in  the  higher  operational  level.  Using  these  con¬ 
cepts  and  procedures,  the  following  behavioral  pattern  for  thicknesses  of 
pavement  layers  to  properly  design  underlying  layers  to  resist  shear  defor¬ 
mation  was  derived: 


(2)  The  pass  concept  is  presented  as  being  representative  of  current 
aircraft  volume  of  traffic  and  related  to  actual  lateral  load  applications  at 
airfield  facilities.  Also,  gross  aircraft  loadings  are  used  to  better  reflect 
the  total  effect  of  the  aircraft  and  to  eliminate  the  necessity  of  assembly 
determination  on  aircraft  sijnilar  to  the  Boeing  7^7  aircraft. 

b.  Rigid  Pavements 

The  following  is  a  list  of  conclusions  based  on  analysis  of  the  data 
collected  during  the  testing  with  the  12-wheel  assembly.  Conclusions  are 
based  on  the  results  of  static  load  and  traffic  tests. 

(1)  The  Westergaard  algorithm  for  determination  of  pavement  deflec¬ 
tion  and  stress  under  static  loading  conditions  appears  to  yield  reasonable 
results  for  the  M7KGL  assemblies;  hovrever,  it  appears  to  become  more  ques¬ 
tionable  as  the  number  and  spacing  of  wheels  increase. 

(2)  The  existing  relationship  for  DF  versus  traffic  volume  (cover¬ 
ages)  developed  by  the  CE  for  edge  loading  is  considered  to  be  valid  for  the 
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initial- crack  failure  condition  but  slightly  unoonsorvative  i* or  the  shattered- 
slab  failure  condition.  Based  on  the  fevr  test  points,  the  DF  versus  coverage 
relation  for  nonrigid  overlays  appears  to  be  reasonable.  '' 

( V]  Keyed  longitudinal  construction  joints  for  pavements  on  low- 
strength  subgrades  are  inadequate  for  the  MWHGL  assemblies. 

(4)  Filter  courses  or  other  positive  protection  devices  are  required 
to  prevent  detrimental  effects  of  pumping  when  rigid  pavements  are  constructed 
for  MiGL  assemblies  on  subgrades  composed  of  fine-grained  material. 

2 .  RECOMMENDATIONS 

a.  Flexible  Pavements 

Based  on  the  findings  of  the  flexible  pavement  MWHGL  study,  the  fol¬ 
lowing  recommendations  are  made: 

(1)  The  procedures  presented  in  this  report  should  be  used  for  the 
development  of  flexible  pavement  design  and  evaluation  criteria,  and  ail 
existing  criteria  should  be  modified  to  reflect  the  results  of  this  study. 

(2)  Further  theoretical  studies  should  be  accomplished  using  the 
MWHGL  traffic  and  instrumentation  data. 

b.  Rigid  Pavements 

Based  on  the  results  of  the  tests  conducted  on  the  rigid  pavement 
portion  of  the  M*fHGL  study,  the  following  recommendations  are  made: 

(1)  For  the  construction  of  new  pavements  that  will  be  subjected 
to  MtfHGL  assemblies,  the  use  of  keyed  construction  joints  is  not  recom¬ 
mended.  Either  thickened- edge  or  doweled  construction  joints  are  recommended. 

(2)  Further  work  to  determine  the  adequacy  of  keyed  construction 
joints  on  medium  and  strong  foundations  should  be  performed. 

(3)  When  pumpable  materials  are  encountered,  the  use  of  a  base 
course  or  stabilized  material  is  recommended,  as  pumping  will  be  quite 
severe  where  traffic  volumes  are  high, 

(4)  Current  evaluation  and  design  criteria  for  rigid  pavement  slab 
thickness  should  not  be  changed. 

(5)  Consideration  should  be  given  to  the  establishment  of  a  maximum 
allowable  permanent  deflection  that  can  be  tolerated  in  a  rigid  pavement 
structure . 

(6)  Consideration  should  be  given  to  increasing  compaction  require¬ 
ments  for  subgrades. 
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(?)  Studies  should  be  undertaken  to  examine  possible  methods  of 
strengthening  keyed  construction  joints  in  existing  airfield  pavements. 

(8)  Studies  should  be  undertaken  to  establish  allowable  soil  stresses 
in  rigid  pavement  structures  by  expanding  the  work  already  done  in  flexible 
pavements  and  in  repetitive  triaxial  soils  testing. 


119 


i 


REFERENCES 


AFWL-TR-70-113 


1.  Peutz,  M.  G.  F.;  Jones,  A.;  van  Kempten,  H.  P.;  layered  Systems  Under 
Normal  Surface  Leads,  Computer  Program,  Ko n i nkli J ke/Shell  Labor  at  orium , 
Amsterdam ,  Hoi land . 

2.  Harr,  M.  E. ;  Foundations  of  Theoretical  Soil  Mechanics,  p  19,  McGraw- 
Hill,  Hew  York,  1555^ 

3.  Michelow,  J.;  analysis  of  Stresses  and  Displacements  in  an  N-Layered 
Elastic  System  Under  a  Load  Uniformly  Distributed  on  a  Circular  Area, 
Computer  Program,  California  Research  Corporation,  Richmond,  Calif., 
September  1963* 

4.  Pletta,  F.5  Engineering  Mechanics,  pp  265-271,  Ronald  Press,  New  York, 
1964. 

5.  Timoshenko,  S.;  Young,  D.  H.;  Advanced  Dynamics,  pp  153-181,  McGraw-Hill, 

New  York,  1948.  ~ 

6.  Boresi,  A.  P.;  Elasticity  in  Engineering  Mechanics,  pp  103-112,  Prentice- 
Hall,  1965. 

7.  Smith,  G.  L.;  Smith,  G.  M.;  Advanced  Dynamics,  pp  117-153,  International 
Textbook  Co..  Scranton,  Pa. ,  i960". 

8.  Kristnamurthy,  N.;  Short  Course  Matrix  Operations  -  Finite  Elements  as 
Applied  to  Structural  Analysis,  Auburn  University,  Auburn,  Ala.. , 

February  1970. 

9.  Wang,  C „  K.;  Matrix  Methods  of  Structural  Analysis,  International  Text- 
-book  Co.,  Scranton,  Pa.,  1966. 

10.  Kristnamurthy,  N.;  Short  Course  in  Finite  Elements,  Auburn  University, 
Auburn,  Ala.,  March  1970. 

11.  Scheid,  Francis;  Theory  and  Problems  of  Numerical  Analysis,  pp  15-21, 
McGraw-Hill,  Hew  York,  1968. 

12.  Banks,  D.  C.;  Palmerton,  J.  B,;  Application  of  Finite  Element  Method  in 
Determining  Stability  of  Crater  Slopes,  Preliminary  Report,  Miscellaneous 
Paper  S-68-3,  U.  S.  Army  Engineer  Waterways  Experiment  Station,  CE, 
Vicksburg,  Miss.,  May  1968. 

13«  Timoshenko,  S.  P.;  Goodier,  J.  N. ;  Theory  of  Elasticity,  pp  235-244, 
McGraw-Hill,  New  York,  1970. 

14.  Duncan,  J.  M. ;  Monismiih,  C.  L.;  Wilson,  E.  L.;  ’’Finite  Element  Analysis 
of  Pavements,"  47th  Annual  Meeting  of  Highway  Research  Board,  January 
1968. 

15.  U.  S.  Naval  Civil  Engineering  Laboratory;  Rational  Pavement  Evaluation- 
Review  of  Present  Technology,  AFWL-TR-69-9*  prepared  for  Air  Force 
Weapons  Laboratory,  Kirtlard  Air  Force  Base,  N.  Mex.,  October  1969. 


120 


AFW  1.-TR-  70- 113 


16 .  W i .1  son ,  E .  L . ;  A  Digital  Computer  Program  for  the  Finite  Element  Analysis 
of  Solids  with  Nonlinear  Material  Properties,  Department  ~  of  Civil  Engi¬ 
neering,  University  of  California,  Berkeley,  Calif July  1 9^5. 

17.  Harr,  M.  E.;  Hosner,  J.  C.;  Theoretical  Study  of  landing  Mat  Behevior, 
Contract  Report  S-69-7,  prepared  by  Purdue  Re*search  Foundation,  Purdue 
University,  for  the  U.  S.  Array  Engineer  Waterways  Experiment  Station, 

CE,  Vicksburg,  Miss.,  June  1970. 

18.  •  Vlasov,  V.  2.;  Leont’ev,  N.  M.;  Beams,  Plates,  and  Shells  on  an  Elastic 

Foundation,  Israel  Program  for  Scientific  Translations ,  Jerusalem ,  19&6 . 

19.  Freudenthal,  A.  M.;  Lorsch,  H.  G.;  "The  Infinite  Elastic  Beam  on  a  Linear 
Viscoelastic  Foundation,"  Proceedings  of  the  ASCE,  Journal  of  Engineer¬ 
ing  Mechanics  Division,  V0.1IB3,  No.  EM  1 ,  January  1957,  pp  1158;  1- 
1158;  22. 

20 .  Barks  dale ,  R .  D . . Elastic  and  Viscoelastic  Analysis  of  Layered  Pavement 
Systems,  Ph.D.  Thesis,  Purdue  University,  Lafayette,  Ind.,  198^. 

21.  Chou,  Y.  T.;  "A  General- Theory  of  Stresses  and  Displacements  in  Elastic 
and  Viscoelastic  Layered  Systems,"  Transactions  of  the  Fifteenth  Con¬ 
ference  of  Army  Mathematicians,  ARL'-D  Report  70-1,  pp  141- l68,  June  1969. 

22.  Hudson,  R.  N.5  Pi scontinubua  Orthotropic  Plates  and  Pavement  Slabs, 
Computer  Program  No.  56-9  >  Center  for  Highway  Research,  University  of 
Texas,  Austin,  Tex.,  October  1967- 

23.  Odemark,  N.;  Investigations  ar  to  Elastic  Properties  of  Soils  and  Design 
of  Pavements  According  to  Theory  of  Elasticity.  Statens  Vaginstitute, 
Meddelande  No.  77,  Stockliolm,  Sweden,  19^9* 

24.  Westergaard,  H.  M.,  "A  Problem  of  Elasticity  Suggested  by  a  Problem  in 
Soil  Mechanics"  and  "Soft  Material  Reinforced  by  Numerous  Strong  Hori¬ 
zontal  Sheets,"  Contributions  to  the  Mechanics  of  Solids,  Timoshenko 
Sixtieth  Anniversary  Volume,  Macmillan,  New  York,  1938. 

25.  Hogg,  A.  K.  A.;  "Equilibrium  of  a  Thin  Plate,  Symmetrically  Loaded,  Rest¬ 
ing  on  an  Elastic  Subgrade  of  Infinite  Depth,"  Philosophical  Magazine, 
Series  7,  Vol  25,  March  1938. 

26.  Holl,  D.  L.;  "Thin  Flates  on  Elastic  Foundations,"  Proceedings,  Fifth 
International  Congress  of  Applied  Mechanics,  Cambridge,  Mass”,  193-8. ~ 

27.  Barenberg,  E.  J.;  A  Structural  Design  Classification  of  Pavements  Based 
on  an  Analysis  of  Pavement  Behavior,  Material  Properties,  and  Modes  of 
Failure,  Ph.D  Thesis,  University  of  Illinois,  Urbana,  Ill.,  1965. 

28.  Peattie,  K.  R.;  "A  Fundamental  Approach  to  the  Design  of  Flexible  Pave¬ 
ments,"  Proceedings ,  International  Conference  on  Structural  Design  of 
Asphalt  Pavements,  pp  403-4ll ,  19^2 . 

29.  Vesic,  A.  S.;  Domaschuk,  L.;  Theoretical  Analysis  of  Structural  Behavior 
of  Road  Test  Flexible  Pavements,  Georgia  Institute  of  Technology,  Atlanta 
Ga.  19 6k,  p 29, 


AFWL-TR-70-113 


30.  Brandt,  E.  W.;  Concentration  Factors  for  Stress  Distribution  in  Soil.*>, 
Sols  Soils  No.  21,  pp  11-20,  I96B. 

'jl.  Shockley,  W.  G. Theory  and  Experimentation  Paper,  unpublished,  U.  S. 

Array  Engineer  Waterways  Experiment  Station,  CE,  Vicksburg,  Miss. 

*  32.  Wesbergaard,  H.  M. ;  "hew  Formulas  for  Stress  in  Concrete  Pavements  of 
■  Airfields ASCE  Transactions,  Paper  Ho.  2340,  May  194?. 

-.33.  "Development  of  GBR  Flexible  Pavement  Design  Method  for  Airfields,"  ‘ 
Transactions  of  an  ASCE  Symposium,  Vol  119,  PP  453-589,  1950. 

343  Department-  of  Defense;  Military  Standard  for  Test  Methods  for  Pavement, 
Subgrade,-  Subbase,  and  Base-Course  "Materials ,  MIL-STD-62LA,  Washington, 

D.  C.  j  December  - 19^43  - 

35.  Porter,  0.  J.;  "Development  of  Original  Method  of  Highway  Design," 

ASCE  Transactions ,  Vol  115,  P  46l,  1950. 

36.  Jurgenson,  Leo; -"Application  of  Elastic  Theory  and  Plasticity  of  Founda¬ 
tion  Problems," -Journal,  Boston  Society  of  Civil  Engineers,  July  1934, 

-  p  24P.  •  ‘  -  . 

-37.  Middlebrooks ,  T.  A.;  Bertram,  G.  E.;  "Adaptation  to  the  Design  of  Air¬ 
field  Pavements,"  ASCE  Transactions,  Vol  115,  P  468,  1950. 

i 

33.  U.  S.  Array  Engineer  Waterways  Experiment  Station,  CE;  Collection  of 

Letter  Reports  on  Flexible  Pavement  Design  Curves,  Miscellaneous  Paper 
No .  4-6l ,  Vicksburg ,  Miss.,  June  1951. 

39 ^  0.  J.  Porter  and  Co.;  Accelerated  Traffic  Tests  at  Stockton  Airfield, 

Stockton,  Calif.,  Stockton  Test  No.  2,  7  Vols,  prepared  for  U.  S.  Engi¬ 
neer  Office,  Sacramento,  Calif.,  May  1948.' 

40.  U.  S.  Army  Engineer  Waterways  Experiment  Station,  CE;  Rigid  Pavement 
Tests,  Marietta,  Georgia,  Construction,  Testing,  and  Data  Report, 
Vicksburg,  Miss.,  August  1Q4^. 

41.  U.  S.  Army  Engineer  Waterways  Experiment  Station,  CE;  Effects  of  Traffic 
with  Small  High-Pressure  Tires  on  Asphalt  Pavements,  Technical  Memorandum 
No.  3-314,  Vicksburg,  Miss.,  June  1950. 

42.  U.  S.  Army  Engineer  Waterways  Experiment  Station,  CE;  Design  of  Flexible 
Airfield  Pavements  for  Multiple-Wheel  Landing  Gear  Assemblies,  Report 
No.  1,  Technical  Memorandum  No.  3-3^9>  Vicksburg,  Miss.,  September  1952. 

43.  U.  S.  Army  Engineer  Waterways  Experiment  Station,  CE;  Design  of  Uppex 
Base  Courses  for  Hi #1- Pres sure  Tires,  Report  No.  1,  Technical  Memorandum 
No.  3-373,  Vicksburg,  Miss.,  December  1953. 

44.  u.  S.  Army  Engineer  Waterways  Experiment  Station,  CE;  Proof -Test  Section, 
Columbus  Air  Force  Base,  Technical  Report  No.  3-4-90,  Vicksburg,  Miss., 
December  1958 . 


122 


4FWL-TR- 70-113 


*' 

1 


i 

I 


I 

| 


§ 


* 


i 


4r>.  U.  S.  Army  Engineer  Waterways  Experiment  Station,  CE;  Investigation  of  the 
Design  and  Conttol  of  Asphalt  Paving  Mixtures,  Technical  Memorandum  No. 
3-254,  Vol  3,  App  D,  Vicksburg,  Miss. ,  May  1948. 

46.  U.  S.  Army  Engineer  Waterways  Experiment  Station,  CE;  Investigation  of 
Effects  oi  Traffic  with  High-Pressure  Tires  on  Asphalt  Pavements,  Tech¬ 
nical  Memorandum  No.  3-312,  Vicksburg,  Miss.,  May  1950. 

4?.  Ahlvin,  R.  G. ;  "Consolidated.  CBR  Criteria,’1  Journal  of  Soil  Mechanics 
and  Foundations  Division,  ASCE,  Vol  84,  No.  SM4,  paper  1825,  pp  1825"; 
1-1825;  lb,  October  195tT 

48.  U.  S.  Army  Engineer.  Waterways  Experiment  Station,  CE;  Developing  a  Set 
of  CBR  Design  Curves,  Instruction  Report  No.  4,  Vicksburg,  Miss., 

November  1959* 

49*  U.  S.  Army  Engineer  Waterways  Experiment  Station,  CE;  Compaction  Re¬ 
quirements  for  Soil  Components  of  Flexible  Airfield  Pavements,  Technical 
Report  No.  3-529>  Vicksburg,  Miss.,  November  1959* 

50.  Headquarters,  Department  of  the  Array;  Standard  Practices  for  Bituminous 
Pavements,  Technical  Manua^  No.  5-822-8,  in  process  of  publication, 
Washington,  D.  C. 

51*  Pickett,  G.;  Ray,  G.  K. ;  "influence  Charts  for  Concrete  Pavements," 

ASCE  Transactions,  Paper  No.  2425,  April  1950. 

52.  Cooksey,  D.  L.;  Ladd,  D.  M. :  Pavement  Design  for  Various  Levels  of 

Traffic  Volume,  AFWL-TR- 70-133,  Prepared  by  U.  S.  Army  Engineer  Water¬ 
ways  Experiment  Station,  CE,  for  U.  S.  Air  Force  Weapons  Laboratory, 
Kirtland  Air  Force  Base,  N.'Mex.,  July  1970. 


53.  Bryant,  E.  C.;  Statistical  Analysis,  pp  113-135  and  198-224.  McGraw-Hill, 
New  York,  i960. 

54*.  Snedeior,  G.  W.;  Cochran,  W.  G.;  Statistical  Methods,  Iowa  State  Univer¬ 
sity  Press,  1967. 

55*  Brown,  D.  N.;  Thompson,  0.  0.;  Lateral  Distribution  of  Aircraft 

Traffic,  Miscellaneous  Paper  in  preparation,  U.  S.  Army  Engineer  Water¬ 
ways  Experiment  Station,  CE,  Vicksburg,  Miss. 


123/124 


Security  Classification 


DOCUMENT  CONTROL  DATA  .  R  &  D 

(Security  classification  of  ifle,  body  o  '  nostrmrt  and  indosing  annotation  must  t  ‘ntered  wb  ->n  tba  overall  report  ta  c fasattiad) 


»  ORIGINATING  ACTIVITY  ( Corpora t a  author-.  l2«.  REPORT  SEC'RITV  Cl  A SSI  FI CATION 

U.  S.  Army  Engineer  Waterways  Experiment  Station  I  UNCLASSIFIED 
Vicksburg,  Mississippi  39181 


J  REPORT  TITLE 

MULTIPLE-WHEEL  HEAVY  GEAR  LOAD  PAVEMENT  TESTS:  Volume  IV,  Analysis  of  Behavior  Under 
Traffic 


4  DESCRIPTIVE  notes  (Typa  of  raport  and  inclusive  dates) 

1  January  1968  through  1  August  1971 


•  •  AuTHCRli)  (First  name,  middle  initial,  last  name) 

George  M.  Hammitt.  II 
Ronald  L.  Hutchinson 
John  L.  Rice 


•  ■  REPORT  DATE 

November  1971 

ta.  contract  or  grant  no. 

MI PR  68-7 

b.  PROJECT  NO.  3224 

C. 

d. 

?».  total  no.  of  pages 

134 


9a.  ORIGINATOR'S  REPORT  NUMBER'S} 

AFWL-TR-70-113 ,  Vol  IV 


#6.  other  REPORT  NOtl)  (Any  other  numbers  that  may  be  assigned 
this  report; 


10.  DISTRIBUTION  STATEMENT 

Distribution  limited  to  U  S  Government  agencies  only  because  of  test  and  evaluation 
(1  Nov  1971) .  Other  requests  for  this  document  must  be  referred  to  AFWL  (DEZ) , 
Kirtland  AFB,  NM  87117.  _  _  _ 


12.  SPONSORING  MILITARY  ACTIVITY 

AFWL  (DEZ) 

Kirtland  AFB,  NM  87117 


Ass  ract  'Distribution  Limitation  Statement  B) 

'Flexible  and  rigid  pavement  test  sections  were  constructed  and  tested  to  obtain  data 
to  validate  present  criteria,  to  establish  modifications  to  present  criteria,  or  to 
develop  new  criteria  for  the  evaluation  and  design  of  airfield  pavements  to  be  sub¬ 
jected  to  multiple-wheel  heavy  gear  loads  (MWHGL) .  The  basic  CBR  design  method  was 
expanded  and  modified  to  obtain  a  method  of  design  for  flexible  pavements  subject  to 
traffic  by  MWHGL.  The  recommended  method  reflects  a  reduction  of  existing  d.  S.  Army 
Corps  of  Engineers  (CE)  thickness  requirements  that  is  especially  significant  for 
multiple-wheel  assemblies  in  the  higher  operational  level.  Current  CE  evaluation  and 
design  methods  for  rigid  pavements  are  based  on  stress  in  the  concrete  pavement  as 
calculated  from  the  Westergaard  analysis ^(extrapolations  to  the  existing  criteria  were 
found  to  be  valid  for  MWHGL  assemblies  insofar  as  pavement  thicknesses  were  concerned. 
Results  of  traffic  testing  of  the  rigid  pavement  indicate  that  current  jointing  recom¬ 
mendations  allowing  keyed  construction  joints  may  be  unconservative  when  MWHGL  assem¬ 
blies  will  be  trafficking  a  pavement  resting  on  a  low-strength  subgrade. 


,1473 


Security  ClLSSificution 


